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Abstract 

We review recent progress on the instabilities of black strings and branes both for pure 
Einstein gravity as well as supergravity theories which are relevant for string theory. We 
focus mainly on Gregory-Laflamme instabilities. In the first part of the review we provide 
a detailed discussion of the classical gravitational instability of the neutral uniform black 
string in higher dimensional gravity. The uniform black string is part of a larger phase dia- 
gram of Kaluza-Klein black holes which will be discussed thoroughly. This phase diagram 
exhibits many interesting features including new phases, non-uniqueness and horizon- 
topology changing transitions. In the second part, we turn to charged black branes in 
supergravity and show how the Gregory-Laflamme instability of the neutral black string 
implies via a boost /U-duality map similar instabilities for non- and near-extremal smeared 
branes in string theory. We also comment on instabilities of D-brane bound states. The 
connection between classical and thermodynamic stability, known as the correlated sta- 
bility conjecture, is also reviewed and illustrated with examples. Finally, we examine 
the holographic implications of the Gregory-Laflamme instability for a number of non- 
gravitational theories including Yang-Mills theories and Little String Theory. 
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1 Introduction 

Gravitational instabilities were already considered around a hundred years ago by Sir 
James Jeans [1] in the context of Newtonian gravity, where it was shown that a static 
homogeneous non-relativistic fluid is unstable against long wavelength gravitational per- 
turbations. Subsequently, with the advent of Einstein's theory of general relativity, the 
perturbative stability of gravitational backgrounds such as black holes became an impor- 
tant issue. More recently, string theory has provided a plethora of higher dimensional black 
brane solutions which are captured at low energies by higher-dimensional (super)gravity 
theories. The question of perturbative stability in these backgrounds was examined by 
Gregory and Laflamme in 1993 in the seminal work [2[|3], where it was shown that neutral 
black strings in more than four dimensions suffer from a long wavelength instability, in 
close analogy to the Jeans instability. The main focus of this review will be the important 
progress that has been achieved over the last few years in understanding various aspects 
of the Gregory-Laflamme (GL) instability. 

The classical stability of black holea^ is related closely to the central question of unique- 
ness in gravity. Black holes in four dimensions are known to be classically stable [H [5] . 
This fits nicely with the fact that in four-dimensional vacuum gravity, a black hole in an 
asymptotically flat space-time is uniquely specified by the ADM mass and angular mo- 
mentum measured at infinity [U [8] [9] . Uniqueness theorems \10\ [TT] for D-dimensional 
(D > 4) asymptotically flat space-times state that the Schwarzschild-Tangherlini black 
hole solution [12] is the only static black hole in pure gravity. However, in pure gravity 
there are no uniqueness theorems for non-static black holes with D > 4$ or for black 
holes in space-times with non-flat asymptotics. On the contrary, there are known cases of 
non- uniqueness. An explicit example occurs in five dimensions for stationary solutions in 
asymptotically flat space-time: for a certain range of mass and angular momentum there 
exist both a rotating black hole with S 3 horizon p3] and rotating black rings with S 2 x S 1 
horizons [15J. 

1 With the term "black hole" we refer in this review to any black object irrespective of its horizon 
topology. 

2 See [13] for recent progress in this direction. 
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One may speculate [IB] that a universally valid version of black hole uniqueness con- 
cerns only those solutions that are perturbatively stable. If this is true, then taking into 
account classical instability could restore the lost black hole uniqueness in higher dimen- 
sions. At present, there is no more than circumstantial evidence for this conjecture. 

Beyond the question of uniqueness, instabilities in gravity are important because they 
are related to time-dependent phenomena and provide a glimpse into the full dynamics of 
the theory. Understanding the full time evolution of a condensing instability of a static 
or stationary solution and its end-point is an interesting and important problem that 
provides information about the full non-perturbative structure of the classical equations 
of motion. In this respect, information about the full structure of the static or stationary 
phases of the theory can provide important clues about the time-dependent trajectories 
that interpolate between different phases. In this review, we will see that black holes 
in spaces with compact directions exhibit a rich phase structure and give strong clues 
about time-dependent processes that involve exciting new dynamics like horizon topology 
changing transitions. In particular, for the Gregory-Laflamme instability of the neutral 
black string, a lot of the current research is inspired by the question [17\ I18j of whether 
the Cosmic Censorship Hypothesis is violated [21 [3] in the condensation process of the 
Gregory-Laflamme instability, a process in which the topology of the horizon is conjectured 
to change. 

While the stability properties of black holes in (super)gravity are interesting in their 
own right, their study is also highly relevant for the stability of black holes and branes in 
string theory. At low-energies closed string dynamics is captured by appropriate supergrav- 
ity theories admitting a large variety of interesting solutions with event horizons that have 
received a lot of attention in the past decades. A particularly interesting class of solutions 
in string theory are the D-brane backgrounds. It has been conjectured [El [201 ED EZ] that 
the near-horizon limit of these backgrounds is holographically dual to a non-gravitational 
theory living on the brane (or configuration of branes). This gauge/gravity (or AdS/CFT) 
correspondence implies a deep connection between the thermodynamics of the near-horizon 
brane backgrounds and the thermodynamics of the dual non-gravitational theories. For 
example, it suggests that the gauge theory dual of a classically unstable brane background 
must have a corresponding phase transition. This hints at an interesting connection be- 
tween the classical stability of brane backgrounds and the thermodynamic stability of the 
dual non-gravitational theories. 

In this review, we will discuss instabilities of black holes and branes in pure gravity 
and string theoryo We will focus mainly on Gregory-Laflamme-type instabilities of black 

3 Perturbative instabilities can occur more generally in string theory in vacua without supersymmetry. 
The instabilities are signaled by the presence of a tachyonic mode in the perturbative spectrum with mass 
that is in general of order string scale. In this review we discuss light tachyons that are visible already in the 
supergravity approximation. In recent years much progress has been made in understanding perturbative 
instabilities in string theory (see [231 124] for a status report on open and closed string instabilities) . In 
string theory there is also an interesting relation between perturbative instabilities and the high energy 
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strings and branes that wrap or are smeared on a compact or non-compact direction, but 
along the way we will also comment on a variety of other related cases. The first part of 
this review (Sections 2-4) is a discussion that involves exclusively pure Einstein gravity in 
higher dimensions, and can be seen as separate from the second part which is more related 
to string theory. In this connection, we also refer the reader to the excellent review by Kol 
[27], which contains discussions on several important topics that we only briefly mention. 
In the second part (Sections 5-7), we turn to manifestations of the Gregory-Laflamme 
instability in supergravity solutions that are relevant for string theory and examine the 
implications for the holographically dual gauge theories and little string theory. 

In the pure gravity case, we begin the discussion with an introduction to the Gregory- 
Laflamme instability of the neutral black string and its various manifestations. The black 
string metric is a solution of Einstein gravity in five or more dimensions that has a factor- 
ized form consisting of a Schwarzschild-Tangherlini black hole and an extra flat direction. 
It has been shown [2j [3] that this background suffers from a long wavelength instability, 
known as the Gregory-Laflamme instability, that involves perturbations with an oscillating 
profile along the direction of the string. We review the precise form of this perturbation in- 
cluding the (time-independent) threshold mode where the instability sets in. New features 
appear when the string direction is compactified. In this case, we find that the neutral 
black string has a critical mass, the Gregory-Laflamme mass, below which the solution is 
unstable. 

The existence of the threshold mode at the critical mass strongly suggests that a new 
static non- uniform black string emerges at that point [281 129] . In recent years, non-trivial 
information about the non-uniform black string solution |29[ [3"01 \5T\ [32l [33] has been 
obtained with advanced numerical methods. One remarkable phenomenon of the non- 
uniform string is that it exhibits a critical dimension pD] (see also [3l] ) at which the phase 
transition of the uniform black string into the non-uniform black string changes from first 
order into second order. See also Refs. |31[ I35 [ [36] for the large dimension behavior of the 
threshold mode and |34[ [37] for the connection between the Gregory-Laflamme instability 
and the Landau-Ginzburg theory of phase transitions. 

The unstable neutral black string is, in fact, part of a larger phase diagram that 
consists of different phases of Kaluza-Klein black holes, which are defined to be solutions 
with an event horizon that asymptote to Minkowski space times a circle. This diagram 
includes for instance the phase of black holes localized in the compact direction for which 
much analytical [381 E3 US El H21 03] and numerical [1H 051 06] progress has been made 
in recent years. All the static, uncharged phases can be depicted in a two-dimensional 
phase diagram |47[ |4"H1 |4"9] parameterized by the mass and tension. Despite the fact that 
we have only one compact direction in this case the phase diagram exhibits a very rich 
structure. For example, contrary to a neutral black hole in Minkowski space, where black 
hole uniqueness restricts to just one solution for a given mass, in the case of Kaluza-Klein 

density of states of the theory [251 126] . 
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black holes one finds a continuous family of solutions [50J. An especially fascinating aspect 
of the phase structure of Kaluza-Klein black holes is the occurrence of a phase transition, 
i.e. the occurrence of a point in a continuous line of phases where the topology of the 
event horizon changes \34 \ \51\ \52 \ I33j. More specifically, in this point, called the merger 
point, the localized black hole phase meets the non-uniform black string phase. 

Many of the insights obtained in this simplest case are expected to carry over as we go to 
Kaluza-Klein spaces with higher-dimensional compact spaces [351 E3| > although the degree 
of complexity in these cases will increase substantially. If there are large extra dimensions 
in Nature [531 [Si] the Kaluza-Klein black hole solutions, or generalizations thereof, will 
become relevant for experiments involving microscopic black holes and observations of 
macroscopic black holes. 

A similar structure occurs for black strings and branes in supergravity and string 
theory. We analyze this aspect of the story in the second part of this review focusing 
mostly on the case of non- and near-extremal D-branes smeared on a transverse circle. In 
fact, we will see that this case is directly related to the neutral Kaluza-Klein black holes, 
discussed in the first part of this review, via a boost/U-duality map [351 EH1 ESI [57] . This 
powerful connection enables us to obtain many interesting results for smeared D-branes 
in supergravity. The results are interesting in their own right as statements about various 
phases of string theory, but are also relevant via the gauge/gravity correspondence for the 
phase structure of the dual non-gravitational theories at finite temperature. For instance, 
it is possible to obtain in this way non-trivial predictions (SHI EH ES] about the strong 
coupling dynamics of supersymmetric Yang-Mills theories on compact spaces and little 
string theory. Interesting generalizations involving D-brane bound states |59[ [601 EI] will 
also be discussed briefly in various places in the main text. 

Objects with an event horizon also exhibit thermodynamic properties, in particular 
one can obtain from the horizon area an expression for the entropy using the Bekenstein- 
Hawking formula. It is then possible to study the thermodynamic stability of the black 
holes and branes under consideration. This raises a question: is there any connection be- 
tween thermodynamic stability and classical gravitational stability? Holography suggests 
that such a relation is likely to exist. For a certain class of branes, this connection has 
been formulated more precisely and is known as the Correlated Stability Conjecture (CSC) 
[SH E3 [Ml EH] . The conjecture has proven to be a rather robust statement whose validity 
has been verified in a large set of examples |591 \60\ [66] which we will review. Recently, 
|67j presented a set of counterexamples, which suggest that the conjecture cannot survive 
in full generality in its current form. We will review the key elements of the existing ar- 
guments in favor of the CSC, what may go wrong in these arguments and suggested ways 
to refine the conjecture. 

There are many topics related to this review which will not be covered fully in the main 
text. For instance, besides the Gregory-Laflamme instability, solutions in (super)gravity 
can exhibit a variety of other classical instabilities. We give a list of known instabilities in 
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the concluding section. 

Although there are many exciting developments based on the material presented here, 
we would like to emphasize and summarize now some of the most prominent ones. The 
study of the Gregory-Laflamme instability of neutral black strings has revealed a rich 
phase structure of black objects in Kaluza-Klein spaces, and is a concrete manifestation 
of the large degree of non-uniqueness in higher dimensional gravity, that would be very 
interesting to understand better. Also, the phase diagram of Kaluza-Klein black holes 
has revealed a concrete example of horizon topology changing phase transitions in gravity. 
Further understanding of the properties of the merger point, where two distinct phases 
meet, would provide us with valuable new insights into the above issues. An interesting 
development in this respect is the connection |68|, \69\ [7U[ [7"T] between the merger point 
transition and Choptuik scaling |72j in black hole formation. 

Another important issue concerns the violation of the Cosmic Censorship Hypothesis 
in the decay of the unstable uniform black string, where a naked singularity may be formed 
when the horizon of the black string pinches [2j E] . The results [301 SSI EH [33] on the non- 
uniform black string phase suggest that for certain dimensions, the localized black hole is 
the only solution with higher entropy than that of the uniform black string (for masses 
where the black string is unstable). On the other hand, it was argued in Ref. [T7] that the 
horizon cannot pinch off in finite affine parameter, thus suggesting that it is impossible for 
the black string horizon to pinch off. A way to reconcile these two results is if the horizon 
pinches off in infinite affine parameter (see also [73]). Recently, the numerical analysis of 
[741 l75| [76] indicates that this indeed is the case. If this is correct it would be interesting 
to examine the implications for the Cosmic Censorship Hypothesis. 

A central question regarding the material presented here concerns the microscopic un- 
derstanding of the entropy [771 EH] of the various phases of black holes and branes on 
transverse circles. In particular, it would be very interesting to pursue further a micro- 
scopic explanation [78] of the new phases. 

Another interesting direction concerns the holographic relation between the phase 
structure on the gravity side and the strongly coupled dynamics of the dual gauge theory. 
By now there are many non-trivial examples where purely gravitational phase transitions 
imply via holography phase transitions in strongly coupled gauge theories, which in some 
cases seem to be continuously connected to phase transitions that appear at weak coupling 
(see e.g. Refs. pU MH ED E21 [83] ) . A case relevant to this review is that of Ref. [56] that 
considered two-dimensional supersymmetric Yang-Mills on a circle, which is dual to near- 
extremal smeared DO-branes on a circle. It is remarkable to see how the phase structure 
of black branes in supergravity can be matched qualitatively to the phase structure of 
weakly coupled gauge theories, c.f. Ref. [56] where a new phase was found in the weakly 
coupled gauge theory that is also present in the holographic dual of the strongly coupled 
gauge theory. 

For a more detailed guide to the topics we discuss in this review, we recommend the 
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reader to examine the table of contents and/or read the introductory paragraphs of each 
section. 

2 Classical stability analysis: black strings in pure gravity 

In this section we take a first look at the Gregory-Laflamme (GL) instability of neutral 
black strings in Einstein gravity. We start with a brief reminder of the Jeans instability, 
which is a long wavelength instability already present in Newtonian gravity, and then 
discuss the GL instability of the neutral black string. We will see that there is a nice 
parallel between the two instabilities. Then we go on to consider the GL instability for 
the neutral black string with the string direction compactified along a circle. An important 
feature of the GL instability is a critical wavelength where the GL unstable mode becomes 
a marginal (threshold) mode. This mode signals the existence of a neutral non-uniform 
black string. The GL mode shows up in many different settings of black hole physics, 
all connected to the question of stability. Finally, we summarize the most prominent 
manifestations of the GL instability and conclude with a brief discussion of the instability 
of the Kaluza-Klein bubble. Further developments concerning the GL instability of neutral 
black strings will be discussed in Section [H 

2.1 Long wavelength instabilities in gravity 

Gravitational instability due to long wavelength modes was first discovered around a hun- 
dred years ago by Sir James Jeans pQ in the context of Newtonian gravity. The salient 
features of this instability are as follows. Consider a gravitational perturbation of a sample 
of static matter with constant mass density p = po, constant pressure p = po and zero 
velocity field v = 0. The equations of motion for non-relativistic isotropic hydrodynamics 
with gravity that govern this system are the equation of continuity 

d t p + V -(pv) = , (2.1) 

the Euler equation 

d t v + (v ■ V)v = --Vp+g (2.2) 
P 

and the equations for the gravitational field 

Vx<? = 0, V • g = -4nGp , (2.3) 

where p is the mass density, p the pressure, v the velocity field and g the external gravi- 
tational field with G being the Newton constant. 
Plugging the perturbation 

po + 8p , po + 5p , Sv , g + Sg (2.4) 
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into the equations of motion we ge 

dt5p + p V • 5v = , d t 5v = -—V5p + 5g, V x Sg = , V • Sg = -A-xGSp (2.5) 

Po 

From these equations we deduce that 

8 2 Sp = V 2 5p + 4ttG Po 5p . (2.6) 

Then, demanding that the perturbation should obey the equation of state Sp = v 2 Sp, we 
obtain for a particular Fourier component dp = j[ e tujt - tk - x the dispersion relation 

to 2 = v 2 s k 2 - AnGpo . (2.7) 

It is now evident that for sufficiently long wavelengths 



A > A, S ^ ,2.8) 

there is an unstable mode with lo 2 < 0. This is the unstable mode of the Jeans instability. 
Qualitatively we find that an object becomes unstable to gravitational perturbations if its 
size becomes equal to or larger than its critical Jeans wavelength Aj. Although originally 
derived in Newtonian gravity, this phenomenon appears in a variety of contexts, among 
these cosmological perturbations in the early Universe. As another example, one can view 
black hole formation as a process that occurs in accordance with the Jeans instability. 
Indeed, in this case we have v s ~ c = 1 and by writing the mass as M ~ pr 3 and the 
Schwarzschild radius as 2MG we find that a black hole will form provided r > 1/y/Gp, 
which is precisely the critical Jeans instability wavelength. 

Here we will be interested in the fact that the Jeans instability has another pendant for 
black hole physics, namely the Gregory-Laflamme instability. We will see below that the 
Gregory-Laflamme instability occurs precisely when the size of the system is larger than 
its critical Jeans wavelength. Another classical analogue in which some features of the 
Gregory-Laflamme instability have been successfully observed J3B] is the Rayleigh-Plateau 
instability of long fluid cylinders. This is briefly reviewed in Section [4.41 

2.2 Gregory-Laflamme instability 

Black holes in four dimensions are known to be classically stable [HE]. When the relevance 
of General Relativity in more than four dimensions emerged, it was natural to ask whether 
there exist any higher-dimensional pure gravity solutions with an event horizon exhibiting 
a classical instability. 

Gregory and Laflamme found in 1993 a long wavelength instability for black strings in 
five or more dimensions [21 [3]. The mode responsible for the instability propagates along 



4 In the present analysis we are neglecting the contributions due to the self-gravitation of the sample 
of matter with p = po, p — po and v = 0. One can include these extra contributions in a more careful 
analysis, but the end result will be the same. 
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the direction of the string, and develops an exponentially growing time-dependent part 
when its wavelength becomes sufficiently long. 

The metric for a neutral black string in D = d + 1 space-time dimensions is 

ds 2 = -fdt 2 + f- l dr 2 + r 2 dfl 2 d _ 2 + dz 2 , f = 1 - (2.9) 

where e?f^_ 2 is the metric element of a (d — 2) -dimensional unit sphere. The metric (|2.9J) 
is found by taking the D — 1 dimensional Schwarzschild-Tangherlini static black holqj 
solution [12] and adding a flat z direction, which is the direction parallel to the string. 
The event horizon is located at r = ro and has topology S d ~ 2 x R. 

2.2.1 The Gregory-Laflamme mode 

The Gregory-Laflamme mode is a linear perturbation of the metric (|2.9p . which we will 
denote as 

v + e V' ( 2 - 10 ) 

g^u stands for the components of the unperturbed black string metric (|2.9p . e is a small 
parameter and h^ u is the metric perturbation 

V = »|expf^ + i^P^| (2.11) 

P tt = -ftp, P tr = 7], P rr = f-\, P spheie = r 2 K (2.12) 

where if), rj, x an d k are functions of x = rk/r$ such that the perturbed metric (|2.10p 
solves the Einstein equations of motion. The symbol 5ft denotes the real part. The resulting 
Einstein equations for the perturbation are analyzed in Appendix[X]and reduce to the four 
independent equations (|A.5|) - (|A.8|) . The gauge conditional on h^ v are the tracelessness 
condition (|A.3j) and the transversality conditions (1A.4|) . Combining the gauge conditions 
with the linearized Einstein equations one can derive a single second order differential 
equation for ip 

+ Q d (x)4>'(x) + V d (x)ij(x) = . (2.13) 

The explicit form of the d-dependent rational functions Qd{x) and Vd(x) appears in Ap- 
pendix [A] This equation depends only on k, and the variable x. A thorough analysis 
produces the curve of possible (0, k) values for which (|2. 13|) has a well-behaved solution. 
We plot a sketch of this curve found numerically in [21 [3] in Figure [TJ 

Figure Q] reveals that the curve of possible values of SI and k intersects the fe-axis in 
two places: k = and k = k c , where k c is a critical non-zero wave- number. The fact 

5 The classical stability of these higher-dimensional black hole solutions was addressed in Refs. [841 1851 

6 Various methods and different gauges have been employed to derive the differential equations for the 
GL mode. See Ref. [87] for a nice summary of these, including a new derivation (see also [88]). 



10 



0.4 0.8 1.2 1.6 2.0 2.4 
k 



Figure 1: The Gregory-Laflamme Q(k) curve in d dimensions reprinted from Ref. [2], 
where -D|fi gure = d used in this review. 

that the curve does not intersect the fc-axis at k < follows from the stability of the 
Schwarzschild-Tangherlini black hole (otherwise we would have an unstable mode with 
17 > and k = 0). On the other hand, the presence of an intersection at k = k c is a 
signature of a long wavelength instability: it means that there is an unstable mode for 
any wavelength larger than the critical wavelength 



The critical wave-number k c marks the lower bound of the possible wavelengths for 
which there is an unstable mode. Therefore, we call the mode with k = k c and = the 
threshold mode. It is a time-independent mode of the form 



As we shall see in the following, the threshold mode is important for several reasons: 

(i) It signals the instability of the black string itself. 

(ii) It can be mapped to unstable modes for other gravity solutions, 
(hi) It suggests the existence of a static non-uniform black string. 

Below, we go through each of these different aspects of the threshold mode. The values of 
k c for d = 4, 14, as obtained in [21 [291 [31], are listed in Tabled) 

2.2.2 The Lichnerowitz operator and the threshold mode 

The fact that the mode (|2.1ip - (|2,12p solves the Einstein equations as a perturbation can 
be formulated in the following way. Define the Lichnerowitz operator for a background 




(2.15) 
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d 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


k c 


0.88 


1.24 


1.60 


1.86 


2.08 


2.30 


2.50 


2.69 


2.87 


3.03 


3.18 



Table 1: The critical wave numbers k c for the Gregory-Laflamme instability. 



metric g pu as 

(Ai)^ = -g w g va D K D K - 2R pupa (2.16) 

Here L refers to Lorentzian signature. In the ensuing it will be convenient to use the 
short-hand notation 

AxV = (AlWV • ( 2 - 17 ) 
Then, the statement that the perturbation h pu of g pv satisfies the Einstein equations of 
motion can be stated as the differential operator equation 

A Z V = ( 2 - 18 ) 

Hence what appears in Appendix El is a check that the perturbation (|2.1in - (12.12j) satisfies 
the Lichnerowitz equation (|2.18j) . with (j2.9f) as the background metric. We can write the 
relevant equations more explicitly as 

{A L ) pv p ° (p p(7 e nt / r °e ikz / r A = (2.19) 



In this expression, P^y has no explicit dependence on t and z and Pzll — for all \x. Since 
the background metric f)2.9p has g zz = 1, we deduce that we can rewrite (12.19P as 

(A L ),/ <T (p pCT e n ^°) = -^P^e^o (2.20) 

where A^ is the dimensionally reduced Lichnerowitz operator arising from the metric in 
one dimension less, i.e. with the z-direction excluded. In particular, for the threshold 
mode with O, = and k = k c 

k 2 

^lPc,^u = -^Pc,»u ■ (2.21) 
r o 

This expression has taken the form of an eigenvalue equation for the dimensionally reduced 
Lichnerowitz operator A^. In general, one could look for all possible solutions to the 
eigenvalue equation 

A L U pu = aU pu (2.22) 

thus obtaining the full spectrum of eigenvalues a of the Lichnerowitz operator. It is 
now clear from (|2.2ip that the existence of the threshold mode is equivalent to the exis- 
tence of a negative eigenvalue for the dimensionally reduced Lichnerowitz operator A^. 
Taking the existence of the threshold mode as a signal for the existence of the unstable 
Gregory-Laflamme mode, we deduce that a negative eigenvalue of the dimensionally re- 
duced Lichnerowitz operator is a signal of a long wavelength instability for the solution. 
This statement will be important in Section [6. 11 
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2.2.3 Comparison with Jeans instability 

Is there any relation or analogy between the Gregory-Laflamme instability of neutral black 
strings and the Jeans instability of a massive object with similar physical properties such 
as dimension and mass? For example, does such an object have a Jeans Instability that 
sets in at approximately the same point as the Gregory-Laflamme instability? To see 
such an analogy, we note that the black string in a d + 1 dimensional flat space-time 
has mass per unit length M/L ~ Tq~ 3 /G. From this we deduce that the mass density 
is po ~ M/(Lrp _1 ) ~ I/^Gvq). Thus, for a system with these physical parameters (|2,8|) 
implies that the Jeans instability should set in at wavelengths larger than the critical 
wavelength 

Aj ~ -=3= ~ r (2.23) 
y/Gpo 

where we take the velocity of sound to be the speed of light. So, indeed there is a Jeans 
instability that sets in when A/ro is larger than some number of order one. This is in good 
qualitative agreement with the Gregory-Laflamme instability, which sets in at wavelengths 
larger than the critical wavelength Aql = 27rro/A; c where k c is a number of order one given 
in Table [TJ 

It is interesting to push this correspondence between the Gregory-Laflamme instability 
and the Jeans Instability even further to the large dimension limit d — ► oo, especially 
since the gravitational field behaves in a Newtonian fashion in this limit. Keeping track 
of the d-dimensional factors in the derivation above, one finds for d > 1 that Aj ~ 
tT 3 / 4 r . As we will see in Section 14. 1\ this critical wavelength does not have the right 
asymptotic behavior compared to the one found for the GL mode. Interestingly, however, 
the (classical) Rayleigh-Plateau instability reviewed in Section 14.41 does predict the correct 
asymptotic behavior |36j. 

2.2.4 Gregory-Laflamme instability for the compactified neutral black string 

We have discussed the neutral black string solution in D dimensional Minkowski space. 
However, for many reasons it is more meaningful to consider instead the neutral black 
string in a space-time with a compact direction parallel to the string. For one thing, the 
question of the end point of the black string instability becomes in this setting a more 
well-defined problem. Moreover, the fact that the string is wrapped on a circle allows us to 
stabilize it by making the circle smaller than the critical Gregory-Laflamme wavelength, 
so that the Gregory-Laflamme mode cannot fit into the geometry. In what follows, we 
briefly analyze the neutral black string in a Kaluza-Klein space-time. Later in Section [3] 
we will put this into the more general context of Kaluza-Klein black holes. 

To this end, consider the neutral black string solution (|2.9p with a periodic z coordinate 
whose period we will denote by L. This solution describes a neutral black string in a 
D = d + 1 dimensional Kaluza-Klein space-time M d x S 1 . Indeed, the black string 
solution (12. 9h asymptotes to Ai d x S 1 as r —* oo. Henceforth we shall denote this solution 
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as the (compactified) neutral uniform black string solution, in order to distinguish it from 
another branch of black string solutions that will be discussed later. The mass of the 
compactified neutral uniform black string is 



M 



-(d-2)r, 



d-3 



(2.24) 



WttG^ -' 

It is useful to define a rescaled mass which is dimensionless, using the circumference L: 

16ttG 



fi = 



L d-2 



M . 



Then for the compactified neutral uniform black string we have 



('/ -2)0,,_ 2 (^ 



d-3 



2.25) 



(2.26) 



We see that the Gregory-Laflamme mode (|2.11[) - (|2.12|) cannot obey the correct periodic 
boundary condition on z if L < Agl 5 with Aql given by (12.141) . On the other hand, 
for L > Aqli w e can fit the Gregory-Laflamme mode into the compact direction with 
the frequency and wave number £1 and k in f|2.11[) - f|2.12|) determined by the ratio tq/L. 
Translating this in terms of the mass of the neutral black string (|2,26|) , we deduce that we 
have a critical Gregory-Laflamme mass given by 

MGL = (d-2)fi d _J^J . (2.27) 

For jj, < jUgl the Gregory-Laflamme mode can be fitted into the circle, and the compactified 
neutral uniform black string is unstable. For \x > /Uql, on the other hand, the Gregory- 
Laflamme mode is absent, and the neutral uniform black string is stable. For /i = //ql 
there is a marginal mode which, as we discuss below, signals the emergence of a new branch 
of black string solutions which are non-uniformly distributed along the circle. In Table [2] 
we record the critical Gregory-Laflamme mass /xql that follows from the numerical data 
of Table Q] (obtained in [H [291 HQ ) using (12371) . 



d 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 




3.52 


2.31 


1.74 


1.19 


0.79 


0.55 


0.37 


0.26 


0.18 


0.12 


0.08 



Table 2: The critical masses //ql for the Gregory-Laflamme instability. 



As we mentioned above, for the compactified neutral uniform black string we can 
meaningfully ask about the end-point of the instability of the string for \x < /j-gl- We will 
discuss this question further in Section 13.21 One of the problematic issues regarding the 
question about the end-point of the Gregory-Laflamme instability in the uncompactified 
uniform black string case is the infinite mass of the black string. If the instability ended, 
for example, with a breaking up of the horizon into disconnected pieces, there would have 
to be an infinite process of localized instabilities happening one after another without end. 
This process is regularized in the compactified case and therefore it seems more sensible 
to consider the instability of the neutral black string in the compactified setting. 



14 



2.2.5 Neutral non-uniform black strings 

Let us re-consider now the threshold mode, i.e. the critical Gregory-Laflamme mode with 
$7 = and k = k c . This mode obeys the equation (|2.2ip . or equivalently the equation 

(A L ) M / CT (P^e^o) = . (2.28) 

P C]/J 2,e* fcc2//r ° is a marginal mode that depends explicitly on the compact direction z. There- 
fore, as noticed in [28L 129[, [30] . the mode corresponds to a new static classical solution that 
can be seen as a neutral black string with a horizon that is non-uniform in the z direction. 
In other words, the mode signals the emergence of a new branch of black string solutions 
which are non-uniformly distributed along the circle. We will call this branch of solutions 
non-uniform black strings, as opposed to the uniform black string solutions (j2.9f) . In sec- 
tion [3] we review the properties of the finitely deformed non-uniform black string solution 
that has been found numerically, and the connection with other branches in the full phase 
diagram of solutions. Note that the topology of the horizon is S d ~ 2 x S 1 for both the 
uniform as well as the non-uniform black string (in the compactified case). 

2.2.6 Manifestations of the Gregory-Laflamme mode 

The Gregory-Laflamme mode discussed above shows up in many different set ups in black 
hole physics, all connected to the question of stability. In the following, we describe briefly 
various implications and manifestations of the Gregory-Laflamme mode and the context 
in which they appear: 

• Neutral non-uniform black string branch. This is described above in this section. 
Here the threshold mode corresponds to a static black string solution which is non- 
uniform along the direction of the string. 

• Instability of static Kaluza-Klein bubble. This will be discussed below in Section [2.31 
In this context the threshold mode is double- Wick-rotated into an unstable mode 
for the static Kaluza-Klein bubble. 

• Semi-classical black hole instability in the canonical ensemble. Consider the 5D 
threshold mode (|2.2ip . involving the dimensionally reduced 4D Lichnerowicz op- 
erator A^. By Wick rotating the mode, we obtain an eigenmode u^ u for the 
Lichnerowitz operator Ae of the Euclidean section of the 4D Schwarzschild black 
hole with eigenvalue —k 2 r§ 2 . Using the relation ro = 2G4M, where G4 is four- 
dimensional Newton's constant and M the four-dimensional mass, we deduce that 
Ufj, u obeys the equation 

A B tv = -0.19(G 4 M)- 2 u^ (2.29) 

where we used the value k c = 0.88 from Table [TJ This is precisely the eigenmode for 
the Euclidean section of the 4D Schwarzschild black hole found by Gross, Perry, and 
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Yaffe in 




|. In [89J it was argued as a consequence of the mode ([2,29[) that the 



Euclidean flat space M 3 x S 1 (hot flat space) is semi-classically unstable to nucleation 
of Schwarzschild black holes. In other words, for any non-zero temperature it is 
thermodynamically preferred for a gas of gravitons in Minkowski space to form a 
black hole. It is easily seen using the higher-dimensional threshold modes that this 
can be extended to higher dimensions, i.e. that hot flat space M. d x S 1 is unstable to 
nucleation of {d + l)-dimensional Schwarzschild- Tangherlini black holes. 

• Local thermodynamic instability. In [M] it was shown, following the conjectures of 
|62(, [63] , that a Euclidean negative mode implies a local thermodynamic instability, 
and vice versa. For the Schwarzschild black hole this is clearly the case since the heat 
capacity is negative. If there is in addition a flat non-compact direction one can Wick 
rotate the Euclidean mode to a threshold Gregory-Laflamme mode, implying that 
classical instability for black branes with non-compact directions is in correspondence 
with local thermodynamic instability. We explain and review this topic in Section 



• Gregory-Laflamme modes for charged branes. Finally, one can also map the neutral 
Gregory-Laflamme mode to Gregory-Laflamme modes for charged branes [561 166] . 
Moreover, in |66j it is shown that the Gregory-Laflamme modes in the limit (k, fi) — > 
(0, 0) become marginal modes for extremal branes uniformly smeared along a trans- 
verse direction. We review and discuss these issues in Section [5j 

2.3 Instability of the static Kaluza-Klein bubble 

Kaluza-Klein bubbles were discovered in [90] by Witten. In [90] it was explained that 
the Kaluza-Klein vacuum Ai 4 x S 1 is unstable semi-classically, at least in the absence 
of fundamental fermions. The semi-classical instability of A4 4 x S 1 proceeds through 
the spontaneous creation of expanding Kaluza-Klein bubbles, which are Wick rotated 5D 
Schwarzschild- Tangherlini black hole solutions. The Kaluza-Klein bubble is essentially a 
minimal S 2 somewhere in the space-time, i.e. a "bubble of nothing". In the expanding 
Kaluza-Klein bubble solution the S 2 bubble expands until all of the space-time is gone. 
However, apart from these time-dependent bubble solutions there are also solutions with 
static bubbles, as we now review. 

To construct the static Kaluza-Klein bubble in d + 1 dimensions we take the metric 
(|2.9p of the neutral uniform black string and make a double Wick-rotation in the t and z 
directions. After a convenient relabeling, we obtain the metric 



7 It is possible to extend this matching to the cases where the unstable mode of the Schwarzschild- 
Tangherlini metric has been calculated (see [64] for details). 



6. 




(2.30) 
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We see that there is a minimal (d — 2)-sphere of radius R located at r = R. To avoid a 
conical singularity we need to make z a periodic coordinate with period 

4vri? 

L = _ . (2.31) 

Clearly, the solution asymptotes to M d x S 1 for r — > oo. Notice that the only free 
parameter in the solution is the circumference of the S 1 . 

The static Kaluza-Klein bubble is classically unstable. This can be seen using the 
threshold Gregory-Laflamme mode (12.2ip . After a double Wick-rotation t — > iz, z — ► it of 
the threshold mode, we obtain a mode U pu obeying AU^ V = —k^R~ 2 U^, where A is the 
Lichnerowitz operator for the static Kaluza-Klein bubble (|2.30j) . Using that g# = — 1 in 
(I2.30P we deduce that 

A M / CT (U pa e ±k ^ rlt ) = (2.32) 

Hence, U pa e ±kcR corresponds to an unstable mode for the static Kaluza-Klein bubble 
(|2.30p , implying that the static Kaluza-Klein bubble is classically unstable to small s-wave 
perturbations. 

The classical instability of the static Kaluza-Klein bubble causes the bubble to either 
expand or collapse exponentially fast|f| For five-dimensional Kaluza-Klein space-times, 
there exist initial data [92] for massive bubbles that are initially expanding or collapsing 
[93] . and numerical studies [M] show that there exist massive expanding bubbles and 
furthermore indicate that contracting massive bubbles collapse to a black hole with an 
event horizon. In the latter case, it is noteworthy that the value of \i in (|3,20p is smaller 
than /xql, the Gregory-Laflamme mass, for 4 < d < 9 (as can be seen by comparing 
in (|3.20p to fiQL in Table [3]) . This implies that the static Kaluza-Klein bubble does not 
decay to the uniform black string. It is therefore likely that the Kaluza-Klein bubble in 
that case decays to whatever is the endstate of the uniform black string (see end of Section 

Finally, we note that in Ref. [95] charged bubble solutions were constructed that are 
perturbatively stable, though it was argued that these are non-perturbatively unstable. 
These bubble solutions play a role in the end state of tachyon condensation of charged 
black strings. 



3 Phases of Kaluza-Klein black holes 

In the previous section we focused primarily on the uniform black string and its GL 
instability. In this section we turn to a more general description of the phases of Kaluza- 
Klein black holes. A (d + l)-dimensional Kaluza-Klein black hole will be defined here as 

8 In Ref. [91] the linear stability of static bubble solutions of Einstein-Maxwell theory was examined. A 
unique unstable mode was found and shown to be related, by double analytic continuation, to marginally 
stable stationary modes of perturbed black strings. 
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a pure gravity solution with at least one event horizon that asymptotes to d-dimensional 
Minkowski space times a circle (A4 d x S 1 ) at infinity. We will discuss only static and 
neutral solutions, i.e. solutions without charges and angular momenta. Obviously, the 
uniform black string is an example of a Kaluza-Klein black hole, but many more phases 
are known to exist. In this section, we present their properties as well as the possible 
relation of these phases to the GL instability (see also the shorter review |96j). 

3.1 Physical parameters and definition of the phase diagram 

In this subsection we present a general method of computing the mass fi and relative 
tension n of a Kaluza-Klein black hole, which will be later used to plot each phase in a 
(fj,, n) phase diagram. We review the main features of the split-up of this phase diagram 
into two regions. Finally, we discuss some general results on the thermodynamics of 
Kaluza-Klein black holes. 

3.1.1 Computing the mass and tension 

For any space-time which asymptotes to M d xS 1 we can define the mass M and the tension 
T. These quantities can be used to parameterize the various phases of Kaluza-Klein black 
holes, as we review below. 

Let us define the Cartesian coordinates for (i-dimensional Minkowski space M d as 
t, x 1 , and the radius r = y^iO 3 ^) 2 - m addition we use a coordinate z of period 

L to label the S . Hence, the total space-time dimension is D = d+ 1. In this notation, 
we have for any localized static object the asymptotics 



9tt - - 1 + ^3 > 9zz ^ 1 + ^ , (3.1) 
as r — ► oo. The mass M and tension T are then given by |47[ [4*8] 



The tension in (|3.2p can also be obtained from the general formula for tension in terms 
of the extrinsic curvature [97] analogous to the Hawking-Horowitz mass formula [98J. 
The mass and tension formulas have been generalized to non-extremal and near-extremal 
branes in [57] . Gravitational tension for black branes has also been considered in [991 1100|, 

mn]. 

In what follows, it will be convenient to define the relative tension (also called the 
relative binding energy) as [37J 

M (d-2)c-c, 

This ratio provides a measure of how large the tension (or binding energy) is relative to 
the mass. It is a useful quantity because it is dimensionless and bounded as [17J 

0<n<d-2. (3.4) 
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The upper bound is due to the Strong Energy Condition whereas the lower bound was 
found in [102] 1103] . The upper bound can also be understood physically in a more direct 
way from the fact that we expect gravity to be an attractive force. For a test particle 
at infinity it is easy to see that the gravitational force on the particle is attractive when 
n < d — 2 but repulsive when n > d — 2. 

It is also useful to define a rescaled dimensionless quantity from the mass as 

16vrG N ,„ _s 

» = L d-2 M ■ ( 3 - 5 ) 

The program set forth in [471 US] is to plot all phases of Kaluza-Klein black holes in a 
(fj,, n) diagram. We shall turn to this in the following. 

3.1.2 The split-up of the phase diagram 

According to the present knowledge of the phases of static and neutral Kaluza-Klein black 
holes, the n) phase diagram appears to be divided into two separate regions [50] : 

• The region < n < \/{d — 2) contains solutions without Kaluza-Klein bubbles, and 
the solutions have a local SO(d — 1) symmetry. We review what is known about 
solutions in this part of the phase diagram in Section [3.21 Because of the SO(d— 1) 
symmetry there are only two types of event horizon topologies: S^ -1 for the black 
hole on a cylinder branch and S d ~ 2 x S 1 for the black string. 

• The region l/(d — 2) < n < d — 2 contains solutions with Kaluza-Klein bubbles. We 
review this class of solutions in Section 13.31 This part of the phase diagram, which 
is much more densely populated with solutions compared to the lower part, is the 
subject of [50] . 



3.1.3 Thermodynamics, first law and the Smarr formula 

For a neutral Kaluza-Klein black hole with a single connected horizon, we can find the 
temperature T and entropy S directly from the metric. Together with the mass M and 
tension 7", these quantities obey the Smarr formula |471 148| 

(d-l)TS = {d-2)M - LT (3.6) 

and the first law of thermodynamics [101] |4"5] |4"9"] 

5M = T5S + T5L. (3.7) 

This equation includes a "work" term (analogous to p5V) for variations with respect to 
the size of the circle at infinity. See e.g. Ref. [49J for a proof of the first law based on the 
Smarr formula (|3.6|) using an ansatz (see (|3.1ip below) for black holes/strings on cylinders 
and a class of static Ricci-flat perturbations. See also Ref. [104] for a more general proof 
based on Hamiltonian methods. 
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It is useful to define the rescaled temperature t and entropy s as 

167tGn 



t = LT , s 



L d-i 



-S . 



(3i 



In terms of these quantities, the Smarr formula for Kaluza-Klein black holes and the first 
law of thermodynamics take the form 



(d — l)ts = (d — 2 — n)fj, , Sfjj = t5s . 

Combining the Smarr formula and the first law, we get the relation 

J logs d — 1 
5 log fj, d — 2 — n ' 



(3.9) 



(3.10) 



so that, given a curve n(fi) in the phase diagram, the entire thermodynamics can be 
obtained. 

From (|3.10p it is also possible to derive the Intersection Rule of [17]: for two branches 
that intersect at the same solution, the branch with the lower relative tension has the 
highest entropy for masses below the intersection point, whereas the branch with the 
higher relative tension has the highest entropy for masses above the intersection point. 
An illustration of the Intersection Rule appears in Figure [2j 



Highest entropy 



Highest entropy 




Figure 2: Illustration of the Intersection Rule. 

It is important to note that there are also examples of Kaluza-Klein black hole solutions 
with more than one connected event horizon [491 \5U[ 1105] . The Smarr formula and first 
law of thermodynamics generalize also to these cases. 

3.2 Black holes and strings on cylinders 

In this subsection we review the properties of neutral and static black objects without 
Kaluza-Klein bubbles. These turn out to have local SO(d — 1) symmetry, which implies 
that the solutions fall into two categories, black holes with event horizon topology iS"^ -1 
and black strings with event horizon topology S d ~ 2 x S 1 . First we discuss the ansatz that 
follows from the local SO(d — 1) symmetry. Then we review the uniform black string, 
non- uniform black string and localized black hole phases. These phases are drawn in the 
(//, n) phase diagram for the five- and six-dimensional cases. Finally we consider various 
related topics, including the existence of copies of solutions, the endpoint of the decay of 
the uniform black string and an observation related to the large d behavior. 
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3.2.1 The ansatz 

As we mentioned above, all solutions with < n < l/(d — 2) have, to our present knowl- 
edge, a local SO(d — 1) symmetry. Using this symmetry it has been shown [511 [49] that 
the metric of these solutions can be written in the form 



ds 2 = _ fdt 2 + lJ 



(27T) 



jdR 2 + dy2 + KR 2 d^l_ 2 



r>d-3 



where i?n is a dimensionless parameter, R and v are dimensionless coordinates and the 
metric is determined by the two functions A = A(R, v) and K = K(R, v). The form (|3.1ip 
was originally proposed in Refs. |38|, 1106] as an ansatz for the metric of black holes on 
cylinders. 

The properties of the ansatz (|3.1ip were extensively considered in [38]. It was found 
that the function A = A(R, v) can be written explicitly in terms of the function K(R, v) 
thus reducing the number of free unknown functions to one. The functions A(R, v) and 
K(R, v) are periodic in v with the period 2ir. Note that R = R$ is the location of the 
horizon in (|3.1ip . 

As already stated, all phases without Kaluza-Klein bubbles have, to our present knowl- 
edge, 0<n<l/(d — 2) and can be described by the ansatz (|3.1ip . In what follows we 
review the three known phases in this class and describe their properties. 

3.2.2 Uniform string branch 

The uniform string branch consists of neutral black strings which are translationally in- 
variant along the circle direction. The metric of a uniform string in d + 1 dimensions is 

^ 2 = -(l-^)^ 2 +^-^y) dr 2 + r 2 ^i„ 2 + dz 2 . (3.12) 

As we can see easily from the metric (13. 12f) using A3. ID and (13. 3p . the relative binding 
energy is n = l/(d — 2) for the whole uniform string branch. Hence, the uniform string 
branch is a horizontal line in the (p, n) diagram. The rescaled entropy is given by the 
expression 

eM = C { t l) ^ , (3.13) 
where the constant c[ 9 ^ is defined as 

C[ q) = 4vr(0 (? _i)~^(g - 1)~«=5 . (3.14) 

The horizon topology of the uniform black string is clearly S* 1 x S d ~ 2 , where the S 1 is 
along the circle-direction. 

The most important physical feature of the neutral uniform black string solution is the 
GL instability reviewed in Section 12.21 As stated there, a compactified neutral uniform 



9 The metric (|3.12p corresponds to A(R,v) = K(R,v) — 1 in the ansatz p. lip . 
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black string is classically unstable for [i < /xgl> ^- e - it is unstable for sufficiently small 
masses. See Table [2] for the numerical values of the critical mass /iGL when d < 14. For 
f 1 > A*GL the string is believed to be classically stable. We comment on the endpoint of 
the classical instability of the uniform black string below. 

3.2.3 Non-uniform string branch 

It was realized in [29ji§ that the classical instability of the uniform black string for fj, < 
implies the existence of a marginal mode at /i = ugl, which again suggests the existence 
of a new branch of solutions. Since the new branch of solutions is continuously connected 
to the uniform black string it is expected to have the same horizon topology, at least when 
the deformation away from the uniform black string is sufficiently small. Moreover, the 
solution is expected to be non-uniformly distributed in the circle-direction z since there is 
an explicit dependence in the marginal mode on this direction. 

The new branch, which we call here the non-uniform string branch, has been found 
numerically in [291 [30] [3Tj . Its most prominent features are: 

• The horizon topology is S l x S d ~ 2 with the S 1 being the circle of the Kaluza-Klein 
space-time M d x S 1 . 

• The solutions are non-uniformly distributed along z. 

• The non- uniform black strings have a local SO(d — 1) symmetry and therefore can 
be written in terms of the ansatz (|3.1ip |51[ I49j . 

• The non-uniform string branch meets the uniform string branch at /i = /Uql on the 
line with n = — 2). 

• The branch has relative tension n < l/(d — 2). 

More concretely, considering the non- uniform black string branch for |/i — ^gl| "C 1 one 
obtains 

n(ji) = ~ 7(M " + 0((fi - /i GL ) 2 ) . (3.15) 

Here 7 is a number representing the slope of the curve that describes the non-uniform 
string branch near \x = /xgl- The numerical data for /xql when 4 < d < 14 are given 
in Table [2] and the corresponding results for 7 in Table [3j These data were computed in 

[a m Ea m eh 

The large d behavior of /iGL was examined numerically in [31] and analytically in [35] . 
The latter will be discussed in Section 14. 1[ 
10 See also [28]. 

11 7 in Table [3] has been found in terms of 771 and 02 given in Figure 2 of |31| by the formula 

_ 2(d-l)(d-3) 2 0-2 1 
7 ~ (rf-2) 2 (r;i) 2 /iGL ' 

r]\ and a"2 are also determined in |29l 130] for d = 4, 5. 
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d 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


1 


0.14 


0.17 


0.21 


0.31 


0.47 


0.74 


1.4 


2.8 


7.9 


-40 


-9.2 



Table 3: The constant 7 which, together with //ql of Table [21 determines the behavior of 
the non-uniform branch for \fx — /xgl| <C 1 • 



The qualitative behavior of the non-uniform string branch depends on the sign of 7. 
If 7 is positive, then the branch emerges at the mass \x = /^gl with increasing \x and 
decreasing n. If instead 7 is negative the branch emerges at fj, = //ql with decreasing \x 
and decreasing n. One can use the Intersection Rule of [47] (see Section 13. lh to see that 
if 7 is positive (negative) then the entropy of the uniform string branch for a given mass 
is higher (lower) than the entropy of the non- uniform string branch for that mass. This 
result can also be derived directly from (I3.15P using (|3.10p . We find 

* M ~ 1 - nfJ {d ~T ^ —{^ - mJ 2 + 0((a* - mgl) 3 ) , (3.16) 



s u (/i) 2(d-l)(d-3) 2 mgl 

from which we clearly see the significance of the sign of the parameter 7. In this expression 
S u (fi) (Snu (/•*)) refers to the rescaled entropy of the uniform (non- uniform) black string 
branch. From the data of Table [3] we see that 7 is positive for d < 12 and negative for 
d > 13 [31]. Therefore, as discovered in [31 j . the non-uniform black string branch has a 
qualitatively different behavior for small d and large d, i.e. the system exhibits a critical 
dimension D = 14^1 



Numerical results 

In six dimensions, i.e. for d = 5, Wiseman found in [30J a large body of numerical data 
for the non-uniform string branch. These data are displayed in the (//, n) phase diagram 
on the right side of Figure [3] [37j . The branch starts at \i = /iQL and the data end around 
\i ~ 2.3 fiQi, ~ 5.3. In [3B] numerical evidence has been found that suggests that the non- 
uniform string branch more or less ends where the data of [30] end, i.e. around fi ~ 5.3 
for d = 5, supporting the considerations of [52]. Recently, numerical results extending the 
non-uniform branch into the strongly non-linear regime were also obtained for the case 
of five dimensions, i.e. for d = 4, in Ref. [32], and for the entire range d < 5 < 10 in 
Ref. [33j. All these non- uniform black string solutions have masses above the GL mass 
and less entropy than the corresponding uniform black string. 

Another interesting feature of the non-uniform solutions is the local geometry around 
the "waist" of the most non-uniform solutions. This is cone-like to a fairly good approxi- 
mation [52j[33], lending support to the ideas proposed in [33J [68] (see also Section I3~3]) . As 

12 The first occurrence of a critical dimension in this system was given in [34] . where evidence was given 
that the merger point between the black hole and the string depends on a critical dimension D = 10, in 
such a way that for D < 10 there are local tachyonic modes around the tip of the cone (the conjectured 
local geometry close to the thin "waist" of the string) which are absent for D > 10. 
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we further discuss below, all these data provide important evidence that the non-uniform 
branch has a topology changing transition into the localized black hole branch. 



3.2.4 Localized black hole branch 

On physical grounds, it is natural to expect a branch of neutral black holes in the space- 
time Ai d x S l . One defining feature of these solutions is an event horizon of topology 
S d ~ l . We will call this branch the localized black hole branch in the following, because 
the S^ -1 horizon is localized on the S 1 of the Kaluza-Klein space. 

Neutral black hole solutions in the space-time .M 3 x S 1 were found and studied in 
[T071 EUSJ [HEJ QTD]. However, the study of black holes in the space-time M d x S 1 for 
d > 4 has begun only recently. The complexity of the problem stems from the fact that 
such black holes are not algebraically special and moreover from the fact that the 
solution cannot be found using a Weyl ansatz since the number of Killing vectors is too 
small. 

Analytical results 

Progress towards finding an analytical solution for the localized black hole was made in 
[38] where, as reviewed above, the ansatz (|3.1ip was proposed. Subsequently it was proven 
in |5H |4*9] that the localized black hole can be put in this ansatz. 

In [39] the metric of small black holes, i.e. black holes with mass /jC 1, was found 
analytically using the ansatz (|3.1ip of [38] to first order in fi. Subsequently, an equivalent 
expression for the first order metric was found in Refs. [401 134] using a different method. 
An important feature of the localized black hole solution is the fact that n —* as fi — > 0. 
This means that the black hole solution becomes more and more like a (d + l)-dimensional 
Schwarzschild black hole as the mass goes to zero. For d = 4, the second order correction 
to the metric and thermodynamics have recently been studied in [42J. More generally, the 
second order correction to the thermodynamics was obtained in Ref. [13] for all d using 
an effective field theory formalism in which the structure of the black hole is encoded in 
the coefficients of operators in an effective worldline Lagrangian. 

The first order result of [5S] and second order result of [33] can be summarized by 
giving the first and second order corrections to the relative tension n of the localized black 
hole branch function of u I 




n = 



(d-2)C(d-2) 

2(d-i)n d _ l 



(d - 2)C{d - 2) 
2(d-l)Q d _! 



H +^(M 3 )- 



(3.17) 



Plugging this expression into (|3.10p one can find the leading correction to the thermody- 



13 



Here C(p) = E 



n 



-p 



is the Riemann zeta function. 
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namics as 



Sbh(At) = C 




2(d-2)n d _i 



C(d-2) 



(rf 2 - 6rf + 7) 
2(d- 1) 



( 



2(d-2)fi d _i 



C(d-2) 



M +0(/^ 3 ) 




where C 1 } is defined in (|3.14[) . This constant of integration is fixed by the physical 
requirement that in the limit of vanishing mass we should recover the thermodynamics of 
a Schwarzschild black hole in (d + l)-dimensional Minkowski space. 

Numerical results 

The black hole branch has been studied numerically for d = 4 in [441 146] and for d = 5 
in [45 1 146] . For small fj,, the impressively accurate data of |46j are consistent with the 
analytical results of [391 US H2] . We have displayed the results of [46] for d = 4, 5 in a 
(fj-,n) phase diagram in Figure [3] 

Amazingly, the work of [46j seems to give an answer to the question: "Where does 
the localized black hole branch end?" . Several scenarios have been suggested, see [49] for 
a list. The scenario favored by [46J is the scenario suggested by Kol [34] in which the 
localized black hole branch meets with the non-uniform black string branch in a topology 
changing transition point. This is strongly implied by the (fj,, n) phase diagram for d = 5 
in Figure [3] Further evidence arises by examining the geometry of the two branches near 
the transition point, and also by examining the thermodynamics [46J. 

Hence, it seems reasonable to expect that the localized black hole branch is connected 
with the non-uniform string branch in any dimension. This means that we can go from 
the uniform black string branch to the localized black hole branch through a connected 
series of static classical geometries. The point in which the two branches are conjectured 
to meet is called the merger point. See Section 14.4.31 for more on the critical behavior of 
the two branches near this point. 

3.2.5 Phase diagrams for d = 4 and d = 5 

In Figure [3] we display the (ji,ri) diagram for d = 4 and d = 5, which are the cases where 
the most information about the various branches of black holes and strings on cylinders 
has been accumulated. For d = 4 we have shown the complete non- uniform branch as 
obtained numerically by [32], which emanates at /xql = 3.52 from the uniform branch 
given by the horizontal line n = 1/2. For d = 5, we have shown the complete non-uniform 
branch, as obtained numerically by Wiseman [30], which emanates at /xgl = 2.31 from 
the uniform branch that has n = 1/3. These data were first incorporated into the (/x, Jt) 
diagram in Ref. [47J. For the black hole branch we have plotted the recently obtained 
numerical data of Kudoh and Wiseman [46], both for d = 4 and d = 5. It is evident from 
the figure that this branch has an approximate linear behavior for a fairly large range of 
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\x close to the origin and the numerically obtained slope agrees very well with the analytic 
result (13171) . 

These data strongly suggest that the localized and non-uniform phases meet in a 
topology changing transition point. Another interesting feature to note is the upwards 
bending, both for d = 4 and d = 5 at the end of the non-uniform branch. 



n 




Figure 3: Black hole and string phases for d = 4 (left) and d = 5 (right), drawn in the 
(/i,n) phase diagram. The horizontal (red) line at n = 1/2 and 1/3 respectively is the 
uniform string branch. The (blue) branch emanating from this at the Gregory-Laflamme 
mass is the non-uniform string branch. For d = 4 and d = 5 this branch was obtained 
numerically in [32] and [30]. The (purple) branch starting at the point (fj,,n) = (0,0) is 
the black hole branch which was obtained numerically by Kudoh and Wiseman [46] . For 
both dimensions the results strongly suggest that the black hole and non-uniform black 
string branches meet. 



3.2.6 Copies of solutions 

In [32] it has been argued that one can generate new solutions by copying solutions on 
the circle several times, following an idea of Horowitz [18]. This works for solutions which 
vary along the circle direction (i.e. in the z direction), so it works both for the black hole 
branch and the non-uniform string branch. Let k be a positive integer. Then if we copy 
a solution k times along the circle we get a new solution with the following parameters: 

A = ^3> n = n, t = kt, s = -^. (3.19) 

See Ref. [49] for the corresponding expression of the metric of the copies in the ansatz 
()3.1ip . Using the transformation (|3.19p . one easily sees that the non-uniform and localized 
black hole branches depicted in Figure[3]are copied infinitely many times in the (fi, n) phase 
diagrams. 
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Beyond these copies of solutions, more general multi black-hole configurations with lo- 
calized black hole solutions can be shown to exist as well [105J. These solutions correspond 
to having several localized black holes of different sizes on the transverse circle. 

3.2.7 The endpoint of decay of the uniform black string 

As we mentioned above, the uniform black string is classically unstable for fj, < fiQL- It is 
natural to ask: "What is the endpoint of this classical instability?" . 

The entropy of a small localized black hole is much larger than that of a black string 
with the same mass, i.e. Sbh(/-*) ^ s u(aO for [i -C 1, as can be easily seen by comparing eqs. 
(I3.13P and (|3.18j) . This suggests that a light uniform string will decay to a localized black 
hole. This is the global thermodynamic argument that appeared already in the original 
work [21 [3]. However, one can imagine other possibilities, for example the uniform black 
string could decay to another static geometry, or it could even keep decaying without 
reaching an endpoint. 

For d = 4, 5 we can be more precise about this issue. Since both cases show similar 
behavior, we focus on d = 5 and have displayed in Figure [5] the entropy s versus the mass 
H diagram for the localized black hole, uniform string and non- uniform string branches. 
We see from this figure that in six dimensions the localized black hole has always greater 
entropy than the uniform strings in the mass range where the uniform string is classically 
unstable. This fact together with the absence of the non-uniform string branch in this 
mass range suggest that the unstable uniform black string decays classically to the localized 
black hole branch in six dimensions (d = 5). The same conclusion holds for five dimensions 
(d = 4), and is expected to hold in any dimension less than fourteen (see below). 




o.9 l u/M-GL 

0.5 1 1.5 2 2.5 

Figure 4: Entropy s/s u versus the mass [i/hgl diagram for the uniform string (red), 
non- uniform string (blue) and localized black hole (purple) branches. 

The viewpoint that an unstable uniform string decays to a localized black hole has 
been challenged in [17] . In this work it is shown that in a classical evolution an event 
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horizon cannot change topology, i.e. it cannot pinch off in finite affine parameter (on the 
event horizon). 

However, this does not exclude the possibility that a classically unstable horizon 
pinches off in infinite affine parameter [73]. Indeed, in [75\ I76j the numerical study [74 J of 
the classical decay of a uniform black string in five dimensions was reexamined, suggesting 
that the horizon of the string pinches off in infinite affine parameter. 

Interestingly, the classical decay of the uniform string is quite different for d > 13. 
As we reviewed above, the results of [31] show that for d > 13 the non-uniform string 
has decreasing mass [i for decreasing n, as it emanates from the uniform string at the 
Gregory-Laflamme point at /i = //gl- This means that the entropy of a non-uniform black 
string is higher than the entropy of a uniform string with same mass. Hence, for d > 13 
we have a certain range of masses where the uniform black string is unstable, and where 
a non-uniform black string exists with higher entropy. This suggests that a uniform black 
string in that mass range will decay classically to a non-uniform black string. This type of 
decay can occur in a finite affine parameter, according to [T7], since the horizon topology 
remains fixed during the decay. 

Note that the range of masses for which we have a non-uniform string branch with 
higher entropy is extended by the fact that we have copies of the non-uniform string branch. 
The copies, which have the thermodynamic quantities given by the transformation rule 
(|3.19p . can easily be seen from the Intersection Rule of [47J (see Section [37TJ) to have higher 
entropy than that of a uniform string of the same mass, since they also have decreasing 
\x for decreasing n. Thus, for d > 13, it is even possible that there exists a non-uniform 
black string for any given \x < /xql with higher entropy than that of the uniform black 
string with mass /i. This will occur if the non-uniform string branch extends to masses 
lower than 2 3_ci /iGL- Otherwise, the question of the endpoint of the decay of the uniform 
black string will involve a quite complicated pattern of ranges. 

We will return to the ci-dependence of the phase diagram in Section 14.11 where we 
consider the large d behavior of a number of relevant quantities. 

3.3 Phases with Kaluza-Klein bubbles 

Until now all the solutions we have been discussing lie within the range < n < X/(d — 2). 
But, there are no direct physical restrictions on n preventing it from being in the range 
l/(d — 2) < n < d — 2. Hence, the question arises about the existence of solutions in this 
range. This question was answered in [50] where it was shown that pure gravity solutions 
with Kaluza-Klein bubbles can realize all values of n in the range l/(d — 2) < n < d — 2. 

In this subsection we start by reviewing the place of the static Kaluza-Klein bubble (see 
Section [2.3|) in the (n,n) phase diagram. We then discuss the main properties of bubble- 
black hole sequences, which are phases of Kaluza-Klein black holes involving Kaluza-Klein 
bubbles. In particular, we comment on the thermodynamics of these solutions. Subse- 
quently, we present the five- and six-dimensional phase diagrams as obtained by including 
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the simplest bubble-black hole sequences. Finally, we comment on non-uniqueness in the 
phase diagram. 

3.3.1 Static Kaluza-Klein bubble 

In Section 12.31 we reviewed the static Kaluza-Klein bubble. Since it is a static solution of 
pure gravity that asymptotes to A4 d x S l the static Kaluza-Klein bubble belongs to the 
class of solutions we are interested in. Thus, it is part of our phase diagram, and using 
flEED, (pOjh (GOD, ([33D and ([23l"D we can read off /i and n as 

(7 g \ d— 3 
-^-j , n = d-2. (3.20) 

These expressions give the static Kaluza-Klein bubble as a specific point in the (/i, n) 
phase diagram. This is consistent with the fact that this solution does not have any free 
dimensionless parameters. Notice that n = d — 2 precisely saturates the upper bound on 
n in (13. 4p . In fact, a test particle at infinity will not experience any force from the bubble 
(in the Newtonian limit). 

3.3.2 Bubble-black hole sequences 

We now have a solution at n = d — 2 with no event horizon, and hence no entropy or 
temperature. But so far in this review we have not mentioned any solutions lying in the 
range l/(d — 2) < n < d — 2. It was argued in [50] that such solutions exist and contain 
both Kaluza-Klein bubbles and black hole event horizons with various topologies. 

For d = 4, 5 Emparan and Reall constructed in |112j exact solutions describing a 
black hole attached to a Kaluza-Klein bubble using a generalized Weyl ansatz, describing 
axisymmetric static space-times in higher-dimensional gravitvF^I For d = 4 this was gen- 
eralized in |114j to two black holes plus one bubble or two bubbles plus one black string. 
There, it was also argued that the bubble balances the gravitational attraction between 
the two black holes, thus keeping the configuration in static equilibrium. 

In [50] these solutions were generalized to a family of exact metrics for configurations 
with p bubbles and q = p,p ± 1 black holes in D = 5,6 dimensions. These are regular 
and static solutions of the vacuum Einstein equations, describing alternating sequences of 
Kaluza-Klein bubbles and black holes, e.g. for (p, q) = (2, 3) there is a sequence of the 
form: 

black hole — bubble — black hole — bubble — black hole . 

These solutions will be called bubble-black hole sequences and we will refer to particular 
elements in this class as (p, q) solutions. This large class of solutions, which was anticipated 
in Ref. [112J, includes the (1,1), (1,2) and (2,1) solutions obtained and analyzed in [112, 
114J as special cases. 

14 See Ref. [1131 for the generalization of this class to stationary solutions. 
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We refer the reader to [50J for the explicit construction of these bubble-black hole 
sequences and a comprehensive analysis of their properties. Here we list a number of 
essential features: 

• All values of n in the range l/(d — 2) < n < d — 2 are realized. 

• The mass [i can become arbitrarily large, and for [i — > oo we have n — > l/(d — 2). 

• The solutions contain bubbles and black holes of various topologies. In the five- 
dimensional case we find black holes with spherical S* 3 and ring S 2 x S 1 topology, 
depending on whether the black hole is at the end of the sequence or not. Similarly, 
in the six-dimensional case we find black holes with ring topology S 3 x S 1 and tuboid 
topology S 2 X S X S , depending on whether the black hole is at the end of the 
sequence or not. The bubbles support the S^'s of the horizons against gravitational 
collapse. 

• The (p, q) solutions are subject to constraints enforcing regularity, but this leaves q 
independent dimensionless parameters allowing for instance the relative sizes of the 
black holes to vary. The existence of q independent parameters in each (p, q) solution 
is the reason for the large degree of non-uniqueness in the (//, n) phase diagram, when 
considering bubble-black hole sequences. 

Another interesting feature of the bubble-black hole sequences is that there exists a map 
between five- and six-dimensional solutions [50]. As a consequence there is a corresponding 
map between the physical parameters which reads 

M (M» = (5 - n&») , = 5ni5D) ~l , (3.21) 



t (6D) _ t (5D) JfiD) _ 4tt (fiD) 

where the superscripts 5D and 6D label the five- and six-dimensional quantities respec- 
tively. This form of the map assumes a certain normalization of the parameters of the 
solution, or equivalently, a choice of units, as further explained in Ref. [50] . 

3.3.3 Thermodynamics 

For static space-times with more than one black hole horizon we can associate a temper- 
ature to each black hole by analytically continuing the solution to Euclidean space and 
performing the proper identifications needed to make the Euclidean solution regular near 
the locations of the Wick rotated event horizons. The temperatures of the black holes need 
not be equal, and one can derive a generalized Smarr formula that involves the tempera- 
ture of each black hole. The Euclidean solution is regular everywhere only when all the 
temperatures are equal. It is always possible to choose the q free parameters of the (p, q) 
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solution to give a one-parameter family of regular equal temperature solutions, which we 
denote (p, g)t- 

The equal temperature (p,q)t solutions are of special interest for two reasons. First, 
the two solutions, (p, g)t and (q,p)t, are directly related by a double Wick rotation which 
effectively interchanges the time coordinate and the coordinate parameterizing the Kaluza- 
Klein circle. This provides a duality map under which bubbles and black holes are inter- 
changed, giving rise to the following explicit map between the physical quantities of the 
solutions 

/jf = nt d "V , n' = - , t' = - , s' = {d ~ 2)n ~\ d - 1 5 . (3.23) 
n t a — 2 — n 

Secondly, for a given family of (p, q) solutions, the equal temperature solution extremizes 
the entropy for fixed mass \x and fixed size of the Kaluza-Klein circle at infinity. For all 
explicit cases considered we find that the entropy is minimized for equal temperatures F^l 
Furthermore, the entropy of the (1, 1) solution is always lower than the entropy of 
the uniform black string of the same mass \x. We expect that all other bubble-black hole 
sequences (p, q) have entropy lower than the (1, 1) solution, and this is confirmed for all 
explicitly studied examples. The physical reason behind the expectation that all bubble- 
black hole sequences have lower entropy than a uniform string of the same mass, is that 
some of the mass has gone into the bubble rather than the black holes, giving a smaller 
horizon area. 



3.3.4 Phase diagrams for d = 4 and d = 5 

The general (/i, n) phase diagrams for d = 4, 5 can in principle be drawn with all possible 
values (p,q), though the explicit solution of the constraints becomes increasingly compli- 
cated for high p, q. However, the richness of the phase structure and the non-uniqueness 
in the (/i, n) phase diagram, is already illustrated by considering some particular examples 
of the general class of solutions, as was done in [50J. As an illustration, we give here the 
exact form of the curve for the (1,1) solution, corresponding to a bubble on a black hole, 



d = 4: n (1)1) (/i) 



I -'- ; m/1 + ^2 



d = 5: n(lil) ( M ) = I + _. (3.24) 



in five and six dimensions respectively. These two solutions are related by the map in 
(|3.2ip and one may also check that these curves are correctly self-dual under the duality 
map dS2B|> . 

In Figured] we have drawn for d = 4 and d = 5 the curves in the (//, n) phase diagram 



lo This is a feature that is particular to black holes, independently of the presence of bubbles. As an 
analog, consider two Schwarzschild black holes very far apart. It is straightforward to see that for fixed 
total mass, the entropy of such a configuration is minimized when the black holes have the same radius 
(hence same temperature), while the maximal entropy configuration is the one where all the mass is located 
in a single black hole. 
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for the (p, q) = (1, 1), (l,2) t and (2, 1) solutions. These correspond to the configurations 



(1,1) : black hole — bubble 

D = 5 S 3 D (3.25) 

D = 6 S 3 x S 1 Dx S 1 

(1,2) : black hole — bubble — black hole 

D = 5 S 3 S 1 x I S 3 (3.26) 

D = 6 S 3 x^ T 2 x I S 3 x S 1 



(2, 1) : bubble — black ring — bubble 

D = 5 D S 2 x S 1 D (3.27) 

D = 6 Dx S 1 S 2 xT 2 Dx S 1 

where the first/second line in each configuration corresponds to the topology in five/six 
dimensions. Here D denotes the disc and I the line interval. 

n n 




Figure 5: (fi, n) phase diagrams for five (left figure) and six (right figure) dimensions. 
We have drawn the (p, q) = (1,1), (l,2)t and (2,1) solutions. These curves lie in the 
region 1/2 < n < 2 for the five-dimensional case and 1/3 < n < 3 for the six-dimensional 
case. The lowest (red) curve corresponds to the (1, 1) solution. The (blue) curve that 
has highest n for high values of [i is the equal temperature (1, 2){ solution. The (green) 
curve that has highest n for small values of fi is the (2, 1) solution. The entire phase 
space of the (1,2) configuration is the wedge bounded by the equal temperature (l,2)t 
curve and the (1, 1) curve. For completeness we have also included the uniform (orange) 
and non-uniform (cyan) black string branch, and the small black hole branch (magenta) 
displayed in Figure [3l 
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3.3.5 Non-uniqueness in the phase diagram 

Finally we remark on the non-uniqueness of Kaluza-Klein black holes in the (//, n) phase 
diagram^ Clearly for a given mass there is a high degree of non-uniqueness. The non- 
uniqueness is not lifted by taking into account the relative tension n, as there are explicitly 
known cases of physically distinct solutions with the same mass and tension. For example 
the (1, 2)t solution and the (2, 1) solution intersect each other in the phase diagram. This 
means that we have two physically different solutions in the same point of the f/x, n) phase 
diagram. Moreover, there is in fact a continuously infinite non-uniqueness^ for certain 
points in the (/x, n) phase diagram. This is due to the fact that the (p, q) solution has q 
free parameters [50J. Hence, for given p > 2 and q > 3, we have q — 2 free continuous 
parameters labeling physically different (p, q) solutions, for certain points in the (/x, n) 
phase diagram. 

4 More on the Gregory-Laflamme instability 

In this section we present further results on the Gregory-Laflamme instability. First we 
discuss the large dimension limit of the GL critical mass and the phase diagram that 
emerges. The case of the GL instability for boosted neutral black strings, which has 
various interesting applications, will be discussed next. We also summarize some of the 
main results connecting the GL analysis to an LG (Landau-Ginzburg) analysis of the 
thermodynamics and what is known about the GL instability for neutral black branes with 
higher-dimensional compact spaces. Finally, we comment on an interesting connection 
between the GL instability and a classical membrane instability, known as the Rayleigh- 
Plateau instability. 

4.1 Large d limit 

In Section [2~T2l we reviewed the GL instability of the neutral black string in d+1 dimensions, 
and in particular presented the differential equations satisfied by the threshold mode. 
Unfortunately, no analytic solution to these equations is known and one has to solve them 
numerically for a given value of d on a case by case basis. This was done for d up to 
50 in Ref. [31]. It is interesting to regard d as a parameter in this system and study the 
behavior of the GL mode as d becomes very large. 

In Ref. [31] it was observed from the data that the critical mass exhibits the exponential 
behavior 



16 Non-uniqueness in higher dimensional pure gravity has also been found for stationary black hole 
solutions in asymptotically flat space-time. In five dimensions for a certain range of parameters there 
exists both a rotating black hole with S 3 horizon [14] and rotating black rings with S 2 x S 1 horizons |15| . 

17 Infinite non- uniqueness has also been found in [115] for black rings with dipole charges in asymptotically 
flat space. 



/x GL ~ 16.21 • (0.686) d 
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at large d. One can verify this result by solving analytically for the threshold mode in 
the large d limit, as was done in Ref. [35]. Combined with the behavior of small localized 
black holes for large dimensions, this result enables us to conjecture the way the localized 
and non-uniform phases appear in the (//, n) phase diagram when the dimension is large. 

To derive the large d behavior of the GL mode, it is convenient to start with the 
linear second order ODE that appears when considering the negative modes of a static 
and spherically symmetric perturbation of the Schwarzschild black hole in d dimensions 
(see (|2.29p ). By taking the large d limit, this ODE simplifies considerably and by dropping 
the exponentially growing solution, one findi^l 

X(r) ~X:d =1 (r)r" i t i (4.2) 

2 

where KL S is the modified Bessel function of the second kind (we remind the reader that 
the function x( r ) 1S related to the h rr component of the perturbation, c.f. eqs. (12. lip . 
(I2.12j) ), The GL wave number A;gl — Vd can be determined by imposing the appropriate 
boundary conditions and the resulting GL mass takes the form [35j 




log/iQL ^ dlogW— , d>l (4.3) 



This analytical result meshes very nicely with the earlier numerically obtained result (|4.ip . 

Moreover, it is interesting to consider at large d the slope of the localized black hole 
branch (|3.17p . measured in units normalized with respect to the Gregory-Laflamme point 
(/•'GL) n GL = l/(d — 2)), with //gl given by the large d formula (|4.3|) . One finds 

n . . jd/2 V 



d d/2 — , d» 1 . (4.4) 

n>GL MGL 

From this expression we see that the slope becomes infinitely steep as d —* oo. This 
suggests that the curve describing the localized black hole and non-uniform black string 
branches should behave as sketched in Figure El In particular, the large d behavior of the 
phase structure is expected to be such that the non-uniform string gets closer and closer 
to having n = ncL = — 2). 

Further evidence for the above picture stems from the realization that the behavior 
(|4.3j) is equivalent to the critical Schwarzschild radius being equal to ro ~ L/(2ir)Vd. This 
shows that the large d critical string is "fat" , and that most probably it cannot decay into 
a black hole since the black hole cannot fit into the compact dimension at the GL point 
|35| . Hence, the uniform black string should decay towards a different end-state, which 
presumably is the non-uniform black string. 

Finally, from the comparison of the GL instability with the Jeans instability in Section 
2.21 we recall that the Jeans instability predicts that the large d behavior of the critical 
wave number is kj ~ d 3 / 4 . Now we see that this is different from the result k c ~ \[d 
above. The classical Rayleigh-Plateau instability, on the other hand, does predict the 
correct large d behavior as we will review in Section 14.41 



18 Equivalently, one can take the large d limit of the second order differential equation (IA.9|l 
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Figure 6: Sketch of the (/i, n) phase diagram for large d. 



4.2 Boosted black strings 

Another interesting direction that has been pursued in [116] is the study of the GL insta- 
bility for boosted black strings. As we briefly review below it was found in |116j that the 
critical dimension (see Section 13. 2p depends on the internal velocity and disappears for 
sufficiently large boosts. Moreover, the five-dimensional black ring [15] reduces in the large 
radius limit to a boosted neutral black string (in five dimensions), so these instabilities are 
also relevant to the stability analysis of black rings (we will briefly return to those issues 
in Section [87L]) . 

Boosted black strings can be obtained by boosting the static black string (|2.9p along 
the z-direction. The resulting solution is 

/ r d ~ 3 \ r d ~ 3 ( r d ^\ 

ds 2 = — I 1 — cosh 2 a °_ 3 j dt 2 + 2 ° d _ 3 cosh a sinh adtdz + I 1 + sinh 2 a ° d _ 3 j dz 2 



„d-3\ 1 



yd— 3 



dr 2 + r 2 (d6 2 + sin 2 6dn 2 d _ 3 ) . (4.5) 



The mass and tension of the boosted black string can be computed from the general 
expressions (|3.2p . (|3.3|) . (13.51) . giving 



„ /ro\' 1 " 3 , „. ,9 „, (d — 3) sinh 2 a — 1 . 

H = Sl±-2[-f) P -3) cosh 2 a + 1] , n=)— — { (4.6) 

V L J (a — 3) cosh a + 1 

In this case the string also carries the dimensionless moment urrF^I 

p = r^_ 2 (<i — 3) cosh a sinh a . (4-7) 

One can now extend the original global thermodynamic argument of Gregory and 
Laflamme to this case, by comparing the entropy of the boosted black string (with horizon 
ro and boost a) to the entropy of a boosted spherical black hole (with horizon r' Q and boost 



19 The momentum is given by P — T tz — j§^§^ c tz with c tz the first correction around flat space in 
g tz — 1 + c tz /r d ~ 3 . We define, in analogy with (|3.5[1 . the dimensionless momentum as p = 1 ^ G " P. 
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a') with the same energy and momentum. Because we compare the solutions at the same 
energy and momentum, one can express (r' ,a') in terms of (ro,a) of the boosted black 
string. Then, equating the two entropies determines a minimal circumference |116j 



mm cosha[( ( i-3)(d-l)+cosh- 2 a](' i - 1 )/ 2 O d _ 2 { ' } 

and hence from (|4.6p one can define a maximum mass /U max = ^L=L min - This mass can be 
seen as a rough estimate of the point where the boosted black string becomes unstable. 
The instability appears for [i < ^ max (or L > L mm ). Note that L m \ n scales like 1/ cosh a, 
so this naive analysis suggests that the instability will persist for large boosts. 

One can further analyze the instability of the boosted black string by appropriately 
boosting the time-dependent GL mode [116J. In this way, one finds that the instabilities 
of the boosted string can be related to instabilities of the static black string with complex 
frequencies. The relation is even more direct for the static threshold mode (j2.15|) . on 
which one can simply apply the same boost that gave (|4.5|) . As a result one finds the 
mode exp(iu;t + ik c ,\>oostZ) with 

fcc,boost = cosh ak c , u) = smhak c (4-9) 

where k c is the threshold wave number for the static black string (see Table [1]). So the 
threshold modes are waves that travel in the z-direction with the same speed as the 
boosted string. Note that the boost of the GL mode in a transverse direction, including 
the time-dependent one, was considered in Refs. [56l [66] and will be reviewed in Section 

Given the above results one can show that the critical dimension (see Section 13.2ft 
depends on the boost parameter, and in fact disappears (for any given d) for large enough 
boost [116] . To see this observe from (|4.9|) that 



-^GL.boost 



cosh a 



fJ-GL 



and hence from (14.6 



cosh d 3 a .. , „. . 2 in / i -i -i \ 

/^GL.boost = d _ 2 [{a - 3) cosh a + 1\hgl (4.11) 

where is the neutral GL mass. Comparing this to the //max 

obtained above (see 

below (|4.8|) ) one finds that for small boosts /i max > AiGL.boost so that there is a regime 
where the boosted black string is stable, but the localized black hole has higher entropy, 
suggesting the existence of an unstable non-uniform black string in between. For zero 
boost this occurs for dimensions D < 13. On the other hand, for sufficiently large boosts 
one has /i max < jUGL,boost- So in this case, there is a regime, ^ max < \i < /JGL,boost, where 
the boosted black string is unstable but the localized black hole does not have a higher 
entropy, suggesting that the end-state is a stable non-uniform (boosted) black string in 
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this range. In other words, for sufficiently large boosts the critical dimension disappears 
and we have a phase diagram that is analogous to that found for large d in the neutral 
case in Fig. [6] 

We will return to some other aspects of boosted black strings in Section 18.11 
4.3 Higher-dimensional compact spaces 

Most of the work on the stability and phases of black objects in spaces with compact 
transverse directions has been performed so far for the simplest case of one compact 
direction. Already for this case one observes a very rich phase structure. This phase 
structure is expected to become even richer as we go to higher dimensional compact 
spaces, such as T n , S n , K3 and CY% manifolds. As the number of moduli increases, the 
structure of the moduli space and hence the phase structure is clearly going to become 
more involved. 

4.3.1 Instability of neutral black p-branes 

The one higher-dimensional case that has been studied |35|. [37] in more detail is that of 
the torus T p , which is the simplest generalization of the circle S . We summarize some of 
the main results here. 

Consider neutral solutions of Einstein gravity with an event horizon, that asymptote 
to M*- 1 ' 1 x JP. Here we have the uniform black p-brane phase which is the obvious 
generalization of the black string, and is just a Schwarzschild-Tangherlini black hole in d 
dimensions, with p flat directions added ds 2 = X^f=i dz l dz l for the p-dimensional torus. 
The torus is defined by p period vectors ei such that any function Y(z + e) = Y(z) with 
e belonging to the integer lattice spanned by the period vectors. Performing a similar 
stability analysis as in Section \2. 2\ one finds that the threshold mode P CtfJiU (r)Y \z) satisfies 
[35] 

A S ch + A T P=0 (4.12) 

where As c h and Atp are the eigenvalues in 

A L PcM r ) = Asch^(r) , -AY(z) = \ TP Y(z) (4.13) 

The only way to satisfy (|4.12p is to take the GPY mode for the Schwarzschild black hole, 
and then it suffices to consider modes which are tensor on the Schwarzschild geometry 
and scalar on the torus. Consequently, we find that the instability of the black brane sets 
in when 

rl = -M- (4.14) 

where k c is the eigenvalue in (|2.21j) and A m j n is the mininimal (non-zero) eigenvalue of 
the Laplacian on the torus T p . In fact, this argument is general and applies to any stable 
compact manifold. For the case of T p this eigenvalue is given explicitly by 

Amin.TP = |femin| 2 (4-15) 
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where k m \ n is the shortest vector in the reciprocal lattice of TP (i.e. the lattice spanned by 
the vectors k with k ■ e = 2ttti, n £ Z). 

It is also physically clear that for a highly asymmetric torus, with e.g. one large di- 
rection and the other ones small, the instability will set in for the mode along the large 
direction. In this sense the stability properties of black branes on such asymmetric tori 
are closely related to those of the black string. Hence, it seems more plausible that new 
effects will show up in the study of the other extreme, namely a square torus, with all 
directions of equal length. Moreover, understanding these two extremes may lead to a 
qualitative understanding of the intermediate region. 

Recently, the Landau-Ginzburg method (see Section 14. 4p . was applied to the case 
of black p-branes on a square torus T p . In this case it is easy to see that p modes 
(corresponding to each of the directions) become marginally tachyonic at the same GL 
point. Ref. [37] then studied the order of the phase transition, i.e. the critical dimension, 
for such tori. It was shown that for the phase transition to be first order, it is enough 
to find a single direction in tachyon space for which this is the case. It was furthermore 
shown that in order to study the T p case one needs to combine the known results for the 
S 1 case together with results for the T 2 . As a result it was concluded that the transition 
order for square tori depends only on the number of extended dimensions d, and not on 
p. In other words, the critical dimension discussed in Section 13.21 for p = 1 is unchanged 
when considering higher-dimensional square tori. Note that in an earlier work [35] (see 
also [117] ) a more naive analysis based on the comparison of the entropy of black holes 
localized in the compact directions and the entropy of black branes seemed to indicate that 
the critical dimension could become lower as the number of compact directions increases. 
As explained in [37J, this conclusion is misleading because the corrections to the entropy 
of localized black holes can be shown to have an important effect in this case, and these 
were ignored in the analysis of [35] . More generally, this shows that one should be cautious 
when using the Schwarzschild black hole entropy as an approximation to make qualitative 
statements, i.e. one should carefully study the effects of corrections due to the presence of 
the compact manifold. 

4.3.2 Phase structure for higher tori 

So far we have just considered the aspect of instability of uniform black branes with 
compact transverse directions, but if we think about possible phases it is not difficult to 
see that this will be very rich. Consider, for example, the case of the two-torus. First of 
all, for any phase of Kaluza-Klein black holes on M} ,d ~ 1 x S 1 we trivially get a phase for 
the case of .M 1,d— 1 x T 2 by adding a flat compact direction. This includes the uniform 
black two-brane phase discussed above, but also e.g. semi-localized black holes of horizon 
topology S d ~ l x S 1 (these are of black string type) as well as bubble-black hole sequences 
with an extra circle direction. There are also going to be entirely new phases each time we 
add a compact direction. For example, in the case discussed above, there is also a phase of 
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black holes that are localized in two of the T 2 -directions, i.e. a phase with horizon topology 
S d . These solutions reduce to the (d+2)-dimensional Schwarzschild black hole in the limit 
of zero mass, and it would be interesting to obtain the first order correction away from 
the zero mass limit using the methods of [38j [39] or [3D] . There are also new bubble-black 
hole sequences for the case of d = 5, 6 (see for example the Conclusion of |50j) that are in 
principle not difficult to construct explicitly. All these will be part of some complicated 
higher dimensional phase structure with a complicated pattern of (in)stabilities. 

4.3.3 Comments on general compact Ricci fiat manifolds 

More generally we may consider KK black hole solutions in (d+p)-dimensional pure gravity 
asymptoting to x Af p , with Af p being a generic p-dimensional compact Ricci flat 

manifold, for example a compact Calabi-Yau manifold. In this case, the phase diagram of 
neutral, static black hole solutions will be parameterized by the mass A4, the components 
of a generalized tension tensor T and the moduli that characterize the manifold M p at 
infinity (minus one for the overall volume of M p that can be scaled away). The components 
of the tension tensor T are straightforward generalizations of the single tension n appearing 
above and can be thought of as the chemical potentials corresponding to the moduli of 
M p . See Ref. [104J for a treatment of the two-torus case, including a derivation of the first 
law of thermodynamics. 

To obtain the uniform black membrane solution we simply tensor the d-dimensional 
Schwarzschild solution with the Ricci flat metric of AC, 



An obvious generalization of the argument appearing around eq. ([4.12ft suggests again the 
presence of a Gregory-Laflamme instability. The discussion of the higher tori case above 
indicates the presence of many more stable or unstable phases in the phase diagram with 
a complicated pattern of (in)stabilities. Beyond the phases of localized black holes, non- 
uniform black membranes and bubble-black hole sequences one can imagine new "uniform" 
phases, where the Af p moduli are radially varying, the tension tensor T ab has values other 
than those of the uniform black brane solution (|4.16p and the metric components have no 
dependence on the coordinates of N p . Black hole solutions of a similar spirit are abundant 
in extended supergravity, for example in d = 4, J\f = 2 supergravity |118[ 1119] (see [120] 
for a recent review). 

In general, many features of the phase diagram are expected to be strongly correlated 
with the topology of the compact manifold M p . It would be interesting to explore this 
connection further. 




(4-16) 
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4.4 Other developments 

We conclude this section by briefly discussing several other developments related to the GL 
instability and the Kaluza-Klein black hole phase diagram. These include an application 
of the Landau-Ginzburg theory of phase transitions, the use of Morse theory and the study 
of the merger point, where the black string and black hole branches meet. While these 
developments are important and very interesting, we choose to be brief here since many 
of these are very nicely explained in the review [27] by Kol. We also comment on a recent 
connection observed in Ref. |36j between the (classical) Rayleigh-Plateau instability of 
membranes and the GL instability. 

4.4.1 LG in GL 

Some of the features of the black hole-black string transition, in particular the order of 
the phase transition, can be nicely studied using the Landau-Ginzburg (LG) theory for 
the thermodynamics, as was done in Refs. [Ml EH EZ] (see also [291 E01 EI]). We briefly 
summarize some of the main results here, referring the reader to the review [27] for a useful 
expose of the general idea and further details, as well as the recent article [37] where a 
new method for determining the order of the phase transition was developed based on the 
LG perspective. 

In statistical physics the Landau-Ginzburg theory [121] determines the order of a phase 
transition from the local behavior of the free energy F near a critical point. For a simple 
system parameterized by a single variable A and a control parameter [i the free energy is 
a function F = F(X;fi). We are interested in systems where a certain parity symmetry 
forces the free energy to be an even function of A, i.e. F = F(X 2 ; fi). For such systems the 
existence of a marginally tachyonic mode at some critical value [i c of the control parameter 
implies that we can expand F near A = 0, jj, = /u c in a power series in A 

F = F 1 {n- fi c )X 2 + F 2 A 4 + 0(A 6 ) , (4.17) 

where without loss of generality we can take the constant F\ to be positive. The two 
possible behaviors of F are summarized in Figures [5] and are distinguished by the sign 
of the quartic coefficient F2 . 

For negative F2 (see Fig. [7J the system undergoes a first order transition. As [i is 
varied the system passes from a symmetric stable phase (A = and fi > /i c ) to an unstable 
symmetric phase. Since F is bounded from below there is a new stable asymmetric vacuum 
at finite A and the system discontinuously jumps to this new stable vacuum when fi < fj> c . 

For positive F2 (see Fig. [8]) the system undergoes a second (or higher) order phase 
transition. As we vary \x from fi > fi c to [i < fi c the system passes continuously from the 
stable symmetric vacuum at A = to a new asymmetric vacuum at A 7^ 0. 

The non-uniform black string exhibits similar behavior with details that depend cru- 
cially on the space-time dimension. In the language of the above one-dimensional example, 
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Figure 7: The typical behavior of the free energy during a first order phase transition. 
The red lines are at constant values of n > \i c and the green lines at constant [i < \x c . 
The thick black line represents a set of unstable asymmetric (A / 0) vacua at \i > jx c . 
The dashed horizontal line represents the symmetric (A = 0) phase of the system, which 
is stable at [i > [i c and unstable at fj, < [i c . 

we can think of A as a parameter that measures how far the black string is from the uniform 
solution. This role can be played by the amplitude of the marginally tachyonic GL mode 
(I2.2ip . i.e. A can be taken as the coefficient of the threshold mode metric perturbation 



The control parameter \i can be taken as the mass of the black brane solution, or better 
yet in the canonical ensemble as the dimensionless ratio (3/L of the inverse temperature 
over the z radius. 



When dealing with gravity one has to consider the free energ; 




as a function on the infinite dimensional configuration space of all metrics with given 
boundary conditions. For black branes smeared on a transverse T p space the correspond- 
ing Landau-Ginzburg analysis around the uniform solution has been performed in [37J. 
In this system one finds a generalized multi-dimensional version of (|4.17p . which arises 
by analyzing the geometry of the non-uniform black brane up to second order in pertur- 
bation theory (for more details we refer the reader to [37]). The order of the transition 
is determined by the sign of an appropriate set of quartic coefficients in the free energy 
expansion. One finds that the result depends crucially on the dimension d = D — p of the 
non-compact part of the space time (see also comments in Section [4.3p . For D < 12 the 

20 The r.h.s. of this equation is the standard Gibbons-Hawking expression for the gravity action. R is the 
Ricci scalar and K the extrinsic curvature on the boundary dM. K° is the same quantity for a reference 



h ~ A e 



,ik c z 



Pc ■ 



(4.18) 




(4.19) 



geometry. 



41 




A 



Figure 8: The typical behavior of the free energy during a second (or higher) order phase 
transition. Again, the red lines are at constant values of \i > [i c and the green lines at 
constant \i < ji c . The thick black line represents a set of stable asymmetric (A 7^ 0) vacua 
at \i < \i c . The dashed horizontal line represents the symmetric (A = 0) phase of the 
system, which is stable at \i > /j, c and unstable at /x < /x c . 

transition is first order and the phase diagram looks like Fig. [3l whereas for D > 13 the 
transition is second order and the phase diagram becomes as in Fig. [6l Hence, there is 
an interesting critical dimension U D* = 12.5" where the behavior of the system changes 
dramatically. 

The order of the phase transition for a black string on M 3 ' 1 x S was first determined 
with a different method by Gubser in [29] and was improved upon by Wiseman in |30j. 
This computation was generalized to a more general background Ai d x S 1 in [31] where the 
critical dimension D* was first observed. In addition, |122j pointed out that the critical 
dimension depends also on the ensemble; in the micro-canonical ensemble the critical 
dimension is one unit higher than that in the canonical ensemble. 

4.4.2 Morse theory 

Another useful method that has been applied to address the issue of the endpoint of the 
decay of the unstable black string is Morse theory [123j . as originally proposed in |34j 
and reviewed in [27]. While Landau-Ginzburg theory can be employed to study the local 
properties of the phase diagram, Morse theory can be used to obtain information about 
the global (topological) structure of the phase diagram. 

In general, when we try to determine the pattern of static solutions in the phase 
diagram of Kaluza-Klein black holes on A4 d x S 1 , we are in essence studying the extrema 
of the gravitational action in the space of metrics with given boundary conditions. Morse 
theory can be useful in this study, since it is, by definition, the topological theory of the 
extrema of functions on a given manifold. For instance, it can be used to put constraints on 
how the extrema of a function behave under continuous deformations. The latter aspect of 
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Morse theory was employed in [M] to argue that a natural completion of the black string 
phase diagram involves a point, called the "merger point", where the non-uniform black 
string branch meets the localized black hole branch. 

The intuitive argument behind this analysis is based on the fact that we know that 
there is a localized black hole phase as well as a non-uniform black string phase, but a 
priori not where these phases end in the phase diagram. Since it is unphysical for a phase 
to end in nothing, it follows by "phase conservation" that the simplest possibility is that 
these two phases meet in a topology changing phase transition, the merger point. 

4.4.3 Merger point 

The prediction of a merger point in the phase diagram of black strings implies an intrigu- 
ing topology changing transition where the black string horizon with topology S d ~ 2 x S 1 
pinches off and becomes a black hole horizon with topology S d ~ 1 . This transition reminds 
of the Gross- Witten transition [124j in two dimensional large N gauge theories. Also cer- 
tain features of the topology change near the merger point resemble the conifold transition 
in the moduli space of Calabi-Yau threefolds [125J. 

To this effect, it was argued in Ref [34] that not only does the horizon topology change, 
but also the local topology in the Euclidean section. This topology change can be modelled 
in analogy to the conifold transition as a pyramid, with the singular topology being the 
cone over S d ~ 2 x S 2 , 

n 2 

ds 2 = dp 2 + -^-j [dn 2 2 + (d- 3)dnj_ 2 } (4.20) 

which was termed the "double-cone". Here, p measures the distance from the tip of the 
cone and the constant prefactors ensure Ricci- flatness. Refs. |52[ 133] show that the local 
geometry around the "waist" of the most nonuniform solutions is indeed cone-like to a 
good approximation. 

Evidence for the merger point transition has been found with a numerical study in 
|52[l33j. Interesting parallels between the merger point transition and the Choptuik scaling 
in black hole formation [72] have been observed in [68], [691 EQl [TT] . More details about the 
merger point can be found in the original work [33] and the review |27j . 

4.4.4 Rayleigh-Plateau instability 

The first law of thermodynamics 

dM = TdA , (4.21) 

A being the horizon area, suggests a natural parallel between the mechanics of black holes 
and hydrodynamics. Indeed, by viewing (|4.2ip as a law for fluids where T is interpreted 
as an effective surface tension one is led to a simple picture of black holes, where the 
event horizon is treated as a kind of membrane with well-defined mechanical, electrical 
and magnetic properties [126] . It has been noticed that the membrane approach works 
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surprisingly well in reproducing various aspects of black hole physics in different contexts 
(see for instance [127} 1128] for an application in the context of AdS/CFT). 

The membrane analogy was applied to black branes in |36[ 1129] , where it was observed 
that there is an impressive correspondence between the Rayleigh-Plateau instability of 
long fluid cylinders and the Gregory-Laflamme instability of black strings |36j. Both 
instabilities are long-wavelength instabilities involving only the axisymmetric s-modes (for 
the black string this statement has been verified in [130] ). The critical wave numbers 
kftp and k c agree exactly at large dimension d - in particular they both scale as yd 
for d > 1 - and there is qualitative agreement between the behavior of the instability 
time scales. Hence, this classical system appears to provide a better description of the 
Gregory-Laflamme instability as compared to the Jeans instability (see Section 12. 2\\ . 

5 Instabilities in supergravity 

Our discussion of the GL instability has so far been entirely within the context of Einstein 
gravity. A natural question is whether GL-like instabilities can also arise in the presence 
of matter fields, such as scalar fields and gauge potentials. This question is especially 
relevant in supergravity theories which appear in the low-energy limit of string theory. In 
this context one encounters a variety of black brane solutions whose stability properties 
are not only important in their own right, but can also have non-trivial consequences for 
their dual non-gravitational theories via the gauge/gravity correspondence (see Section [7|). 

Gregory and Laflamme already considered instabilities of charged branes in j3[ 1131] . 
Further work in this direction appeared in [J32J ES [631 El E3 ESI EH1 CESS UM [1371 

The investigation into the classical stability of charged branes began with the study of 
non-extremal singly-charged black p-branes [3j I131[ |64"1 1135[ I136[ I138[ 1141] . In this case, 
there is a GL instability when the charge is sufficiently small compared to the total mass 
of the brane, whereas the brane is stable for large charge, i.e. close to extremality. This 
has been shown using numerical methods. In Section [6] we comment on how this is related 
to the thermodynamics of these branes. In this review we consider instead the case of 
smeared non-extremal singly-charged black p-branes. The reason is that in this case we 
can be more precise, since we can explicitly find the GL instability by mapping it from 
the GL instability of the neutral black string [56, 66j. 

Before turning to the GL instability of non- and near-extremal smeared branes, we 
review the boost/U-duality map between phases of neutral Kaluza-Klein black holes and 
this class of branes [57J (see also [SSJ [2JJ 1122] ) . This map suggests the presence of a GL 
instability in the case of non- and near-extremal smeared branes, which we then proceed 
to verify more explicitly. We also review recent work on GL instabilities of D-brane bound 
states that appeared in Refs. [59 [ [60l [6Tj . 
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5.1 New phases of non- and near-extremal branes 

The starting point of our discussion is the set of 1/2 BPS branes in type IIA/B string 
theory and M-theory with a circle in the transverse space. We would like to consider 
the class of solutions that arise as thermal excitations of these branes. We will refer to 
the branes in this class as non-extremal branes on a circle. The near-extremal limit of 
these solutions will be defined in a way that keeps the non-trivial physics related to the 
presence of the circle. We will refer to this particular class of near-extremal branes as near- 
extremal branes on a circle. The dual non-gravitational theories living on these branes 
will be reviewed briefly in Section [7J 

In analogy with the neutral Kaluza-Klein black holes of Section El one can define phase 
diagrams for the non- and near-extremal branes on a circle. See Ref. [57] for details and 
conventions. For non-extremal branes on a circle the dimensionless asymptotic quantities 
are the mass p,, the relative tension in the circle direction n and the charge q. For near- 
extremal branes on a circle there are two quantities, the energy above extremality e, 
rescaled to be dimensionless, and the relative tension r which is the ratio between the 
tension of the brane along the transverse circle and the energy. The tension is measured 
using the general tension formula found in [97] that holds for non-asymptotically flat 
spaces. Hence, the resulting phases can be plotted in a two-dimensional (e, r) phase 
diagram. 

In this section we review the following: 

• One can transform any static and neutral Kaluza-Klein black hole to a non-extremal 
brane on a circle with a combined boost and U-dualitJ^I transformation. For a given 
charge q, this gives a map from points in the (//, n) phase diagram to points in the 
non-extremal (/2, n) phase diagram. 

• By taking a particular near-extremal limit of the map to non-extremal branes, one 
can transform any static and neutral Kaluza-Klein black hole to a near-extremal 
brane on a circle. In particular, this provides a map relating points in the (/x, n) 
phase diagram to points in the (e, r) phase diagram. 

• The above maps can be applied to the known phases of neutral Kaluza-Klein black 
holes, from which we generate (new) phases of non- and near-extremal branes on a 
circle. 

5.1.1 Extremal p-branes 

We consider here a class of branes that are thermal excitations of the 1/2 BPS branes of 
type IIA and IIB string theory and M-theory. The singly-charged dilatonic p-branes in a 

21 See e.g. Ref. |142] for a review on U-duality. 
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D-dimensional space-time are solutions to the equations of motion of the action 

'° - isk I d ° x ^ ( R ~ ~ ra e "* (F " +2 » )2 ) • <5 ' 1) 

where cfi is the dilaton field and F(p+2) is a (p+2)-form field strength with -F( p +2) = ^(p+i) 
and -A( p +i) the corresponding (p + l)-form gauge field. For D = 10 the action (|5.ip is the 
bosonic part of the low energy action of type IIA and IIB string theory (in the Einstein- 
frame) when only one of the gauge- fields is present. For D = 11 it is the low energy action 
of M-theory for suitable choices of a and p. 

The extremal 1/2 BPS p-brane solutions in string and M-theory take the well-known 
form 

d-2 



ds 2 = H~ — 



-dt 2 + Yji^ 1 ) 2 + Hds d 



(5.2) 



e 2<t> = H a 5 = ( H -l _ y dt A du l A . . . A du p ; ^ 5 3 ) 

with V 2 H = (away from the sources), ds^ denotes the metric and V 2 the Laplacian 
of the (f-dimensional transverse space. The solutions (|5.2p . (|5.3p represent the 1/2 BPS 
extremal p-branes of String/M-theory when D = 10, 11 and 

a . = «_ 2 C+l>C-»>, ( 5 .4) 
D-2 v ' 

With (|5.4p obeyed we get for D = 10 the D-branes (a = (3 — p)/2), NS5-brane (a = 1) and 
the F-string (a = —1) of type IIA and IIB string theory, and for D = 11 the M2-brane 
and M5-brane of M-theory (a = 0). Note that the string theory solutions are written here 
in the Einstein frame. 

In the following we place the extremal solutions (|5.2p - (|5.3p on a transverse space of 
the form x S 1 . The thermal excitations of these solutions corresponds to the non- 

extremal solutions that we find below. 



5.1.2 Charging up solutions via U-duality 

Consider a static and neutral Kaluza-Klein black hole solution (see Section [3]), i.e. a static 
vacuum solution of (d+ l)-dimensional General Relativity with at least one event horizon 
that asymptotes to M d x S 1 . We can always write the metric of such a solution in the 
form 

L 2 

ds 2 d+1 = -Udt 2 + ——V ab dx a dx b , (5.5) 

where a, b = 1, d and V a b is a symmetric tensor. U and V a b are functions of x 1 , .., x d . In 
the asymptotic region U -> 1 and L 2 (2Ti)~ 2 V ab dx a dx b describes the cylinder R^ 1 x S 1 . 
The factor L 2 /(2tt) 2 , where L is the circumference of the circle in J\A d x S 1 , has been 
included in (|5.5p for later convenience. 
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By adding p + 1 (where p = 9 — d) flat directions, this solution is trivially uplifted to a 
solution of eleven-dimensional (super) gravity. Denoting these directions by z% (i = 1 . . .p) 
and y, we perform a Lorentz boost in the y direction 



^new \ I cosh a sinh a 
Vnpw / \ sinh a cosh a 




(5.6) 



The resulting boosted metric has an isometry in the y-direction so we can make an S- 
duality to obtain a non-extremal DO-brane solution of type IIA string theory. For d < 9 
the DO-brane is uniformly smeared along an M 9_d space and has transverse space x S 1 . 
One can use further U-dualities to transform the solution into a Dp-brane solution, an F- 
string solution or an NS5-brane solution of type IIA/B String theory, or to an M2-brane or 
M5-brane solution of M-theory, depending on what dimension d we start with. In this way 
we obtain non-extremal singly-chargec@ p-branes on a transverse circle in D dimensions 
with D = d + p + 1 (for String/M-theory we have D = 10/-D = 11)@ The resulting 
non-extremal solution takes the form 

ds 2 = (-Udt 2 + V (du 1 ) 2 + t-^— # V ab dx a dx b ) , (5.7) 

V ti ^ J 

H = 1 + sinh 2 a (1 - U) , (5.8) 



e 24> _ R a ^ A^i) = coth a (H' 1 - l)dt A du 1 A • • • A du p . (5.9) 

Notice that U — > near a event horizon in the metric (|5.5p for the neutral solution. 
This translates after the map into an event horizon in the non-extremal solution (|5.7p - 
(|5.9p . Since the harmonic function (|5.8p stays finite and non-zero for U — > 0, we deduce 
that the source of the electric potential ^4( p +i) in (|5.9p is hidden behind the event horizon. 
Note also that if there is a number of event horizons defined by U = 0, the chemical 
potential v = —A tu i... uP \u = Q = tanha is the same for each disconnected component of the 
event horizon. In the following we assume that we have at least one event horizon present. 

Using the map from (|5.5|) to (|5.7p - (|5.9|) a Kaluza-Klein black hole solution in d + 1 
dimensions can be mapped to a corresponding solution of supergravity describing non- or 
near-extremal p-branes on a circle. The dimensions d and p are related by the equation 
D = d + p + 1 with D = 10 for string theory and D = 11 for M-theory. 

The map from (|5.5|) to ([5.7j) - (|5.9|) is a development of an earlier result in [3H]. In [55] 
it was pointed out that one can map any neutral solution into a corresponding non- and 
near-extremal solution in supergravity for the class of Kaluza-Klein black holes that fall 

22 By applying further boosts it is also possible to add two or three charges. For the case of d — 4 
this is considered in Ref. [77] . which considers the corresponding phases of two- and three-charge brane 
configurations on a circle. 

23 One could also use the U-duality to get branes that are smeared uniformly in some directions, but we 
choose not to consider this option here. 
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into the SO(d — l)-symmetric ansatz (|3.11|) of Ref. |38| . While this observation was made 
at the level of the equations of motion, the combined boost and U-duality transformation 
reviewed above reveals the physical reason behind the existence of this map and moreover 
shows that it works for any static and neutral Kaluza-Klein black hole. 

Here we present a summary of the map between physical quantities of the non-extremal 
branes and physical quantities of the neutral Kaluza-Klein black holes (we refer the reader 
to [57] for more details on the boost/U-duality map). In terms of the definitions of Ref. [57J 

b 2 q , d — 2 — n fi 

~Tj~q 



/i 



(d + n) 



n 



n 



(5.10) 



2(d-i)" 1 i + VTT¥ ' " '"' "~ 2(d 

where (n,n) (see (13. 5p . (13.31) ) are the dimensionless mass and relative tension of the neu- 
tral Kaluza-Klein black hole and (jl, n, q) the dimensionless mass, relative tension and 
dimensionless charge of the non-extremal brane. For the temperature and entropy one 
finds 



t 



cosh a 



, s = cosh as, , cosh a(fi,n,q) 



1 + b + Vl + b 2 



2Jb{l + VTTb^) 



(5.11) 



where (t, s) and (t, s) are the temperature and entropy of the neutral Kaluza-Klein black 
hole and non-extremal brane respectively. 

5.1.3 Phases of non-extremal branes 

We can now apply the above boost/U-duality map to the known phases of Kaluza-Klein 
black holes reviewed in Section [3j We restrict ourselves here to the lower part of the 
phase diagram n < l/(d — 2), i.e. the uniform black string, the non-uniform black string 
and localized black holes. Recall at this point that all these phases have SO{d — 1) 
symmetry and have a metric that can be taken to be of the form of the ansatz (|3.1ip . 
Hence the corresponding non-extremal solutions are of the form (|5.7p - (|5.9p with the ansatz 
(|3.1ip substituted. The phases obtained from the bubble-black hole sequences [50] will be 
considered elsewhere [143J. 

Applying the map to the uniform string branch we recover the known phase of non- 
and near-extremal branes smeared on a circle, which we call the uniform phase. For the 
sake of clarity we write the solution explicitly here 



ds' 



d-2 



-fdt 2 + J^(A/) 2 + H 

8=1 



L 2 



27T 



-dR 2 + dv 2 + R 2 dn 2 d _ 2 



(5.12) 



e 2<t> _ R a ^ A _ CQth a _ A du l A . . . A du p 



(5.13) 



f = l R ° 



d-3 



d-3 



H = 1 + sinh 2 a ° „ . 

R d-3 R d-3 



(5.14) 
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In Section 15,21 we demonstrate that this solution suffers from a GL instability which we 
relate to the GL instability of the uniform black string. The thermodynamic quantities of 
the solution will be given in Section \6 . 2 1 when we discuss thermodynamic stability and the 
correlated stability conjecture. 

A strong indication of a possible instability can be deduced already by applying the 
map to the non-uniform black string branch from which we generate a new phase of non- 
extremal p-branes on a circle with a connected but not translationally invariant horizon 
around the circle, i.e. a horizon which is non-uniformly distributed along the circle. We 
denote this phase as the non-uniform p/mse@ 

Using the map (|5.10p from the neutral case to the non-extremal case we see that this 
new branch emerges out of the uniform phase (with n = l/(d — 2)) at the critical mass 

, (rf-1) b 2 c _ (d- 3) MGL 

2{d-2) 1 + ^1 + 62 2(d-2)q 

We can furthermore use the map (|5.10p on Eq. f)3. 15|) . This gives 

n{ji\q) = j— ^ - j(q)(fi - ft c ) , 0</i-/2 c <l, (5.16) 



7(g) = 7 



d-1 7/^gl b c (d-l)(d-3) + (d-2) 7 /x GL 



(5.17) 



2(d-2) 2{d-l) yiTbf 2(d-l)(d-2) 

with 7 given in Table [3l Eq. f|5. 16j) captures the behavior of the non-uniform phase 
of non-extremal branes on a circle for masses slightly higher than Jl = Jjl c . Note that 
i(q) ~^ 7 f° r Q ~~ * 0) as expected. We see from (|5.16p that the non-uniform phase starts 
at (p,,n) = (fie, — 2)) and continues with increasing p, and decreasing fi with a slope 
given by j(q). It is not difficult to show [57] from the above results that the uniform phase 
has higher entropy than the non-uniform phase for the same mass p,, just as in the neutral 
case. 

For those dimensions d where further numerical data are available on the neutral non- 
uniform branch, one can use the map to derive the corresponding data for the non-uniform 
phase of non-extremal branes. In this way one can obtain, for a given charge q, diagrams 
for any thermodynamic quantity one may wish to consider. 

In Section \5. 2 1 we show explicitly that p, c is a critical GL mass for the uniform phase of 
non-extremal branes on a circle. This is natural in view of the fact that the non-uniform 
phase implies the presence of a marginal mode at that mass. 

Finally, the map can be applied to localized black holes on a circle to generate a phase 
of non-extremal p-branes that are localized on a circle. This phase will be called the 
localized phase. Ref. [57] gives the background to first order in (fi — q), obtained from 



24 In [144] a new non-uniform phase of certain near-extremal branes on a circle was conjectured to exist 
for small energies. There does not seem to be any direct connection between the branch that we find here 
and the one in 144 . We comment further on 11441 in the conclusions in Section \8. 2 1 
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(|5.7p - (|5.9|) and the neutral solution of Ref. For the relative tension we have in this 
case 

fi(ft9 ) = ^A i( ,- 9 ) +0 (( P - 9 f), ^ % 2 %-J (5.18) 

In the (p,, n) phase diagram, the localized phase starts at the extremal point (q, 0) and 
goes upwards with a slope that can be read off from (|5.18p . The corresponding entropy 
and temperature of the branch can also be computed from the neutral thermodynamics. 
Again, for those dimensions d where numerical data are available for the neutral phase, 
one can use the map to obtain corresponding diagrams for any thermodynamic quantity 
of the non-extremal localized phase. 

5.1.4 Phases of near-extremal branes 

Given the boost/U-duality map from neutral and static Kaluza-Klein black holes to non- 
extremal branes on a circle, we can go one step further and take the near-extremal limit. In 
this way, we obtain a map that takes a Kaluza-Klein black hole solution to a near-extremal 
brane solution on a circle. 

The definition of the near-extremal limit of a non-extremal p-brane on a circle is 

g^oo, L^O, g = d __ 2 fixed , I = Lyfq fixed , x a fixed (5.19) 

where V p is the (spatial) volume of the p-brane, L the circumference of the circle, q 
the rescaled charge 25 ! and x a the transverse directions parameterizing x S 1 in the 

asymptotic region. This definition ensures that the energy above extremality /2 — q is finite 
and sensitive to the non-trivial physics associated to the compact direction. 

The near-extremal branes defined in this way can be characterized by two dimensionless 
quantities (e, r) which are respectively rescaled versions of the energy above extremality 
and the tension in the circle direction (see [57] for the precise definitions). Applying 
the limit (15. 19H to the general solution (|5.12[ )- (j5.14|) we obtain the following result for 
near-extremal branes on a circle 

ds 2 = i-Udt 2 + ^{du 1 ) 2 + HV ab dx a dx b \ , (5.20) 



e 



where 



2(f> = H a , A (p+1) = H~ l dt A du l A • • • A du p , (5.21) 
- t 1-U t l 2 (2Tr) d - 5 

H = h d —, h d = { > . 5.22 
c t {d - 3)i2 d _ 2 

The fields in (|5.20p - (|5.2ip have been written in units of L/(2n), i.e. we have rescaled the 

fields with the appropriate powers of L/(2n) to get a finite solution. As a trivial check, 



25 In terms of the charge Q and the ratio g, the rescaled charge is defined as q = gQ. 
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one can show that the near-extremal p-brane solutions generated in this way correctly 
asymptote to the near-horizon limit of the extremal p-brane on a transverse circle (15. 2j) . 
(|5.3p . which is taken as the reference space when calculating energy and tensions. 

To summarize, we have seen that there is a direct map from any Kaluza-Klein black 
hole in d + 1 dimensions to a near-extremal p-brane on a circle {D = d + p + 1). The 
corresponding map between thermodynamic quantities gives the near-extremal quantities 

Ml 

d + n (d — l)n % , rr 5 



/x , r = 2^- '— , t = tVts, $ = -=, (5.23) 



2(d-ly' d + n Vis 

where (fj,, n, t,s) and (e, r, t, s) are the (rescaled) mass/energy, relative tension, tempera- 
ture, entropy of the neutral Kaluza-Klein black hole and near-extremal brane respectively. 
The latter quantities are defined as 

e = gE, r=^, 1 = 1 f, a = |fi (5.24) 

where E, T, T, S are the energy, transverse tension, temperature, entropy of the near- 
extremal brane, and / and g appear in (|5.19p . With the use of f|5.23j) we can obtain the 
near-extremal (e, r) phase diagram from the (/x, re) phase diagram of neutral Kaluza-Klein 
black holes. 

In the above language the near-extremal Smarr formula and first law of thermodynam- 
ics take the form 

is = 2 d ~ 2 ~ r e , 5e = ids (5.25) 

As a consequence, knowing a curve r(e) in the near-extremal phase diagram determines 
the entire thermodynamics. 

We can apply the above map to all the known phases of Kaluza-Klein black holes as 
we did with the non-extremal branes on a circle . Here we restrict ourselves again to the 
uniform, non-uniform and localized phase. Note that since these solutions fall within the 
ansatz (|5.5p . we know that the corresponding near-extremal branes will be of the form 
(|5.20p - (|5.22p with the ansatz substituted. Hence we get the following three phases of 
near-extremal branes on a circle: 

• The uniform phase. In this phase a near-extremal brane is uniformly smeared on 
a transverse circle. The solution can be obtained with the map from the uniform 
black string and takes the well-known form 

6,-1 ( p 

ds 2 = -fdt 2 + ^(du*) 2 + H 



i=i 



jdR 2 + dv 2 + R 2 dn 2 d _ 2 



(5.26) 



e 



2<t> = H a , A (p+1) = H~ l dt A du 1 A • • • A dvP , (5.27) 



#n~ 3 - h d - / 2 (2vr) d - 5 
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Setting n = l/(d — 2) in (|5.23|) we find the relative tension r = 2/{d — 1). The 
thermodynamics of the uniform phase i, 



d-l 



4-7T 1 / 9. e \ 2(d-3) 

In Section 15.21 we demonstrate that the uniformly smeared near-extremal p-brane 
solution has a GL-instability. 

• The non-uniform phase. This is a phase with a configuration of near-extremal branes 
that are non-uniformly distributed around a circle. This phase is obtained by apply- 
ing the above near-extremal map to the non-uniform black string branch reviewed 
in Section [3l It emerges from the uniform phase given above at the critical energy 

d-l . 

VGL (5.31) 



2(d-2 

which follows from the map (|5.23p . The thermodynamics of this phase near the 
critical point is 

S(e) t (d-l) 2 7 



s u (e) 4d(d - 3) 2 e c 



(e-e c ) 2 + 0((e-e c ) 3 ) . (5.32) 



with e c the critical energy (|5.3ip . t c and s c the corresponding critical temperature 
and entropy, and c the (rescaled) specific heat 

d{d - iy5c (5.34) 



d(d - 5) + 2(d - l)7e c 



The numerical values e c , 7, t c , s c and c, of the energy, slope, temperature, entropy 
and specific heat at the critical point are listed in Table [4] [57 . 

The localized phase. This is a phase of a near-extremal brane localized on a trans- 
verse circle. It is obtained by applying the map to the black hole on cylinder branch. 
Since the black hole on cylinder branch starts at (p, n) = (0, 0), we deduce from the 
map (|5.23|) that the localized phase starts at (e,r) = (0,0). The vanishing of the 
tension along the transverse circle is expected since the brane becomes completely 
localized on the circle in the limit e — > 0. Furthermore we can use the analytical 



3 f is the rescaled free energy denned as f = e — is. 
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d 


4 


5 


6 


7 


8 


9 


f „ 


2.64 


1.54 


1.09 


0.71 


0.46 


0.31 


/ 


0.25 


0.39 


0.55 


0.88 


1.42 


2.33 


tc 


0.75 


1 


1.07 


1.05 


0.97 


0.89 


5c 


2.33 


1.54 


1.22 


0.91 


0.68 


0.53 


C 


-433 


6.46 


3.06 


1.77 


1.14 


0.79 



Table 4: In this table we list the critical energy e c and the constant 7 determining the 
non- uniform phase of near-extremal branes on a circle for < e — e c <C 1. In addition we 
list the critical temperature t c and the critical entropy s c , as well as the heat capacity c 
of the non-uniform phase at the critical point. These numbers are obtained in Ref. |57j . 



results of Ref. [39) to explicitly compute the first correction to the solution and ther- 
modynamics of the non- and near-extremal branes localized on a circle (see |57j for 
more details). The corresponding small black hole thermodynamics is 

™ = W« t *{i+ lj ^&+<>&)> <"«> 



1 



^M^Q)- ( , 36) 



and 



w rf) „ 2(d-2) / 2(d- l)C(d - 2)K\ d) ,.^z21 / .. 2(d-2) x\ 

f(t) = -< } t— f 1+ 1 ^Jh^_J_l t — +o((t— ) 2 )j , (5.37) 



2(d-2) 



*? - *iW* (^) " ■ (-38) 

One may also compute the second order corrections to these expressions using the 
second order correction in f|3. 17j) . We leave this as an exercise for the reader. 

From the analysis of the non-uniform phase we observe that, in parallel with the non- 
extremal case, we also have a map of the Gregory-Laflamme mass /xgl of the uniform 
black string branch to a critical energy e c . This is the energy where the non- uniform 
phase connects to the uniform phase of near-extremal branes on a circle. The existence 
of the non-uniform phase suggests that near-extremal branes smeared on a circle have a 
critical energy below which they are classically unstable. This is indeed the case, as we 
show in Section [5.2.21 

5.2 GL-instability of smeared non- and near-extremal branes 

In the previous subsection we reviewed the boost/U-duality transformation that takes a 
neutral black string to a Dp-brane smeared on a transverse direction. Following Ref. [66] 
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we describe in this subsection how to use this transformation to map the unstable GL 
mode of the neutral black string to the unstable mode of the non-extremal and near- 
extremal smeared Dp-branes. As we take the extremal limit the unstable GL mode con- 
nects smoothly to the marginal modes of the extremal smeared brane. 

5.2.1 Non-extremal smeared Dp-branes 

The boost/U-duality transformation of the previous subsection converts the uniform neu- 
tral black string solution (|2.9p to the smeared Dp-brane solution 



ds 2 = H- 1 / 2 



-fdt 2 + ^2dx 2 



+ H 1 / 2 [f- l dr 2 +dz 2 +r 2 dn 2 __ p ] , (5.39) 

e 2<t> = H (3- P )/2 ( Aqi p = cotha( #-i _ !) 5 f = 1 _ r Al ) H = l + sinh 2 a(l - /) . 

This is the same solution as (j5. 12|) - ()5. 14j) now written in the string frame without the 
■k- rescaling i.e. with dimensionful radial coordinate r and S 1 coordinate z ~ z + L. 

Z7T ° 

Applying the transformation to the neutral GL mode (|2.11j) . (j2. 12 j) will give (before the 
type II reduction) an eleven dimensional background with a non-normalizable exponential 
dependence on the M-theory circle direction y (the critical case k = k c is an exception). It 
was pointed out in [SB] that this problem can be avoided by considering a complex metric 
perturbation and finding a complex transformation that i) gets rid of the exponential y 
dependence and ii) still has the same effect on the zeroth order part of the metric, i.e. 
on the neutral black string metric. In the end one is instructed to take the real part of 
the transformed perturbation. This last step is allowed because the perturbed Einstein 
equations of motion are linear. 

The precise form of the transformation is as follows. First define the coordinates y' 
and z by making a complex rotation of the y and z coordinates 

(coshw — isinhtt;\ 
(5.40) 
i sinh w cosh w I 

Recall that the y and z directions correspond to the M-theory circle and transverse circle 
respectively. Notice that the neutral black string metric (|2.9p embedded in eleven dimen- 
sional gravity is invariant under (|5.40|) since (dy') 2 + (dz) 2 = dy 2 + dz 2 . At the next step 
we supplement the transformation (I5.40p with a boost that takes the (t, y 1 ) coordinates to 
the (t, y) coordinates given by 

A ( | A — | cos ' ia snina | ^ 4-q 

V y' J V sinh a cosh a J 

We choose the rotation parameters so that 

n n . 

sinhw = — tanha , k = k cosh w , £1 = — - — . (5.42) 
k cosh a 
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Then we can write the real part of the boosted perturbation of the eleven dimensional 
metric (obtained as described around (15.6H ) as 

V = Re | ex P + i^j j • ( 5 - 43 ) 

Using the relation P = A _1 M - 1 P(M~ 1 ) T (hr 1 ) T with P given in (pJ2|) we find 

°yy = -/V' sinh 2 a , P ty 



Pit = — ftp cosh 2 a , -P«y = —/-(/> sinh 2 a , Pt y = fij) sinh a cosh a 



P tr = 77 cosh a , P yr = -77 sinh a , P rr = f 1 x , P sp here = r 2 K 



(5.44) 



where we made the relabelings t t and y — > y. Notice that the complex rotation 
(|5.40p precisely ensures that the exponential factor does not have a y-dependence. We 
can therefore apply the same U-duality transformations on (|5.43p - (|5.44p as on the boosted 
neutral black string. After these U-duality transformations, consisting of one S-duality 
and p T-dualities, we conclude that the perturbed non-extremal Dp-brane smeared on a 
transverse direction can be written as I66r^l 

ds 2 = H- 1 ' 2 
+H}' 2 



v 

—f c dt 2 + dx 2 + 2r] cosh aSdtdr 
»=l 

f-\l + x£)dr 2 + dz 2 + r 2 (l + k£ )dQ^_ 



e 20 = ff (3-p)/2 ) Aoi p = cotha (#-i _ !) ( A rl ... p = -H~ l rj sinh a£ , 



6— p 
-6-p 



/ = 1 -J^> ff = l + sinh 2 a(l-/) , / c = /(l + , # c = l + sinh 2 a(l-/ c ) , 

£ = cos ffcrjj" 1 ;^ exp ffir^Hj . 

(5.45) 

The perturbed Dp-brane background appears here in one compact expression. From this 
one can easily deduce the explicit form of the perturbation by expanding to first order. 
The unstable mode appearing in [56] is equivalent to (|5.45p for p = 0, though written in 
a different gauge. 

Note that the functions t/j, 77, x an d « are still solutions to Eqs. (|A.5[) - (jA.8|) with / 
given by (|A.2p . In particular tp is a solution of Eq. (|2,13p . All these functions depend on 
the variable x = rfc/rn just as in the case of the neutral black string perturbation. 

Furthermore, we see from (|5.42|) that k 2 = k 2 + Q 2 tanh 2 a and £l = Q/ cosh a. Hence, 
we can obtain as a function of k by using the functional dependence £l(k) for the 
neutral black string (as sketched for p = 0, ...,5 in Figure [H where p = 9 — d). We 
note that the critical point (k,Cl) = (k c ,0) is mapped to the critical point (k,Q) = (k c ,0), 
which corresponds to the marginal mode of the Dp-brane with a transverse direction. This 



27 As usual, in the special case of D3-branes the gauge field strength is self dual, so that we have 



F 5 = (dAt+*dA*)/V2. 
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marginal mode appears at the origin of the non-uniform phase of non-extremal Dp-branes 
with a transverse direction [57] . From Figure [T] we also deduce that k < k c for any k < k c 
and the function Ct(k) never exhibits time-dependent modes with wavelength smaller than 
the critical one. 



5.2.2 Near-extremal smeared Dp-branes 

It is straightforward to take the near-extremal limit f)5 . 19j) of the perturbed non-extremal 
solution (|5.45p . In order to make the connection with gauge theory quantities more trans- 
parent later on we find it convenient to define here the 't Hooft and Yang- Mills couplings 

as 

A = 5ym^, 5ym = (27rf- 1 ^r 2 , (5.46) 

where g s is the T-dual string coupling. It can be expressed in terms of the circumference L 
of the T-dual circle (on which the D(p + l)-brane is wrapped) with the use of the T-duality 
relations 

LL = (2ttI s ) 2 , g s = g s ^ . (5.47) 
In these parameters it is convenient to re-express the near-extremal limit f|5. 19[) as 

T Z 

l s , u = j2 > ^ = 72' 9ym , L fixed . (5.48) 

Respectively, the uniform near-extremal smeared Dp-brane solutions take the form 

p 

l~ 2 ds 2 = H^i-fdt 2 + dxj) + ff 1/2 (/ _1 du 2 + dz 2 + u 2 dn 2 7 _ p ) , 

i=\ _ ' (5.49) 

e^ = H^, A Q1 .., = H-\ H = JL, / = 1"^, 
where K is defined as 



r 6 ~ p cosh a sinh a = Klf A ~^ , K = A ( 27r ) ? 2p . ( 5 . 50) 

(6 - p)tt 7 _ p 



K is related to the near-extremal parameters I and in (|5,28|) via the relation 

/ L\ 4 ~ p I 2 

Applying the limit (|5.48p to the background (|5.45p we obtain the following perturbation 
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of the near-extremal Dp-brane (|5,49|) |66| 



l~ 2 ds 2 = H- 1 / 2 



-fcdt 2 + dx i + 2r]VKu P /2 ~ 3 £dtdu 



i=i 



+ x£)du 2 + dz 2 + u 2 (l + K£)dn 



7-p 



p 2<t> _ M3-p)/2 a _ fr-1 A 



-U ' 



-VK- l l 2 ul /2 -\S 



(5.52) 



/ = 1 



6-p 

l _0 

,6-p 



/ c = /(l + ^) , ff c 



K 



u" 



u 



6— p 



It, 



6-p 



l v>£ 



f - cos f x + Q 2 -^ ) exp ( u Q ^ P K- l l 2 £lt 



We note that the near-extremal limit (|5.48p keeps k and $7 fixed, and moreover keeps tfi, 
rj, x and At fixed as functions of the variable x which is now given as 



uk 
u 



(5.53) 



Therefore, the functions tp, r], x an d k are still solutions of the Eqs. (|A.5|) - (|A.8|) with / 
given by (IA.2jl . In particular tp is still a solution of Eq. (|2.13|) . 

The above map shows that the near-extremal limit of the time-dependent mode of 
smeared branes is well-defined and hence that near-extremal smeared Dp-branes are also 
classically unstable. 

5.2.3 Relation to the marginal modes of extremal smeared branes 

In general, the solution of extremal Dp-branes distributed along a single flat direction z is 
given as 



ds 2 



#-1/2 _ d f + dx 2 + H l ' 2 (dr 2 + dz 2 + r 2 dn 2 _ p ) , 



i=l 



(5.54) 



= 2</> 



H^' 2 , A 01 ... p = H^ 1 — 1 , 



where the harmonic function H (r, z) obeys the differential equation 



d 2 + ^—^dr + d 2 ) H(r, z)=0 



(5.55) 



away from the source distribution. With appropriate boundary conditions the general 
solution is 



H(r,z) = l+ / 
J— t 



dz'- 



p(z') 



( r 2 + ( z _ z /)2)(7-p)/2 ' ( 5 - 56 ) 

where the charge distribution function p(z) is arbitrary. The extremal uniformly smeared 
Dp-brane case corresponds to a constant density p. However, it is clearly possible to add 



57 



to the constant charge distribution a single mode of any wave-number q, in which case the 
harmonic function takes the form 

Kl 8 s ~ 2p 

H = H — h m(qr) cos(qz) . (5.57) 

K is the parameter we defined in (|5.50p and the function m{qr) solves the differential 
equation 

m"(y) + — —m'(y) - m(y) = . (5.58) 

y 

The solutions of Eq. (|5.58p are of the form m(y) oc y^( 6 ~p)/ 2 ' } JC^_ p y 2 {y) > &s(y) being the 
modified Bessel function of the second kind. The modes given by (|5.57p are marginal, 
time-independent modes. Hence, we recover the well-known statement that the extremal 
smeared branes exhibit marginal modes of any wavelength. 

It is interesting to consider what happens to such marginal modes if one perturbs 
the extremal brane to a non-zero temperature. Prom the above study of non-extremal 
Dp-branes distributed on transverse directions one does not expect the existence of static 
solutions for arbitrary charge distribution p(z). In particular, it is worth exploring whether 
there is any relation between the extremal brane modes (|5.57p and the GL modes of the 
near-extremal branes. For that purpose, we need to consider briefly the near-horizon limit 
(I5.48P of the extremal Dp-brane background given by (|5.54p and (|5.57p . One easily finds 



l- 2 ds 2 = H- 1 / 2 \-dt 2 + V dx 2 ) + H 1 ' 2 {du 2 + dz 2 + u 2 dn 2 _ p ) , 

V i=i / (5.59) 

e 2<P = £(3-p)/2 ; Aqi p = h-\ H=J£-+ m (qu) cos(qz) . 

Note that the function m(qu) has been rescaled appropriately, and that m(y) still obeys 
Eq. (15381) . 

Consider now the extremal limit of the perturbed near-extremal solution given by 
(I5.52p . The extremal limit corresponds to uq — > with u and K kept fixed. Finiteness 
of the argument of the cosine factor implies that we also need to keep Vk 2 + Q 2 /uq finite 
in the limit. In terms of the dispersion diagram of Figure [1] we move closer and closer to 
the point (k, £1) = (0, 0) on the left part of the curve. In the process, the time-dependent 
exponential factor in the function £ in (|5.52p disappears. A closer look at the solution 
and the equations of motion of appendix [X] reveals the scaling [66] 



k ip(x) r){x) x( x ) k{x) 
uq -> with u , K , -g^ , , , kept fixed . (5.60) 



""0 ""0 U Q ^ Uq F Uq u,q 



Applying this limit to the perturbed near-extremal solution (|5.52p . we get the following 
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background 



l~ 2 ds 2 = HZ 1 ' 2 



dt 2 + ^2 dx 1 



i=l 



+ 



Hi 12 



du 2 + d'z 2 + u 2 dQ 2 _ p 



e 20 = £(3- P )/2 A 



01. -p — H c , H c 



K 

u 6 ~p 



■£ , 



(5.61) 



£ = cos z 



We see that out of the four functions ip, r/, x and k that make up the near-extremal GL 
mode only ip survives. Moreover, the equation (|2.13p for ip(x) becomes 



iP"(x) + 7 —ip'(x) 



1 + 



k 2 



ip(x) 



. 



(5.62) 



We deduce that the background (|5.6ip corresponds precisely to the background ([5.590 
describing an extremal smeared Dp-brane perturbed by a marginal mode with the wave- 
number 

Vk 2 + n 2 



(5.63) 



In particular, one can check explicitly that the ip equation (j5.62p turns into Eq. (|5.58[) 
by identifying m(qu) = ip(ku/uo). This means that the extremal limit of the GL mode 
becomes stable as also observed for unsmeared black branes in Ref. [131J. 

One can get marginal modes of any wave-number in this way. For small k we have 
= 7&, where 7 is an appropriate number. Hence, q = yl-f ^k/uo, and by tuning 
k/uo we can get any value of q we want. To summarize the picture, we see that the 
GL modes of near-extremal smeared branes become marginal modes of extremal smeared 
branes in the extremal limit. To put it differently, when we deform an extremal smeared 
brane to make it non-BPS, the marginal modes of the supersymmetric configuration turn 
into near-extremal GL modes (|5.52p . In general, one can imagine extremal branes with a 
charge distribution of the form 



P = PO + ^2 Pn COS (lri^ 



(5.64) 



where p n are small compared to po. In the presence of a non-BPS perturbation the solution 
turns into a sum over GL modes of the form 



3-1 
7^0 



--q n t 



(5.65) 



with uq small. In this way we can find the non-BPS continuation of any extremal pertur- 
bation of the extremal smeared branes, since we can Fourier decompose the perturbation 
and put it into the form (|5.64p . 
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5.3 GL-instability of D-brane bound states 

So far we have considered only the case of singly-charged D-brane configurations that 
arise as thermal excitations of half BPS branes in type IIA/B string theory and M-theory. 
Another interesting class of examples in type II string theory and M-theory are non- 
extremal D-brane bound states. In this subsection we review known facts about the 
stability properties of a class of D-brane bound states in type II theory and mention some 
interesting possible generalizations. 

General spinning D-brane bound state solutions of type II supergravity appear in [145J, 
|146j . [147] . For completeness, we summarize here the general solution and then proceed to 
discuss some special cases in more detail. The general solution represents a spinning bound 
state of D(j> — 2fc)-branes with k = 0, m and non-zero S-field. The background depends 
on the non-extremality parameter tq, the charge parameter a, the angular momenta l{ 
(i = 1, 2, n = and the angles 6k, k = 1, 2, m. The metric in the string frame is 

(m p \ 

k=l i=2m+l / 

+H 1 ' 2 (f-'Kg^dr 2 + A a0 dr, a dr,^ 



+H~ 1/2 -^ 1 ^^ I E litjrfiUjdtketyj - 2 cosh a ^ litfidtdfo ] (5.66) 

'" \i,j=l i=l 



with a dilaton of the form 

m 

e 2 * = ff( 3 -ri/ 2 JjD fc . (5.67) 

k=l 

The NSNS -B-field has rank 2m < p and is given by the relation 

B 2k -i,2k = tan6 k (H- 1 D k -l) , k = l,...,m. (5.68) 

The non-zero RR gauge potentials, whose explicit form can be found in Appendix B of 
[146] . are A p _ 2 k+i with k = 0, ...,m. 

The functions entering in this background are 

1 rl _p sinh 2 a , , N _i , 

H = 1 + 7tF 7 p » Afc = sin 2 9kH + cos 2 fc , 5.69 

1 r n 7 - p , 1 r, 7 " p -5- ' " 2 



with 



f=i -w,fri> /=i -i^- ^=n( 1+ ^) (5J0) 



i=l 



W p = K 9 „ p L 9 _ p . (5.71) 
Kg-p and A Q( g are functions appearing in the flat transverse space metric 

9-p 

^2(dx a ) 2 = K 9 - P dr 2 + Aapdrfdrf . (5.72) 

0=1 
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Explicit expressions for Kg_ p , A Q( g and ^ can be found in Appendix B of Ref. |145j . 
It should be pointed out that the following condition holds as a consequence of charge 
quantization of the Dp-brane 

7—p i ■ i (27r) 7 -PNg s ll- p " j 

r n ^ cosh a sinh a = — — — I (cos 6^ ) . (5.73) 

(7-p)0 8 _ p AJ. 

5.3.1 The D(p - 2)-Dp system 

A direct analysis of the classical stability properties of the non-rotating D0-D2 bound 
state has been performed by Gubser in [59]. The non-extremal D0-D2 solution is a special 
case of the above general background, whose explicit form is 

ds 2 = H-^i-fdt 2 + Dddx 1 ) 2 + (dx 2 ) 2 )) + H 1 ' 2 (^dr 2 + r 2 dn 2 

A\ = cotha sin# ^1 — — \ dt , 

A 3 = coth a sec 6 ^1 - ^\ dt A dx 1 A dx 2 , (5.74) 



e 20 



1/2 L> , 5 2 = tan# ( 1 - ) dx 1 A cix 2 



sinh 2 a 1 

r s sm z 9 + cos 2 r° 

The thermodynamic quantities of the solution will be given in Section 16.21 when we discuss 
thermodynamic stability and the correlated stability conjecture. 

The above static solution describes a configuration of uniformly smeared DO charge 
inside the worldvolume of a D2-brane. As in the case of non-extremal D-branes smeared 
uniformly on a transverse circle one might expect that this configuration exhibits a GL 
instability towards a non-uniform redistribution of the DO charge. Gubser analyzed an 
important aspect of the perturbative stability of (|5.74p by looking for an s-wave marginal 
threshold mode that depends only on x 1 and r. Reducing the full ten-dimensional action 
on the 5 6 , t and x 2 directions he obtained a complicated two dimensional system whose 
dynamics he was able to analyze with a combination of analytical and numerical tech- 
niques. The results exhibit the presence of a boundary of stability in the phase diagram 
°f (Qq/M,Q2/M) charges, which separates an unstable from a stable phase. Moreover, 
it was shown that the presence of a GL instability is in very good agreement with expec- 
tations based on a local thermodynamic stability and the correlated stability conjecture. 
The precise form of the boundary of stability and its connection with thermodynamics 
will be discussed in more detail in Section [6j 

Once a marginal threshold mode is obtained for the D0-D2 system a corresponding 
threshold mode can be deduced for the more general non-rotating D(p — 2)-Dp bound 
state with T-duality. It would be interesting to proceed further analyzing the full phase 
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diagram of the D(p — 2)-Dp bound state in analogy with what has been achieved for the 
singly-charged non-extremal Dp-branes on a transverse circle. 

5.3.2 Other bound states and further generalizations 

One can generalize the above discussion to other non-extremal bound states of D-branes. 
The classical stability of T)(p — 4)-Dj> and Fl-Dp bound states has been considered with 
general arguments based on boost/U-duality transformations in |60] and from the point 
of view of thermodynamics and the correlated stability conjecture in [60, 61J (see Section 
[6] for further details on this aspect). 

We have seen the presence of a GL instability in a variety of black systems. In general, 
the task of determining the existence of an instability and the full phase structure of the 
system requires a very complicated analysis in classical gravity coupled with a set of matter 
fields. The number of charges and parameters of the system has a non-trivial effect on 
the properties of the GL instability, which one would like to explore in as many different 
cases as possible. Some interesting generalizations include the above-mentioned bound 
states with rotation, black objects with multiple (more than two) charges, non-extremal 
intersecting branes [148J etc. 

6 Thermodynamics and the correlated stability conjecture 

In this section we review the connection between classical stability and thermodynamic 
stability of black branes in supergravity. This connection is known as the correlated 
stability conjecture, which we first present in some detail. As an illustration, we consider 
the conjecture in two examples, namely that of smeared Dp-branes and D-brane bound 
states, which were the subject of the previous section. Subsequently we discuss the support 
of the conjecture based on a general semiclassical argument. We end by commenting on 
known counter examples to the conjecture and how the conjecture may be refined or 
restricted to avoid violation. 

6.1 Connecting thermodynamic and classical stability 

We have seen that the Gregory-Laflamme instability is a classical instability of black brane 
solutions in gravity where the mass tends to clump together non-uniformly. It is natural 
to ask if there is a simple intuitive understanding for the existence of such instabilities. 
Using the entropy to say something about the dynamical properties of classical gravity is 
natural from the point of view of the second law of thermodynamics. Already in [2], a 
global thermodynamic argument was presented for the instability, namely that for small 
masses the entropy of an array of black holes is higher than that of the black brane 
with the same mass. This argument suggests that the black branes might be unstable. 
Nevertheless, a global thermodynamic argument like the one above can sometimes be 
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misleading. Given the existence of a higher entropy configuration the second law permits 
a classical instability, but does not require it. For a black string, in particular, this means 
that we cannot deduce the endpoint of the Gregory-Laflamme instability solely on the 
basis of global thermodynamic arguments like the above entropy argument @ It is clear 
that one should look for a principle that relates two local properties on the configuration 
space of gravity. 

As a refinement of the entropy argument Gubser and Mitra proposed in \62\ [63] a 
natural conjecture that relates classical and local thermodynamic stability. The conjecture, 
which from now on will be referred to as the Correlated Stability Conjecture (CSC), 
postulates that: 

Gravitational systems with translational symmetry and infinite extent exhibit a Gregory- 
Laflamme instability if and only if they have a local thermodynamic instability. 

Local thermodynamic stability is defined here as positive-definiteness of the Hessian matrix 
of second derivatives of the mass with respect to the entropy and any other charges that 
can be redistributed over the direction in which a Gregory-Laflamme instability can in 
principle occur. The additional assumption of translational symmetry and infinite extent 
has been added to insure that finite size effects do not spoil the thermodynamic nature of 
the argument and to exclude well-known cases like the Schwarzschild black hole which is 
classically stable yet has negative specific heat and is therefore locally thermodynamically 
unstable. 

The CSC provides an exciting link between perturbative classical stability and local 
thermodynamic stability bypassing the obvious problems with using global entropic argu- 
ments. The importance of such a connection on a practical level is obvious given the usual 
complexity of the classical stability analysis. In what follows, we outline the main aspects 
of the CSC, review its validity in a series of examples and discuss its limitations, known 
counterexamples and previous attempts to prove it. 

For a system with n conserved charges Qi the criterion of local thermodynamic stability 
translates (after making the appropriate choice of ensemble) into a positivity criterion 
for the (n + 1) x (n + 1) Hessian matrix Hm = ( gg ) , which involves the second 
derivatives of the mass M with respect to the entropy S and the charges Qi. The charges 
Qi entering this definition are determined by the choice of thermodynamic ensemble, 
which is crucial, because it affects the discussion of local thermodynamic stability and 
ultimately the validity of the conjecture as such. The usual practice is to consider gravity 
in the canonical ensemble, where the temperature T is kept fixed and the partition function 
of the system is expressed as a function of T. With more degrees of freedom the choice 

28 In fact, the viewpoint that an unstable uniform black string decays to a localized black hole has been 
challenged in [T7] where it was argued that the black string horizon cannot pinch off in finite "horizon 
time" and that the endpoint of the Gregory-Laflamme decay should be some type of non-uniform black 
string. See Sections 12.21 13.2.71 for a discussion of the endpoint of the GL instability and Section 18.21 for 
some further comments. 
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of ensemble for a general system can become less obvious. To illuminate this point it 
is instructive to consider the class of smeared or unsmeared, magnetically or electrically 
charged black brane solutions of the previous sections. 

Consider then a general (non-extremal) black Dp-brane background of type II string 
theory. Let {zi} denote the set of non-compact directions along which the brane back- 
ground exhibits translational symmetry. The brane may be charged electrically or mag- 
netically under a (p + l)-form gauge field, with a corresponding charge Q p . We say that 
an electrically charged brane is smeared along a direction Zi when it is not charged along 
this direction. For magnetically charged branes, we follow the opposite convention and 
say that the brane is smeared along Zi, whenever it is charged along z,-P^| 

We can now formulate for the CSC in a precise manner which thermodynamical en- 
semble we should work in. If a brane is not smeared in any of the non-compact directions 
with respect to a charge Q p then we should work in the canonical ensemble with respect 
to that charge, i.e. we should keep Q p fixed. On the other hand, if a brane is smeared in a 
particular direction Z{ with respect to the charge Q p we should work in the grand canonical 
ensemble with respect to that charge, i.e. we should instead keep the chemical potential 
v p corresponding to Q p fixed. This rule for the ensemble used for CSC applies equally well 
to electrically or magnetically charged branes as expected by any sensible formulation of 
the CSC that is invariant under the electric-magnetic duality. 

As a simple illustration, consider a non-extremal DO-brane solution smeared on a 
transverse direction z in type IIA theory. The corresponding supergravity solution in the 
string frame takes the form 

ds 2 = -H- 1 l 2 fdt 2 + HV 2 {f- 1 dr 2 + r 2 dn 2 7 + dz 2 ) , e 2 * = H 3 ' 2 , 

6 (6.1) 
A = coth a(H~ l - 1) , f(r) = 1 - -§ , H(r) = 1 + sinh 2 a(l - /) . 

In this example the DO-brane charge Q is smeared along the z-direction and can be redis- 
tributed there freely. Hence, in the context of the CSC it is appropriate to consider local 
thermodynamic stability in the grand canonical ensemble, where Q is allowed to vary. 

With a T-duality transformation along z the background (|6.ip turns into the Dl-brane 
solution 

ds 2 = H- 1 ' 2 (-fdt 2 + dz 2 ) + H 1 ' 2 {r l dr 2 + r 2 dfl 2 7 ) , (6.2) 
e 24> = H , A 0z = cothaiH- 1 - 1) . 

The Dl-brane is now charged along the z-direction and the corresponding charge is a 
fixed quantity. Accordingly, we should now consider local thermodynamic stability in the 
canonical ensemble. The application of the CSC with this choice of thermodynamics yields 
a picture of stability that fits very nicely with the classical stability properties of these 
solutions as we see explicitly below. 

29 This definition is also very natural for D3-branes which are self-dual under electric-magnetic duality. 
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As a final comment, it is well-known in thermodynamics that one can move back 
and forth between the canonical and grand canonical ensembles with a simple Legendre 
transform. The above example suggests that within the context of the CSC it is natural to 
associate a Legendre transform with a T-duality transformation in supergravity. It should 
be noted, however, that the appearance of the Legendre transform is an essential feature 
of supergravity that treats momentum and winding modes asymmetrically. In the full 
string theory, where momentum and winding are exchanged by T-duality the Legendre 
transform would be unnecessary. For instance, a momentum instability for a smeared 
brane transforms under T-duality into a winding instability for a wrapping brane which 
is invisible in supergravity. 

6.2 Examples 

The first evidence for the validity of the CSC was given in the original papers [62\ [63] 
from the perspective of the AdS/CFT correspondence. These papers studied a set of 
Ad^-Reissner-Nordstrom black hole solutions in J\f = 8 supergravity in the large mass 
limit where the black holes are expected to exhibit similar behavior to black branes. Up to 
a small discrepancy arising from a numerical error they found that a classical instability 
appears precisely when a local thermodynamic instability sets in. 

Further progress was made in the paper by Reall |64j where a semiclassical proof was 
found of the CSC (see Section I6.3P in the case of singly-charged black branes. In this 
case the local thermodynamics in the canonical ensemble is stable for large charges and 
unstable for small charges. This fits precisely with the results of the numerical study of 
the classical stability for these branes [3l H5I1 H551 fTM H3S1 H3Ij . 

Below we consider instead another class of brane backgrounds, in which the charges 
can be smeared in one or more directions. 

6.2.1 Non-extremal smeared black branes 

In Section [5] we saw that the non-extremal smeared Dp-branes exhibit a Gregory-Laflamme 
instability below a critical mass. According to the CSC a corresponding instability should 
be present in thermodynamics in the grand canonical ensemble [BB]. In verifying this 
statement it is instructive to consider both the canonical and grand canonical ensembles. 

In the canonical ensemble, the temperature and charge are kept fixed and the appro- 
priate thermodynamic potential is the Helmholtz free energy 

F(T, Q)=M-TS , dF = -SdT + udQ . (6.3) 

The condition for thermodynamic stability requires that the specific heat Cq is positive, 
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On the other hand, in the grand canonical ensemble the temperature and chemical 
potential are kept fixed so the appropriate thermodynamic potential is the Gibbs free 
energy 

G{T,v) = M -TS - vQ , dG = -SdT - Qdv . (6.5) 
The condition for thermodynamic stability is in this case 

dM\ ( dv\ 



Cqs \W) q >0 - cs {oq) t >0 - <66) 

where Cq is again the specific heat and c is the inverse isothermal electric permittivity. 
The condition (|6,6p follows from demanding that the Hessian Hq of the Gibbs free energy 
is negative definite. Because of the matrix relation Hq = —H^, where Hm is the Hessian 
of the mass M(S,Q), this is equivalent to demanding that Hm is positive definite, as 
required by local thermodynamic equilibrium. 

The thermodynamics of the non-extremal Dp-branes (|5.12p - (|5.14p smeared on a trans- 
verse circle of circumference L is given by 

T= IftrG (2^-P R ° [7"P+(6-p) S mh a], - = - ^^ R, cosh a , 

Q n 7 _ p V p L e ~P 6 6-p 

L= 16vrG (2vr)6-^o (6 - p) smh a cosh a , = 2LRq cosh a ' ^ = tanh «' 

(6.7) 

where V p is the world- volume of the brane and the extensive quantities, M/L, Q/L and 
5/ L correspond respectively to the mass density, charge density and entropy density along 
the transverse z direction. T is the temperature and v the chemical potential. These 
quantities satisfy the usual first law of thermodynamics dM = TdS + vdQ. 

With the use of the definitions (|6.6p and the explicit thermodynamic quantities above, 
we find for the non-extremal smeared braneJ^I 



7 — p + (8 — p) sinh 2 a 
— 1 + (4 — p) sinh 2 a 



cosh 2 a[l — (4 — p) sinh 2 a] TS 



(6.8) 



We see for p < 3 that the requirement of positive specific heat translates to 



Cq > 44> a > arcsinh(l/ v /4 - p) , (6.9) 

showing that there is a lower bound on the charge of the branes in order to be thermody- 
namically stable in the canonical ensemble. On the other hand, we have 



c>0 44> < a < arcsinh(l/- v /4 - p) , (6.10) 

which is incommensurate with the condition in (|6.9p . For p = 4,5 we have Cq < 
and c > 0. Hence for any p = 0,1, ...,5 it is impossible to satisfy the requirement 



30 These expressions are most easily computed using the identity that §f^fyU(r,s) = a(r's) l a(r's) ] ^ 
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of thermodynamic stability in the grand canonical ensemble. The presence of a local 
thermodynamic instability in the grand canonical ensemble verifies the validity of the 
CSC in this case and fits nicely with the general criteria for the choice of thermodynamic 
ensemble presented in the previous subsection. Note, however, that it is not possible to 
use the CSC to determine the exact value of the Gregory-Laflamme threshold unstable 
mode. 

It is also worth pointing out that the above thermodynamic quantities are invariant 
under T-duality. For example, the thermodynamics of a Dp-brane smeared on a transverse 
circle is identical to that of a D(p + l)-brane wrapped around the T-dual circle. The 
wrapped non-extremal D(p + l)-brane solution (for 1 < p + 1 < 4) is known to be stable 
in supergravity, which is reflected nicely in the thermodynamics if we choose to consider 
the CSC in the canonical ensemble. 



6.2.2 The near extremal limit 

The near-extremal limit of smeared Dp-branes presents an interesting subtlety. Since 
we send the charge to infinity, it seems that in the near-extremal limit there is only one 
relevant ensemble, namely the canonical ensemble. The corresponding Helmholtz potential 
is given by 

F(T) = E — TS , dF = -SdT, (6.11) 

and thermodynamic stability requires positivity of the specific heat C = dE{T) / dT . 

The thermodynamic quantities of near-extremal Dp-branes (I5.49|) smeared on a trans- 
verse circle are given by 31 l 

1 6- p 8-p 4vr 6 _p 1/2 Q-p(ul~ p \ 1 

E= g U ° S=-uoK K) , r = — ^— j . (6.12) 

Here, E is the energy above extremality defined by E = lim(M — Q) in the near-extremal 
limit (|5.48[) . whereas K (defined in (|5.5(jp ) and Q given by 

^£zE = ^C 2(6 - P) , \ - (2nf^LV p n 7 ^ , (6.13) 

are both finite in the limit. The quantities in (16.120 satisfy the usual first law of thermo- 
dynamics dE = TdS. 

From (|6.12p one obtains for the specific heat the simple relational 

C=^S. (6.14) 
4 — p 

Hence, for any non-zero temperature, the specific heat is positive for near-extremal smeared 
Dp-branes when p < 3, which is a well-known result. Since the charge goes to infinity in 



31 In terms of dimensionless quantities this thermodynamics was given in (|5.29p . 

32 An alternative way to obtain the same result is by taking the near-extremal limit (a — > oo) of the 
expression f|6.8|) for Cq. 
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the near-extremal limit and the chemical potential goes to one, it seems that these quan- 
tities cannot be varied anymore. This suggests that it does not make sense to consider 
the grand canonical ensemble for near-extremal branes. 

To put it differently, one could consider what happens to the quantity c in (|6.8p when 
we take the near-extremal limit. As a consequence of the limit a. — * oo we deduce that 
c — > 0. This would seem to imply that, infinitesimally close to the near-extremal limit, non- 
extremal branes are marginally thermodynamically stable in the grand canonical ensemble 
in clear contradiction with the classical stability analysis of Section 15.21 where we saw 
explicitly the presence of a Gregory-Laflamme instability in the near-extremal limit. It 
would appear that the near-extremal limit is an example where the CSC is manifestly 
violated. 

Fortunately, there is a simple resolution to this discrepancy [66j@ It turns out that 
in the near-extremal limit the right way to compute the quantity c in the grand canonical 
ensemble is not by using the non-extremal chemical potential i>, but by using a new chem- 
ical potential u, which can be understood as the "chemical potential above extremality" . 
This is in analogy to the fact that we do not compute Cq using the mass M (which goes 
to infinity in the near-extremal limit), but using the energy E above extremality, which is 
finite. 

To derive an expression for we note that the first law of thermodynamics for non- 
extremal branes takes the form 

dE = TdS + {v-l)dQ , (6.15) 

when written in terms of E = M — Q. Moreover, one can easily see from (|6.7|) . that the 
chemical potential v approaches 1 in the near-extremal limit a — * oo, and is not a free 
parameter anymore. It is therefore natural to define the rescaled chemical potential v and 
the corresponding rescaled charge Q as 

v = -2— lim - 1) , Q = s-ym Bm liQ , (6.16) 

where the factor of l s has been chosen such that i> is finite in the near-extremal limit. The 
finite factor (jy M has been inserted for later convenience. With the definitions ()6.16j) . we 
then find that (|6.15p can be written as 

dE = TdS + i>dQ, (6.17) 

which is a well-defined differential relation in the near-extremal limit. We have thus 
obtained a modified version of the first law of thermodynamics for near-extremal branes 
involving an extra term containing the rescaled chemical potential v and rescaled charge 
Q. 

33 Note that this also resolves the contradiction found in the numerical investigations of Ref. |140j . 
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Consequently, by using the non-extremal quantities (|6,7p and the near-extremal limit 
(|5.48p in the definitions ()6.16|) we find for near-extremal smeared Dp-branes the expressions 

K u 6 ~ p 
Q = (6-p)-, fi = —2—, (6.18) 

where K, Q are defined in (|5.50p . (|6.13p respectively. The near-extremal Gibbs free energy 
is then naturally defined by G(T, u) = E — TS — vQ, and the condition for thermody- 
namic stability of near-extremal smeared branes in the grand canonical ensemble follows 
immediately by analogy with (|6,6p . namely 

C6-(£) >0. »-(») >0. (6.19) 



- \STJq ' KdQJr 

Computing these quantities for near-extremal smeared Dp-branes, we find that the quan- 
tity Cq is identical to the one computed in (|6.14j) . On the other hand, for c we find using 
(IBTT8D and ([5^9j) thafl 

g u 6 ~p 

e = "4 M i^(4-p)(6-p) • ( 6 - 2 °) 

which is finite and negative for p < 3. Hence we conclude that near-extremal smeared Dp- 
branes with p < 3 are thermodynamically unstable in the near-extremal grand canonical 
ensemble defined above, in full agreement with the classical stability analysis and the CSC. 

6.2.3 The D(p - 2)-Dp system 

Another interesting class of examples in type II string theory involves systems of D-branes, 
for example D-brane bound states, which preserve various amounts of supersymmetry in 
the extremal case. For instance, consider the bound state of a D(p — 2)- and a Dp-brane, 
where the D(p — 2)-brane is embedded in the Dp-brane and smeared along its two trans- 
verse directions. In the extremal case, this system is a half BPS state. The corresponding 
supergravity solution appears in Section T5.31 The classical stability properties for the D0- 
D2 bound state were discussed in [SUJ [BD], where it was found that the system exhibits a 
Gregory-Laflamme instability in accordance with the local thermodynamic stability analy- 
sis and the CSC. Here we review the stability pattern that emerges from thermodynamics. 

The thermodynamic quantities of the D(p— 2)-Dp bound state can be found for instance 
in |146j . For zero angular momentum we have 

^=^gfe £ -0~ P (8-p+(7-p) S inh 2 a), T = , 
S = ^^^r^-P cosh a , Q= {7 ~ p £*£ 8 " p r 7 ~ p sinh a cosh a , v = tanh a , (6.21) 

z^sin# , Q p = QcosO , Q„-2 = Qs'mO . 



34 0ne can compute c directly from the near-extremal thermodynamics using c = gTp^ [ q^'k) > 

-8„-4 



Alternatively, one can show from (I6.16P that c = limZ s g YM c. We obtain the same result by using in this 
expression the non-extremal relation (|6.8p for c and then taking the near-extremal limit. 
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Figure 9: Sketch of the (Q p -2/M,Q p /M) boundary of stability for p = (red), p = 1 
(green), p = 2 (blue), p = 3 (yellow), p = 4 (pink). The bound state is unstable in the 
region below the boundary of stability. 



To determine the local thermodynamic stability of this system, we analyze the thermody- 
namics in a mixed canonical/grand canonical ensemble, where we allow Q p -2 to vary while 
keeping Q p fixed. In other words, we consider the smeared D(p — 2)-brane in the grand 
canonical ensemble and the Dp-brane in the canonical ensemble adhering to the general 
criteria of the previous subsection. 

Much like in the previous example of the non-extremal smeared branes one can show 
that the determinant of the Hessian matrix Hm = (^ dcfdqp ) (with q a 6 {S, Qp-2}) is 

d(T,v p -2) 16 G sech 4 a . , 9 2 /, -,\ 

det = m0^> = V ? al_ r rr" ((9 - P ) cosh' a - 1) ((5 " r) ^ "™ " ~ ^ ' 

(6.22) 

Hence, the condition of thermodynamic stability translates to 

det H M >0 ^ csch 2 a < (5 - p) cos 2 6 , (6.23) 

which can be satisfied only if p < 4. The boundary of stability, where det Hm = 0, can be 
written more compactly as 

^ 2 -2 + (6-p)^ 2 = l • (6.24) 
The stability pattern is summarized in Fig. [9j 

6.2.4 The D(p - 4)-Dp system 

Similar statements apply to the T)p-D(p — 4) bound states, which are quarter BPS in 
the extremal case. The supergravity solution for these states appears already in Section 
15.31 Again, one can apply the CSC to find the stability pattern for these states with a 
simple calculation in thermodynamics. As above, we are instructed to work in a mixed 
canonical/grand canonical ensemble, where we allow Q p -n to vary holding Q p fixed. 
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In a suitable parametrization (see [UT]) the thermodynamic variables of the D(p— 4)-Dp 
system take the form 

M = W I V(2 + (7 - p)(cosh2a p + cosh2a p _ 4 )) , T = ^ cosh 7 ~ p cosh ap _ 4 , 

S = Vp ^Q~ p Tq~ p cosh a p cosh a p _ 4 , v p = tanha p , v p -4 = tanha p _ 4 , (6.25) 

Qp = iJzJ $Mr= Z rl P sinh 2a p , Q p _ 4 = r 7 ~ p sinh 2a p _ 4 . 

Consequently, the determinant of the Hessian matrix becomes 

d{S,Q p -4) 



J^L r o P 8) sech 2 « p sech 4 g ip _4 ( _ _ 

K 2 f}2 cosh2a p _ 4 (7-p + 2cosh2a p ) + (7 - p) cosh 2a p v ' 



In this expression all the factors have a positive definite sign except for the last. 

There are three values of interest for p. p = 4 represents the D0-D4 bound state, for 
which we find an instability when 

cosh 2a 4 < 3 . (6.27) 

The existence of a corresponding threshold unstable mode in the classical stability analysis 
of this bound state was argued in [BD] . For p = 5 and p = 6 there is no stability boundary 
and the bound state is unstable all the way up to extremality. 

6.2.5 Other examples 

Other examples of bound states can be obtained from the above in the following manner. 
For instance, by S-dualizing the D1-D3 bound state we can obtain the F1-D3 bound state. 
This state exhibits the same thermodynamics as the D1-D3 system and therefore also the 
same stability pattern. We can also lift the D2-D4 bound state to the M2-M5 system. 
By further reducing the latter on a different circle we obtain the F1-D4 bound state. All 
three of these bound states have the same thermodynamics and the same stability pattern. 
Unfortunately, we cannot use T-duality to extend this list even further. As already pointed 
out, T-duality acts non-trivially on the CSC. On the level of the classical stability analysis 
the translation invariance is broken explicitly by the threshold unstable mode along the 
direction of the Gregory-Laflamme instability thus making it impossible to perform a T- 
duality transformation along these directions. On the level of thermodynamics, the same 
effect reappears in some cases as a change of the ensemble. 

One can also consider independently the general Fl-Dp bound state or three-charge 
bound states like the D0-D4-F1 system. The former was considered in [61] and the latter 
in [60]. 

Beyond these cases it would be interesting to consider the CSC in a wider class of 
examples, for instance in cases that exhibit rotation, such as spinning bound states. 
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6.3 Towards a proof of the conjecture and known counterexamples 

The correlated stability conjecture has proven to be a rather robust statement whose 
validity has been verified in a large class of examples. A semi-classical proof of the con- 
jecture for magnetically charged black brane solutions was given in [M] by Reall (see also 
|65j). Elements of this proof for more complicated systems, like that of electrically charged 
smeared black branes, were discussed in |61| and [66] (see also [149] for another approach 
to a proof based on holography). Recently [67] presented a set of counterexamples in 
solutions with scalar hair, which suggest that the conjecture cannot survive in its current 
form as a statement that applies to the most general case. To see what may go wrong, it 
is instructive to review the key elements of the existing arguments in favor of the CSC, 
the weak points of these arguments and the known counterexamples. 

To expedite the discussion let us denote again by {z 1 } the infinite extent directions 
along which a Gregory-Laflamme instability can occur and by x^ all the other space- 
time directions. For compactness, it is also convenient to denote the field content of the 
theory collectively by the fields ^ I . J is a multi-index label that summarizes all the 
metric components and any other fields including scalar fields (e.g. the dilaton) and gauge 
fields. The central object in the semi-classical argument of Reall [63] is the Euclidean 
path integral dimensionally reduced along z % . The Euclidean path integral can be written 
compactly as 




(6.28) 



where /[^• 7 ] is the Euclidean action and the path integral is taken over the fields Vl/ in 
the appropriate thermodynamic ensemble specified (for a class of examples) in Section [6.11 
For gravity this is always the canonical ensemble and the path integral is taken over all 
the smooth Riemannian geometries that are asymptotically S 1 x Ai. S 1 is the Euclidean 
time circle with circumference /3 = A and M the asymptotic spatial manifold. 

The Euclidean path integral is well-defined in the semi-classical approximation, where 
one looks for saddle points of the Euclidean action. We can expand the fields ty 1 around 
the saddle p pintail 

¥ J = ¥g + $¥ J (6.29) 

and write the perturbed action as 

= I [¥] + 6 2 I[¥] . (6.30) 

The first order perturbation vanishes by the equations of motion and the second order 
perturbation can be written in terms of the Euclidean Lichnerowicz operator Al as 

(6.31) 

= fd*x y /to6* I (x)(A L ) I j6* J (x') . 
35 In what follows, a subindex will denote the saddle point expressions. 
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The Euclidean path integral Z receives two contributions: one from the saddle point 
configuration and another from the determinant of the Lichnerowicz operator, so 

Zoc (detA L )- 1 / 2 e- / °[*o] . (6.32) 

When the Lichnerowicz operator has a negative eigenvalue, in other words when the system 
of differential equations 

(A L ) U 5* J = \Qu8V J (6.33) 

has a solution with A < 0, the Euclidean path integral Z receives an imaginary contribu- 
tion, which can be seen as a signal of thermodynamic instability. Qjj is a natural metric 
on the space of perturbations defined implicitly by (|6.33p . 

The existence of a negative eigenvalue for the dimensionally reduced Euclidean Lich- 
nerowicz operator provides a natural bridge between classical and local thermodynamic 
instability. The connection to a classical instability of the Lorentzian theory arises in the 
following way. As we found in previous sections, a clear signal of a dynamical instability 
is the existence of a threshold unstable mode. Let us denote this threshold mode as Sip 1 . 
We can separate the z l dependence of this mode by setting 

Sip 1 = Re(V^ 2 W) . (6.34) 

Plugging this ansatz into the second order perturbation of the equations of motion of the 
theory leads to a system of differential equations of the form 

(A L )u6^ J = -k?gu5* J , (6.35) 

where (A.l)ij is the dimensionally reduced Lorentzian Lichnerowicz operator (see (|2.2ip ). 
After a trivial Wick rotation equation (|6.35p becomes (|6.33[) with A = —k 2 . Hence, we see 
that there is a close relation between the negative modes appearing in the dimensionally 
reduced Euclidean path integral and the momentum of the threshold unstable modes of 
the classical stability analysis. 

On the other hand, one can argue for a natural relation between the semi-classical 
path integral and local thermodynamic stability. To make this relation more precise Reall 
suggested an argument, based on an earlier construction given in [15011151] . The goal of the 
argument is to relate the Hessian g^ig^j of the action I with the Hessian Hm of the local 
thermodynamic stability analysis. For that purpose, consider a system with temperature 
T, entropy S and n additional charges Q^. The n corresponding chemical potentials will 
be denoted as V{. In terms of these variables the first law of thermodynamics takes the 
form 

dM = TdS + "idQi ■ (6.36) 

i 

In the grand canonical ensemble, the idea is to construct an n + 1-parameter family of 
configurations with Euclidean action 

n 

ui, ....,u n ) = f3M(x , ...,x n ) - S(x , ...,%) - py)ujQj(x , ...,x n ) . (6.37) 
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The generic point in this family is an off-shell background, which does not satisfy the 
Einstein equations of motion, but satisfies the appropriate Hamiltonian constraints. For a 
particular value of the n+1 parameters (ceo, x n ) = (a?o(A ^1, ^n), x n ((3, v\, u n )) 
the background becomes an exact solution of the equations of motion. At this point the 
functions M, S and Qi are precisely the energy, entropy and charges of the corresponding 
gravitational system. For other values of the parameters (xq, ...,x n ) the interpretation of 
the functions M, S and Qi is not important. 

Given the existence of such a family of configurations, we can immediately draw a 
number of interesting conclusions. First of all, since the configurations become an exact 
solution of the equations of motion at the special point 

P = (x ,....,x n ) = (x (f3,vi,...,v n ),...,x n (f3,vi,....,v n )) 



the action / is extremized there yielding 

/ dl 



\® X 3/ P,V!,...,v n 



0, j = 0,..,n (6.38) 



or equivalently the standard thermodynamic relations 



1 fdM 



P \dS 



(x ,...,x n )=P 



dM 



(x ,...,x n )=P 



l,...,n. (6.39) 



The quadratic perturbation of the action I around the special point P involves the 
Hessian matrix J( 2 ) = dx-dx k ' wm ch can be recast in terms of the Hessian matrix Hm as 

(l ( 2)) ab = M a a (H M 1 )^Ml. (6.40) 

In this relation, which can be deduced easily with the use of the first order equations (j6.38[) . 
M a a is the Jacobian matrix of the transformation from the variables (xq, x n ) to the 
thermodynamic variables fx, v n ). Therefore, assuming M to be an invertible matrix, 
we can use equation ([6.401) to achieve the desired relation between the Hessian matrix /( 2 ) 
of classical stability and the Hessian matrix Hm of local thermodynamic stability. 

The construction of the n+1 parameter family of configurations is a non-trivial exercise 
that has not been performed for the most general system. So far the construction has 
been performed explicitly for magnetically charged black brane solutions in [64J and for 
electrically charged smeared black branes in [66J. 

To complete the proof one has to verify an additional set of properties. For instance, 
one should verify that the norm of the on-shell perturbations on the space of theories 
is positive definite. This ensures that the analysis is restricted to normalizable on-shell 
perturbations excluding for example any non-physical negative-norm conformal perturba- 
tions. The positivity of the norm 



||^|| 2 = J d d x S&GijSV- 
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also ensures that the action I decreases, and therefore an instability exists, precisely when 
the Hessian t^tx^j f aus to be positive definite (see equations (|6.3ip . (|6.33p ). 

Finally, one should demonstrate that there is sufficient overlap between the off-shell 
deformations (xq, ...,x n ) with the actual on-shell perturbations St/) . In [M] it was pointed 
out that a path in the family of off-shell configurations is not related directly to an eigen- 
function of the Lichnerowicz operator, but rather some linear combination of eigenfunc- 
tions. This suggests that when the action decreases along this path, at least one of the 
eigenvalues of the Lichnerowicz operator must be negative and therefore an actual on-shell 
instability should exist. It is not immediately obvious, however, that the converse is also 
true. To complete the proof one should demonstrate that the off-shell deformations and 
the actual on-shell perturbations cover the same linear space. This is a rather subtle point 
that has not been demonstrated fully in any known example. It is natural to expect in 
general that there will be sufficient overlap between the eigenfunctions of the Lichnerowicz 
operator and the family of off-shell deformations for systems that are specified uniquely 
by a complete set of conserved charges. Systems with hair do not fall into this category 
and are therefore the prime candidates as systems that violate the CSC in its current for- 
mulation. In such situations the full set of conserved charges is not enough to determine 
uniquely the classical solution. As we saw the conserved quantities are the only input of 
the local thermodynamic stability; a dynamical instability, however, is a property of the 
entire solution. 

This crucial point was put forward in [BTJ , which presented a set of theories where the 
exact black brane solutions are specified by a set of conserved charges and the asymptotic 
values of a set of scalar fields. It was found that these systems can develop a classical 
instability associated with the dynamics of the scalar fields despite the presence of local 
thermodynamic stability. 

The violation of the CSC observed in [67] is reminiscent of an earlier counterexample 
proposal in [152]. The central objects in that paper are three-charge D1-D5-P spinning 
black brane solutions in type IIB supergravity compactified on T 4 . This class of solu- 
tions contains gyrating black strings [153] carrying the same charges as ordinary spinning 
black strings, but a different value of an extra parameter, which measures the angular mo- 
mentum carried by gyrations of the string (as opposed to spin) . [152] argued that a small 
perturbation of a non-gyrating black string could grow to become large for sufficiently large 
angular momentum. This new instability was proposed as a possible counterexample to 
the CSC. 

Given the explicit violation of the CSC in the above examples, is it still possible to find a 
suitable refinement or restriction of the conjecture that can survive as a correct statement? 
The above counterexamples suggest two natural possibilities. If we choose to leave the 
definition of local thermodynamic stability unchanged in the formulation of the conjecture, 
then it is natural to add a further assumption and conjecture that the CSC works provided 
that there is a unique solution with a spatially uniform horizon and specified conserved 
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charges [67 1. No counterexamples of this restricted version of the CSC are known. An 
alternative, more ambitious, refinement of the conjecture has been suggested by B. Kol 
(see [67J). We can redefine the meaning of local thermodynamic stability by enlarging the 
Hessian matrix to include derivatives with respect to those quantities that characterize the 
asymptotics of the scalar fields or any other quantity that determines uniquely a solution. 
In that case, local thermodynamic stability translates into a positivity criterion for the 
new enlarged Hessian matrix. It seems quite likely that with this refinement the CSC 
survives all tests known up to date. 

7 Holographic implications 

Holography provides a non-trivial relation between a gravitational theory in the bulk and 
a non-gravitational theory (e.g. a Yang-Mills theory) on the boundary. A well-known 
example of this relation is the AdS/CFT correspondence between the M = 4 Super- 
Yang- Mills theory and type IIB string theory on AdS 5 x S 5 [EJ [20l EU |80] . In general, 
holography relates gravity in the near-horizon background of a brane configuration in 
string/M-theory and the non-gravitational world-volume theory living on the branes in a 
low-energy decoupling limit |22jP^l 




In this section we discuss gauge/gravity duals in which the gravity side exhibits a 
Gregory-Laflamme instability. The presence of this classical gravitational instability in 
the bulk has important consequences for the gauge theory on the boundary, which we 
review following [561 E3 1155} 166"! [58] . Focusing on a set of examples that are closely 
connected to the material of this review we will consider first D-branes on a circle, where 
the relevant gauge theory is super- Yang-Mills theory on a circle. This example can be 
generalized with D-branes on higher-dimensional torii, which we discuss briefly at the end 
of the first subsection. Other interesting examples include the M2- and M5-brane on a 
circle, where the relevant non-gravitional theory is (2+l)-dimensional super- Yang-Mills 
and Little String theory. We review the basic features of both cases giving more emphasis 
to the second example, which exhibits new properties not present in the other cases. We 
conclude with a brief discussion of the implications of Gregory-Laflamme instabilities for 
non-commutative Yang-Mills and non-commutative open string theories, which arise as 
part of the worldvolume dynamics of D-brane bound states. 

7.1 D-branes on a circle 

Consider the supergravity solution of a large number of coincident D(p — l)-branes on 
a transverse circle. In Section 15.11 we found that there is a map between the phases of 
neutral Kaluza-Klein black holes (see Section [2D and the phases of non- and near-extremal 
D-branes on a transverse circle. Furthermore, we saw in Section [5.21 that we can map the 

36 See also the review [154] and references therein. 




76 



Gregory-Laflamme mode of the neutral black string to a Gregory-Laflamme mode of non- 
and near-extremal D-branes smeared on a transverse circle. In this section we will discuss 
the consequences of this unstable mode for the holographically dual theories living on the 
D-branes. 

The low-energy theory on N coincident Dp-branes in type II string theory in flat 
space is given by the {p + l)-dimensional U(N) super- Yang-Mills theory with sixteen 
supercharges (on IR^' 1 ). According to holography |19|, I22j. the large N, finite temperature, 
SU(N) version of this theorvP^l is equivalent to the near-extremal limit of the supergravity 
solution of N coincident Dp-branes. The perturbative Yang-Mills description is valid at 
small 't Hooft coupling and the supergravity description at large 't Hooft coupling. 

Similarly, since a near-extremal D(p — l)-brane on a transverse circle is T-dual to a 
Dp-brane wrapped on a circle we conclude that the gauge theory dual of N coincident near- 
extremal D(p — l)-branes on a transverse circle is (p+ l)-dimensional SYM on M^ -1 ' 1 x S 1 
with gauge group SU(N) and sixteen supercharges |156|, I157| 11581, I160| H6lj . 

For example, this implies a correspondence between the phases of near-extremal D2- 
branes on a transverse circle and the phases of N = 4 SYM theory on M 2 ' 1 x S 1 . By 
Wick rotating to Euclidean space, where the non-zero temperature translates to a compact 
direction, we get M = 4 SYM on R 2 x T 2 , the two-torus being T 2 = Si x S 1 . Two other 
interesting examples are near-extremal DO- and Dl-branes on a transverse circle. These 
are respectively dual to (1 + l)-dimensional SYM theory on T 2 = Si x 5* 1 and to (2 + 1)- 
dimensional SYM on R x T 2 (in Euclidean space). In what follows we will restrict the 
discussion to the cases of DO-, D1-, D2- and D3-branes on a circle, but one should keep in 
mind that the results of this section can be easily extended to include the F-strings and 
the remaining D-branes. 

The main point of the discussion here is that we can use the gravity results of Sections 
15.11 and 15.21 to obtain predictions for the thermodynamics of the dual strongly coupled 
gauge theories described in the previous paragraphs [561 E3 ESI EH]. We will see, for 
example, that the localized phase of the supergravity phase diagram corresponds to the 
low temperature/low energy regime of gauge theory, while the non-uniform phase corre- 
sponds to a new phase that emerges from the uniform phase. The latter describes the 
high temperature/high energy regime of gauge theory. The Gregory-Laflamme instability 
in supergravity has a natural counterpart in gauge theory as an instability of the high 
energy/temperature phase when the energy/temperature becomes sufficiently low. 

The dictionary between the supergravity parameters of the near-extremal Dp-branes 
on a (transverse) circle and the parameters of SYM theory has been established in Section 
[5j We summarize it here for the convenience of the reader. On the gravity side, L and 
L denote respectively the circumference of the transverse circle and its T-dual. g s and g s 
denote the string coupling in the theory before and after T-duality. These quantities are 
related to each other via the relations (|5.47p . On the gauge theory side, ./V is the rank 

37 The (7(1) part of the gauge group decouples in holography [154] . 



77 



of the SU(N) gauge group and A = g\ M N = (2ir) p ~ 2 g s % 3 the 't Hooft coupling. The 
energy and temperature of the field theory will be denoted by E and T respectively. The 
definition of the near-extremal limit (|5,19p on the gravity side involves two parameters g 
and I, which are kept fixed. These parameters can be re-expressed in terms of the SYM 
variables in the following wajf] 

9 = tA^^- f AL 3 ~ p ) 2 , I = ^=a/aI^ . (7.1) 



(27t)3at2 Vp ^ v ; ' ^ 

From g, I and T we can form two independent dimensionless parameters, l/g and IT, which 
will play an important role below. These parameters are functions of the dimensionless 
temperature and coupling if = LT, A' = AL 3 ~ P . 

It should be noted that stringy corrections to the supergravity description of coincident 
Dp-branes wrapping a circli 3 ^! 



cw = R2 ' 



Rq p \ ,.o dz 



dt z + 



R 7 -p / (2tt) 



+ ' " dfl" + sJh d R-dQl_ p 



are negligible when the string scale is small compared to the length scale generated at 

p-3 1 I — I 

the horizon ijj ~ \fa'h l / 4: Rq 4 L~ 2 , i.e. when £jj v / o'Cj This is equivalent (up to 
a constant factor with powers of ir that we drop) to the condition LT -C (AL 3- ^)^^ 

R {?-p)/2 L 1/2 

On the other hand, the typical mass of the winding modes is M w ~ ,, — , so the 



effects of the winding modes can be neglected when £hM w S> 1, i.e when Rq 1. This 
is equivalent to LT » (AL 3 ~ P )~5. Hence, when p = 1,2 and A' = AL 3 " P » 1 the 

supergravity solution 17.21 is valid for a large range of temperatures t' . For p = 3 the 

— - — - —l 

supergravity solution is valid when A' 2 <C t -C 1 and for p = 4 when A' 2 <C t' <C A' 

The latter condition implies A' <C 1. 

In a similar fashion, the T-dual background of smeared D(p — 1) branes (see Eqs. 

(I5.26p - (|5,28p ) is valid (i.e. stringy and winding mode corrections can be ignored) when 

t' «; A'3^, t' <C 1. 

Having said this, we are now ready to use the supergravity results of Section [5] to make 
definite predictions for the phase diagram of the dual gauge theory at strong coupling. 



38 We use the notation and conventions of Ref. [145] . 

39 The metric is presented here in the string frame. Compare with Eqs. (|5,26|1 - (|5,28H in Section [5. II where 
the T-dual solution is written in the Einstein frame. The full background also includes a non-trivial dilaton 
and a p-form RR potential which can be deduced easily from the corresponding expressions in Section [5] 

40 The extra factor L~ x l 2 in the definition of £h follows from the extra rescaling that we do when we 
take the near-extremal limit (see sentence below (15.22[) ). 

41 To derive this result we use the definition of hd in (|5.28[) together with (|7.1|) and the relation 

iy _ 2(2^;f 
(6-p)V(6-p)n 8 - P ' 
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For that purpose we recast the thermodynamics of Section in terms of gauge theory 
variables. In (15. 24ft we defined and computed the rescaled energy, temperature, entropy 
and free energy quantities e, t, s and f = e — ts. In terms of these variables we obtain the 
dual gauge theory entropy and free energy S and F by reinstating the dimensions 

S(E) = -s(gE) , F(f) = -f(lf) . (7.3) 
9 9 

In the ensuing we concentrate mostly on the canonical ensemble and summarize the phases 
of (p+ l)-dimensional SYM theory on HP" 1 ' 1 x S 1 that follow from the uniform, localized 
and non-uniform phases of near-extremal D(p — l)-branes in supergravity. 

7.1.1 Uniform phase: High temperature phase of SYM on a circle 

The uniform phase of near-extremal D(p — l)-branes on a circle corresponds to the high 
temperature phase of (p + l)-dimensional SYM theory on R p— 1,1 x S . This is physically 
sensible because high temperatures correspond to short distances. In this regime the 
details of the compact direction in RP -1,1 x S 1 are invisible and the circumference L 
appears only as a trivial proportionality factor in the free energy. This is precisely what 
we find in the uniform phase that corresponds to a uniformly smeared near-extremal 
D(p — l)-brane (recall that this solution has the same thermodynamics as a near-extremal 
Dp-brane wrapping the T-dual circle). We can check this statement explicitly using the 
thermodynamics (|5.30[) and the dictionary (|7.3p . The resulting free energy is 

^SYM( P+ i)Cn = % P -i)(T) = -k^LN'X-^f 2 ^ , (7.4) 

where the coefficients k p are given in Table [5j The expression (|7.4p is indeed the familiar 
result for the free energy of the Dp-brane theory, with V p = V p ^\L. 

7.1.2 Localized phase: Low temperature phase of SYM on a circle 

The localized phase corresponds to the low temperature phase of (p+ l)-dimensional SYM 
theory on MP" 1 - 1 x S . Low temperatures correspond to large distances in the SYM theory, 
so in this regime the presence of the circle enters in a non-trivial manner and the whole 
tower of Kaluza-Klein states can contribute. We will see that as we send the temperature 
to zero, the leading order thermodynamics is captured by the thermodynamics of the near- 
extremal D(p— l)-brane, which is dual to the strongly coupled p-dimensional SYM theory 
on RP -1 ' 1 . The next-to-leading order corrections to the thermodynamics of the localized 
phase (see Section I5.ip will lead to gauge theory corrections that capture the contribution 
of the Kaluza-Klein modes on the circle. 

Using the expressions for the corrected free energy in (|5.37p together with the dic- 
tionary (|7.1|) . (|7.3|) we deduce that the free energy of the localized phase of thermal 
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(2 21 3 2 5 7 7" 19 ^ 14 ) 1 /5 
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(2 13 3 5 5- 13 ^ 8 ) 1 /3 
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2 5 3- 7 vr 2 



Table 5: Coefficients for the free energy of Dp-branes. 



(p + l)-dimensional SYM theory on MP 1,1 x S 1 is 



TiloC 

SYM(p+l) 



(7.5) 



4— p 
6— p 



,8-p 



T 6=P ^ 1 + 



2(9-p)C(8-p)fc p - 
(6-p) 2 (2vr) 3 ^ 



9-p 



6— p 



The coefficients fc p are the same as above and are listed in Table [5T 12 I 

Notice that the leading term in (|7.5p is in perfect agreement with the expected result 
for the strong coupling limit of the p-dimensional SYM theory. In particular, the 't Hooft 
coupling appears through the combination A/L, which is precisely what we expect to find 
by compactifying the (p+l)-dimensional theory with 't Hooft coupling A. The second term 
in (|7.5p gives a quantitative prediction for the first correction in terms of the dimensionless 
parameter 

5 = LfV AL 3 -p . (7.6) 

o 8-p 

The next corrections will be of order [8 6 -p) . It would be very interesting to reproduce 
these corrections from the field theory side. We refer the reader to [57] for the explicit 
expressions of (|7.5p in the cases p = 1, 2, 3 and 4. 

Similar results in the microcanonical ensemble can be easily obtained using the cor- 
rected entropy (|5.35p of the localized near-extremal branes and (|7.ip . (|7.3p . 



7.1.3 Non-uniform phase: New phase in SYM on a circle 

The non-uniform phase of near-extremal D(p — l)-branes on a circle (see Section 15. ip 
predicts the existence of a new phase of (p + l)-dimensional SYM theory on RP -1 ' 1 x S 1 
at intermediate temperatures. With the use of (|5.33p and the dictionary (|7.ip . (|7.3p we 
can also deduce the first correction to the thermodynamics around the point where the 
non-uniform phase connects to the uniform phase. One finds the free energy 
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Z7 nu 

*D(p-l) 



(T) 

(27r) 3 Vp^iiV 2 AL 3 - 2p 



5 — p 



(7.7) 



9-p 

The critical temperature in this expression is given by 



1 + 



t c c / T 



2vr 



L V AL 3 "P 



(7.8) 



These coefficients are related to k[ d) defined in {ppl via the relation iff p, (27r) 2; s-; = fc 



r>(9-p) , 



S(J 



and the values for t c , e c , s c and c can be read off from Table 0] by using d = 10 — p. 
Notice that the first term in the expression (17.7P corresponds to Eg YM( ^ p+1 ^(T c ), i.e. the 
free energy (|7.4p of the uniform phase evaluated at the critical temperature. The second 
term is — $sym(p+i)( T C )(T — T c ) and therefore involves the entropy of the uniform phase 
evaluated at the critical temperature. The third term contains the departure of the free 
energy from the uniform phase into the non-uniform phase. It would be very interesting 
to understand the results (|7.5p and (|7.7p for thermal {p + l)-dimensional SYM theory on 
R x x 5 1 from the gauge theory side. 

Again, similar results for the non-uniform phase in the microcanonical ensemble can 
be obtained using the corrected entropy (|5,32p of the localized near-extremal branes and 

(HD, 

Finally, since both the localized and non-uniform phase of near-extremal branes have 
copies, we also have copies of the thermal (p+ l)-dimensional SYM phases obtained above. 
See Ref. [57] for further details. 

7.1.4 The expected T(E) diagram for D-branes on a circle 

As we discussed above, the thermodynamics of near-extremal Dp-branes on a transverse 
circle, is related to the thermodynamics of phases of Kaluza-Klein black holes in 10 — p 
dimensions. In particular, for the cases of interest in this section, namely DO, Dl, D2 and 
D3-branes, we need the phases of Kaluza-Klein black holes in 10, 9, 8 and 7 dimensions. 
However, at present we do not have complete analytical or numerical data (see, however 
|33j) for the localized and non-uniform phases in those dimensions. For this reason, our 
results so far (c./. the free energies (|7.5p and (|7,7p for the localized and non- uniform phase 
respectively) have been leading order results around some specific point in the phase 
diagram. For the localized phase this point is the origin (large radius limit) and for the 
non-uniform phase it is the GL point. For the uniform phase we have of course the exact 
result flU). 

Nevertheless, based on the expectation that the (yu, n) phase diagram for neutral 
Kaluza-Klein black holes in those dimensions has a similar form as those in 5 and 6 
dimensions (see Fig. [3J we can use the map (|5.23p to make some general observations 
|58j . For example, the map implies that the (E,r) phase diagram for the near-extremal 
Dp-branes has the same qualitative form as that in Fig. From this we get a prediction 
for the qualitative form of the T{E) diagram, because given the curves in the (E, r) phase 
diagram the entire thermodynamics is determined. 

For instance, consider the case of N near-extremal D2-branes on a transverse circle, 
which is dual to finite temperature N = 4 supersymmetric SU(N) Yang-Mills on M 2 x T 2 . 
A similar qualitative picture is expected also for the other Dp-brane cases. In the uniform 
phase, which consists of near-extremal D2-branes uniformly smeared on a transverse circle, 
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the relation between temperature and energy is 

T{E) oc , (7.9) 

We argued that the uniform phase dominates at high energies. In the localized branch 
which dominates at low energies we know [57] that the temperature/energy relation is 

T(E) oc ET5 (1 + aE) , E < 1 , (7.10) 

with a < a known constant. Finally, there is a non-uniform phase which emanates from 
the uniform phase at the critical energy E c = 0.71 • (2-7r) 3 ^iV 2 /(A-L 3 ). At the critical 
point, the ratio of the specific heat of the non-uniform phase to the specific heat of the 
uniform phase is [57] 

— = \ = 0.41. (7.11) 

c u l + f7£ c 

The ratio is less than one, a feature that holds for any of the Dp-brane cases mentioned 
above. 

The salient features of the emerging picture in the T(E) diagram are the following 

• Both the uniform and localized phases begin at the origin (with infinite slope) and 
at very low energies the uniform phase has higher temperature than the localized 
phase (at a given energy). As the energy is increased the localized phase is expected 
to cross the uniform phase at relatively low energies below the critical energy where 
the non- uniform phase emerges. 

• The slope in the T{E) diagram is the inverse of the specific heat. Therefore we 
deduce from (|7,lip that the slope of the non-uniform phase at the critical point 
(E, T) = (E c , T c ) is larger than that of the uniform phase. 

• It is natural to assume that the localized and non-uniform branches meet at a merger 
point in a qualitatively similar way to the one observed explicitly in the NS5-brane 
case in [58] (see Fig. [T2|) . 

The emerging picture is summarized in Fig. [TUl 

We notice that as we increase the energy above extremality in the localized phase we 
first reach a point where the specific heat becomes infinite. This is the point where the 
localized phase attains its maximal temperature. Increasing the energy further we reach a 
point where the specific heat vanishes. This is the point where the localized phase attains 
its maximal energvF 3 ] Despite the presence of a maximal temperature in the localized 
phase, the uniform phase exists for any temperature and therefore no limiting temper- 
ature exists for the system of near-extremal branes on a circle. Hence, for sufficiently 



43 The conditions for infinite and zero specific heat are in this case Er'(E) = (5 — p)/2 — r(E) and 
l/r'(E) = respectively. 
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Figure 10: The qualitative picture of energy versus temperature for the phases of near- 
extremal Dp-branes on a transverse circle. One can see the uniform (red), non- uniform 
(blue) and localized (magenta) phases. The dashed curve connects the localized and non- 
uniform phase. 

high temperatures the system will necessarily find itself in the uniform phase. At low 
temperatures the localized phase will be favored. 

The above behavior is the same in both the canonical and microcanonical ensembles. 
In Section 17.2.21 we will see that the case of NS5-branes is very different in this respect. In 
that case the uniform phase has a constant Hagedorn temperature and there is a striking 
difference between the canonical and microcanonical ensembles. 

7.1.5 From strong to weak gauge coupling 

So far we have obtained a set of predictions for strongly coupled (p + l)-dimensional 

SYM theory on a circle by analyzing the supergravity phases of near-extremal black hole 

solutions on a transverse circle. It is interesting to ask what happens to all these phases 

in gauge theory when we continue the theory from strong to weak 't Hooft coupling A. 

If there is a line of phase transitions at strong coupling, this line is not expected to end 

abruptly in the (A,T) parameter space as we change the coupling A. It may either go to 

zero or infinite temperature at a finite value of the coupling, or it may continue smoothly 

to weak coupling, where one can analyze it with standard perturbative methods directly in 

gauge theory. An example of the latter possibility is the Hawking-Page transition [50] 

in AdSs, which is dual to a thermal deconfinement transition in (3 + l)-dimensional M = 4 

SYM theory on 5 3 . This transition can indeed be continued to weak coupling where it can 

be analyzed explicitly |81[ [82] [83] . A similar continuation to weak coupling also occurs 

for the case of interest in this section, namely the Gregory-Laflamme instabilitjQ and the 

associated transition in the dual SYM theories on a circle. This has been shown for the 

44 In another recent example, Ref. [162] found the weak coupling version of the Gregory-Laflamme local- 
ization instability of the small AdS^ black hole. 
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(1 + l)-dimensional maximally supersymmetric YM theory in [SB]. We review the main 
features of this example in the ensuing. 

The Euclidean (1 + l)-dimensional SU(N) SYM theory on the torus S\ x S 1 has two 
non-contractible Wilson loops, the loops running around the time r and circle z. The 
spatial Wilson loop 

Vz = Pe i§dzA z ^ 12 ) 

is of special interest here, because it serves as an order parameter for the weak coupling 
manifestation of the Gregory-Laflamme transition. Indeed, the eigenvalues of V z represent 
the positions of T-dual DO-branes on a circle and their distribution can distinguish between 
the three different phases. The uniform phase corresponds to a uniform eigenvalue distri- 
bution, the non-uniform phase corresponds to a non-uniform eigenvalue distribution that 
breaks the translational invariance on the transverse circle and the localized phase maps to 
an eigenvalue distribution, which is localized on the transverse circle. [56j showed that the 
eigenvalue distribution of the Wilson loop V z is sharply localized at small temperatures 
and weak coupling. As we continue the theory to higher temperatures the eigenvalues 
spread out and fill out the whole circle. There is a critical temperature T c ~ 1.4/(AL 3 ) 
where the theory undergoes a Gross- Witten like phase transition [124 J from the local- 
ized distribution to the uniform distribution. For more details we refer the reader to the 
original discussion of the Refs. [56 | 1163] . 

7.1.6 Implications of the GL instability for higher tori 

Above we considered the implications for Yang- Mills theory on a circle due to the Gregory- 
Laflamme instability of D-branes smeared on a transverse circle. Similar results can be 
obtained for Yang-Mills theory on a torus T k by analyzing the Gregory-Laflamme insta- 
bility of D-branes smeared on a transverse torus T k . 

As a special case, consider the J\f = 4 SYM theory on M x T 3 . The holographic dual of 
this theory arises from DO-branes smeared on the three-torus T 3 . This brane configuration 
has Gregory-Laflamme instabilities in the periodic directions when the radii are sufficiently 
large. This implies an interesting phase structure which generalizes the phase structure 
of non-extremal D-branes on a circle presented in Section 15.11 First, for T 3 = S 1 x T 2 
we can obtain several phases for the DO-branes on T 3 from the phases of D2-branes on a 
circle with T-duality, imposing that there is no dependence on the T 2 directions. Thus, we 
inherit all the phases of the D2-brane on a circle. Moreover, it is conceivable that there are 
new localized phases, for example, phases in which the DO-brane horizon is topologically 
S 3 x S 5 corresponding to having the DO-branes localized in all three torus directions, or one 
could have that the DO-brane is localized in two directions giving an 5 1 x 5 2 x S 5 horizon. 
In addition, one can well imagine new non-uniform phases, for example in which the DO- 
branes are not localized in any of the three directions but nevertheless are non-uniform in 
all three directions. This would correspond to a T 3 x S 5 topology of the horizon. Therefore, 
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we see that the phase structure of DO-branes on T 3 is even richer than that of D-branes 
on a circle. This is in accordance with the phase structure of neutral p-branes on tori as 
discussed in Section [4.31 In particular, via the boost/U-duality map discussed in Section 
15.11 one can directly relate the phases of neutral Kaluza-Klein black holes on .M 1,6 x T 3 
to those of non- and near-extremal DO-branes smeared on T 3 . Through holography, the 
richness of the phase structure for DO-branes on T 3 implies an equally rich phase structure 
for M = 4 SYM theory onlxT 3 , We refer to [163] for more considerations on the phase 
structure of gauge theories on tori. 

7.2 M2- and M5-branes on a circle 

In this section we consider separately the cases of near-extremal M2- and M5-branes on 
a transverse circle. By IIA/M S-duality, the former is dual to the near-extremal D2- 
brane and the latter is dual to the near-extremal NS5-brane. Hence, it is conjectured 
that N 3> 1 coincident near-extremal M2-branes on a circle are dual to (2+l)-dimensional 
SYM theory with sixteen supercharges and gauge group SU(N) [T5JI22]. Similarly, the 
system of A ;§> 1 concident near-extremal M5-branes on a transverse circle is dual to the 
six-dimensional (2, 0) LST. In both cases, the transverse circle is the M-theory circle and 
the dual non-gravitational theory lives on the non-compact space W' 1 . p = 2 for the case 
of the M2-brane and p = 5 for the case of the M5-brane. 

7.2.1 (2 + l)-dimensional SYM theory 

In this case we place N near-extremal M2-branes transverse to the M-theory circle, which 
we take to have circumference L = 2irg s l s . Using the quantization condition that comes 
from the M2-branes together with (|5.19p and (|5.22p we find that we can recast the super- 
gravity parameters g, I as 

(2vr) 2 A 3 A 3 / 2 , . 

9 = -vt^> / = 27r ~ ' (7 - 13) 

where A = <?ym 2 A = g s l~ 1 N. As in the cases of the previous subsection, these relations 
provide a dictionary between the supergravity results of the near-extremal M2-branes on 
a circle and the dual (2 + l)-dimensional SYM theory. Since there are many similarities 
between this case and the case of Dp-branes on a transverse circle of the previous subsec- 
tion, we will provide here only a brief summary of the main results and refer the reader to 
[57] for a more detailed discussion. Qualitatively similar results can be obtained in both 
the microcanonical and canonical ensembles. 

The uniform phase of near-extremal M2-branes on a circle corresponds to the high 
temperature/high energy phase of thermal (2 + l)-dimensional SYM. The localized phase 
corresponds to the low temperature/low energy phase. At leading order this phase is 
captured by the near-extremal M2-brane solution and is dual to the infrared fixed point of 
(2 + l)-dimensional SYM which is a superconformal field theory with SO (8) R-symmetry. 
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Using the results of Section and the dictionary of eqs. (|7.3|) and (|7.13|) we can compute 
the free energy for both of these phases. In the uniform phase the result reproduces the 
free energy of the near-extremal D2-brane theory. In the localized phase the leading order 
result reproduces the free energy of N near-extremal M2-branes. The small black hole 
corrections capture the non-trivial effects of the theory as we move away from the infrared 
fixed point. 

The phase diagram of the near-extremal M2-brane on a circle also predicts the existence 
of a new non-uniform phase of uncompactified (2 + l)-dimensional SYM theory on R 2,1 
at intermediate temperatures. This phase emerges from the uniform phase at a critical 
temperature and presumably connects to the localized phase at a merger point. The 
first correction to the thermodynamics around the critical point of the uniform phase can 
be computed from the available supergravity data as above. It would be interesting to 
understand these results from the gauge theory side. Finally, since both the localized and 
the non- uniform phase of near-extremal branes have copies, there are similar copies of the 
thermal uncompactified (2 + l)-dimensional SYM phases discussed above. 

7.2.2 Little String Theory 

In the decoupling limit of N coincident NS5-branes in type IIA string theory the string 
length l s is kept fixed and the string coupling g s is sent to zero [1641 1165] . In this limit 
the dynamics of the theory is believed to reduce to a string theory without gravity called 
Little String Theory (LST), or more precisely (2,0) LST of type A N _i [Ml Q2EZ1 QES 
11691 1 170 J O In this subsection (see also [1551 EH] ) we wm consider the thermodynamics 
of the near-extremal NS5-brane in type IIA string theory using as a starting point the 
S-dual description as a near-extremal M5-brane on a transverse circle and the input of 
the phases obtained from six-dimensional neutral Kaluza-Klein black holes in Section [3j 
Other interesting work on the thermodynamics of LST can be found in [1711 11641 11721 

The thermodynamics of the NS5-brane exhibit a variety of new features compared to 
the other cases analyzed earlier in this section. The usual near-extremal NS5-background 
(dual to M5-branes uniformly smeared on a transverse circle) has entropy as function of 
energy S(E) = T^E, where Th g is a fixed Hagedorn temperature. In [1721 1173] the idea 
was put forward that the singular thermodynamics of the NS5-brane background could 
be resolved by computing a string one-loop correction to the background. This was done 
in [174] . where it was shown that the corrected background has T > Thg and negative 
specific heat. In a subsequent analysis the results of [174] were interpreted to indicate 
that one does not have a Hagedorn phase transition and that the Hagedorn temperature 
is a limiting temperature [179] . 

45 A similar decoupling limit of NS5-branes in type IIB string theory gives the (1, 1) LST of type An-i- 
In this review we will focus on the type IIA NS5-branes, which are related more directly to the transverse 
M5-branes of M-theory. 
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Here we will see, following [58J, that in the canonical ensemble the usual near-extremal 
M5-brane background is thermodynamically subdominant at tree-level to a background 
of near-extremal M5-branes localized on the transverse circle. The new localized phase 
exhibits a maximal temperature and it remains unclear what happens if one tries to 
heat up the type IIA NS5-brane to a temperature higher than this maximal temperature. 
A related additional new feature of the NS5-brane is the fact that the canonical and 
microcanonical ensembles seem to be inequivalent at large energies. In what follows we 
will review separately the basic features of each ensemble. 

Canonical ensemble 

We begin with the canonical ensemble, since this is the ensemble that exhibits the most 
interesting results. The free energy versus temperature is displayed in Fig. [TT1 




Figure 11: Free energy versus temperature for the near-extremal NS5-branes. Displayed 
are the non-uniform (blue) and localized (magenta) phases. The uniform phase is located 
at the point (T,F) = (T^ g ,0). 



Before considering the new data for the localized phase, we list the main properties of 
the three different phases (already stated in [S7]) that arise from the corresponding phases 
of neutral Kaluza-Klein black holes @ 

• Uniform phase. This phase has zero free energy and a fixed temperature T = T^ g 

46 Note that the non-uniform and localized phase have copies with free energy F' = F/k 2 and temperature 
T' = T (k — 2,3,...). We do not consider the copied phases here since they are subdominant in this 
ensemble. 
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and captures the well-known Hagedorn-like behavior of LST. It corresponds to a 
single point in the free energy versus temperature diagram in Figure [TTJ 

Non-uniform phase. This phase emerges out of the point (T,F) = (T^ g ,0) and 
exhibits increasing temperature and decreasing free energy. It has positive specific 
heat. For temperatures near the Hagedorn temperature < T — Th g <C Th g , one 
finds 

T \ ( T 



F m (T) 



-2irV 5 N 3 Tf 



1.54- 



1 +3.23 



i 



lis 



(7.14) 



This expression has been computed in [57] using results of [3T] . 



Localized phase. The localized phase begins at the point (T,F) = (0,0), which 
corresponds to the (2, 0) SCFT, the infrared fixed point of the type IIA near-extremal 
NS5-brane. The first correction to the thermodynamics as we raise the temperature 
and move away from the infrared fixed point is [57] 

2 6 7T 3 
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/kg, 



(7.15) 



This has been computed using the results of [39J . 

One can say more about the localized phase by using the numerical data of [46] . The 
localized phase begins at the point (T, F) = (0, 0) and has decreasing free energy until the 
maximum temperature T+ = 1.177 -T^g. At this point the curve reverses direction and has 
decreasing temperature and increasing free energy. 

The most intriguing part of the localized phase is evidently the behavior near the 
maximal temperature T*. One way to understand the behavior around this point is to 
consider the temperature versus energy diagram in Fig. [12j In this diagram the localized 
phase is apparently a smooth curve. Hence, the entropy is continuous at T ~ T*. However, 
since the temperature is maximal at T = T* the (E, T) curve has a horizontal tangent at 
that point. This implies that the localized phase has positive specific heat as T increases 
from zero to T± which becomes infinite at the maximal temperature T*. Continuing along 
the same curve in the (E,T) diagram we begin with minus infinite specific heat and 
continue with negative specific heat until the curve develops a vertical tangent. At that 
point the specific heat becomes zero. For the remaining part of the curve the specific heat 
is positive! 47 ! 

Let us consider more closely the thermodynamics around the critical temperature T*. 
From the available data we get 

^ 2 ^ = 3.60 • {2tt)-H- s 

V 5 



T~T*-c(E- E* 



0.038 



(7.16) 



47 In the [E, r) phase diagram the point with infinite specific heat is characterized by the value of r where 
Er'(E) = 1/2 — r(E). Zero specific heat occurs when r'(E) is infinite. It is not difficult to see from the 
phase diagram in Fig. 1 in [58] that for each of these two equations there exists a particular solution of r 
on the localized branch. 
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Figure 12: Temperature versus energy for the near-extremal NS5-branes. Displayed are 
the uniform (red), non-uniform (blue) and localized (magenta) phases. 



From this expression we can obtain the free energy as a function of temperature near the 
singularity 

F(T) ~F*-S*(T- r*) ± —^=—(T - T,) 3 / 2 , 
£ = 3.5 ■ (2^)- 4 ViV , -i = - SM = -0.53 • (2vr)~ 5 C 6 • 

In the first equation the ± sign corresponds to the two branches on either side of the 
maximum temperature point. The minus branch is the one that has positive specific heat 
(for this branch E < E+) and the plus branch is the one that has negative specific heat 
(for this branch E > E+). It is immediately clear from (|7.17p that the entropy is indeed 
finite and continuous at the critical point, whereas the heat capacity diverges and changes 
sign. 

Also it is clear from Fig. [11] that there are temperatures where more than one phases are 
available. The dominant phase will be the one with the lowest free energy@ It corresponds 
to the localized phase, which lies within the energy range < E < E+. Hence, in the 
canonical ensemble we conclude that the system is (globally) thermodynamically unstable 
in the uniform phase (at zero free energy and temperature T = T^g) and will decay to 
the localized solution at the same temperature. In the dominant, localized phase the heat 
capacity will increase as we increase the temperature and eventually will become infinite 



48 Indeed, for two phases 1 and 2 at the same temperature with free energies Fj < Fi the Euclidean path 
integral Z will be dominated by phase 1, since Z ~ e _/3Fl + e _/3F2 = e _/3Fl (1 + e - {F2 ' Fl) ) 
large N limit. 
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at the maximal temparature T = T+. This implies that T+ is a limiting temperature. Note 
that this limiting temperature exceeds the temperature T^ g associated with the uniform 
phase. 

Microcanonical ensemble 

We now turn to the microcanonical ensemble. The entropy of the three different phases 
as a function of energy is displayed in Fig. [13] The main properties of these phases arJ^l 

• Uniform phase. This phase corresponds to the horizontal line in Fig. [13] and exhibits 
Hagedorn behavior. The entropy as a function of energy (see (|5.29p ) is 

Su(E) = A . (7.18) 
J-hg 

• Non-uniform phase. This phase emerges at the point (E,S) = (E C ,S C ), and for 
energies close to the critical energy E c we have |57| 



S(E) 
S U (E) 



1-0.05- 



<&Lp - Ec )l 2 . £ = 1.54-(27r)~ 5 /~ fi , < «^ « ! . 

(7.19) 

This expression has been computed using the results of [31j . The entropy function 
for the entire non-uniform phase, as displayed in Fig. [131 was computed in [57J using 
the data of 1301. 



Localized phase. The localized phase begins at the origin (E,S) = (0,0). At very 
low energies we obtain, using the results of [39], the analytical result [57] 



9 (F) ~ 27rl/ % 6 )5/6 V 5 ( (2 , 5 l°sE\ 5/6 ( C(3) (27r)5 l 6 s E\ l 6 s E 

(7.20) 

The entropy function for the entire localized phase, as displayed in Fig. [T3l has been 
computed using the data of [46J and the map of [ST] . 

In the microcanonical ensemble, the phase with the highest entropy dominates. Hence, 
from Fig. [13] we conclude that for energies below E#/V 5 = 2.51 • (2vr)- 5 /7 6 , the localized 
phase is preferred in the canonical ensemble, whereas for energies above that value the 
system is in the uniform phase. In the range E# < E < E max , with E max /V5 = 3.92 • 
(27r) _5 /j 6 , the entropy of the uniform phase is greater than that of the localized phase. 
For E > E max only the uniform phase is available. Note that E# < E* < E mSuX . The 
non- uniform phase has always lower entropy than the other two phases. 

The point of the localized phase where E = E max is special. The heat capacity is 
zero there, which corresponds to a singular behavior in the microcanonical ensemble and 



49 The non-uniform and localized phase have copies with entropy S' = S/k 2 and energy E' = E/k 2 
(k — 2, 3, ...). We will not consider the copied phases here since they are subdominant in this ensemble. 
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Figure 13: Entropy versus energy for the near-extremal NS5-branes. On the vertical axis 
we plot the ratio S(E)/S n (E) and on the horizontal axis E/E c . S U (E) denotes the entropy 
function of the uniform phase of near-extremal NS5-branes. Displayed are the uniform 
(red), non-uniform (blue) and localized (magenta) phases. 

prompts the interpretation of the energy ^/ max 3S ci limiting energy for the localized phase. 
It is also interesting to consider what happens to the system with energy E < E# when 
we throw in a finite amount of energy AE. There are three possibilities. If E + AE < E# 
or E# < E + AE < E max the system will remain in the localized phase. In the second 
case, however, the entropy of the uniform phase at the new energy E + AE will be greater 
so the system will probably go through a phase transition to end up in the uniform phase. 
Finally, if E + AE > E mSjX the only available phase is the uniform phase and the system 
presumably has to end up in that phase. 

Comparing the ensembles 

Comparing the canonical and micro canonical ensembles, we see significant differences in 
the qualitative behavior of the near-extremal NS5-brane thermodynamics. In the canonical 
ensemble, the localized phase is always the dominant phase in the interval E < E+, and has 
a positive specific heat and a limiting temperature T*. In the micro canonical ensemble, 
the localized phase is again the dominant phase for energies E < E#. However, the 
qualitative behavior is very different for energies E# < E < E max . In this energy range 
the uniform phase is dominant. Hence, we find different preferred phases in the two 
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ensembles. For energies E > E mSi ^ there is an even more dramatic difference between 
the two ensembles since the only available phase is the uniform phase in the canonical 
ensemble. Consequently, it seems that the two ensembles are completely inequivalent at 
large energies. 

Since we discuss a system without gravity one would expect that it is always possible 
to bring it in contact with a heat bath of a certain temperature. Therefore, the inequiv- 
alence of ensembles could suggest a subtlety in the definition of the canonical ensemble. 
Alternatively, it could point towards the existence of a new hitherto unknown phase of 
near-extremal NS5-branes. This could potentially resolve the issue arising from the ab- 
sence of any stable phases with temperatures greater than T*. See the conclusions of 
Ref. [58] for further discussion on these issues. 

7.3 D-brane bound states 

In the above discussion we considered the implications of the Gregory-Laflamme instability 
when there is a singly-charged brane with a compact direction. It is also interesting to 
consider brane bound state backgrounds without compact directions. In this case the 
CSC conjecture still applies, as reviewed in Section [6l thus giving a connection between 
the Gregory-Laflamme instability, the local thermodynamical stability and the stability of 
the dual non-gravitational theories living on the brane bound states. 

In [59] the implications of the Gregory-Laflamme instability for the phase-structure 
of the non-commutative Yang-Mills (NCYM) theories were described. The dynamics of 
open strings ending on a Dp-brane are described at low energies by super- Yang-Mills 
theory in p + 1 dimensions. For a D-brane in the presence of an NS-NS B-field the same 
dynamics is captured by a non-commutative gauge theory [181} \182\ 11831 1184] . One can 
then obtain the holographic dual of the NCYM theory by taking the near-horizon limit of 
./V parallel D-branes in the presence of an NS-NS B-field. This implies, in particular, that 
the near-extremal limit of the thermal D1-D3 brane bound state is holographically dual to 
a four-dimensional NCYM theory with a non-zero temperature [1851 1186| 11461 [TSTl 188J. In 
|59j it was argued that the D(p — 2)-Dp brane bound state is unstable in the near-extremal 
limit for all temperatures above zero. This can be seen using the condition (|6.23[) for 
thermodynamic stability of the D(p — 2)-Dp brane bound state derived in Section [6.21 

Ref. [61] discussed similar implications for the non-commutative open string (NCOS) 
theories living on the brane bound states. The effective low energy dynamics of open 
strings ending on an Fl-Dp brane bound state are described by a Non-Commutative 
Open String (NCOS) theory, which is a theory of open strings on a target space in which 
space and time does not commute [1891 1190] . Taking the near-extremal limit of the non- 
extremal Fl-Dp bound state we obtain a gravity dual for the {p + l)-dimensional NCOS 
theory |147j . In [61] it was argued that the near-extremal Fl-Dp bound state is unstable 
above a certain critical temperature. 

Similar considerations also hold for other non-gravitational theories on brane bound 
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states, such as the OM theory on the M2-M5 brane bound state [191, 192J and the ODp- 
theories on the Dp-NS5 brane bound state |19R I193J. 



8 Discussion and outlook 

In this review we summarized the current state of knowledge on the classical stability of 
black branes in pure gravity and supergravity. We focused mainly on the case of static 
black hole solutions that asymptote to a Kaluza-Klein space of the form M d x 5 1 . Although 
this is the simplest possible case in a large class of solutions that asymptote to Kaluza- 
Klein spaces of the form Ai d x J\f, where J\f is any compact Ricci flat manifold, we have 
seen that it exhibits a profoundly rich structure with many features that are expected 
to be part of the more general story. We would like to conclude with a summary of 
other known classical instabilities in gravity and a list of important questions that deserve 
further study. 

8.1 Other classical instabilities 

Ultra-spin instability 

Rotating black holes in higher dimensions with arbitrary angular momentum were discov- 
ered twenty years ago by Myers and Perry p3j. The absence of an upper bound on the 
angular momentum (in six and higher dimensions) should be contrasted with the situation 
in four dimensions, where the Kerr black hole obeys the famous Kerr bound J < GM 2 , 
and with the similar bound J 2 < 32GM 3 /(27-7r) for the five dimensional Myers- Perry black 
hole. Nevertheless, it was argued in pT93j that in six or higher dimensions the Myers-Perry 
black hole becomes unstable above some critical angular momentum thus recovering a 
dynamical Kerr bound. The instability was identified as a Gregory-Laflamme instability 
by showing that in a large angular momentum limit the black hole geometry becomes 
that of an unstable black membrane. This result is also an indication of the existence of 
new rotating black holes with spherical topology, where the horizon is distorted by ripples 
along the polar direction. Furthermore it suggests that a possible endpoint of the decay 
process in ultra-spinning black holes is the fragmentation of the original black hole into 
multiple black holes. Ref. [194] verified that this possibility is consistent with a global 
thermodynamic argument. An alternative conjecture for the evolution of the ultra-spin 
instability is that the decay proceeds via emission of gravitational radiation. In the pro- 
cess, the black hole sheds a fraction of both its spin and mass and settles down in a range 
of parameters within the stable regime. Comments on this alternative possibility can also 
be found in [194] . 

It should be emphasized that up to date the complete analytic theory of perturbations 
around the higher-dimensional rotating black holes has not been developed. Progress in 
this direction has been achieved in special cases (odd number of space-time dimensions, 
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equal angular momenta) in [195J. 

Although the arguments presented in [194] apply only to six or higher dimensions an 
instability in five dimensions is also expected [15] . Entropy arguments suggest that in this 
case the endpoint of the instability is the five dimensional black ring. 

Gyration instability 

In the context of D1-D5-P spinning black brane solutions in type IIB supergravity com- 
pactified on T 4 it has been argued in [152] that a small perturbation of a spinning black 
string could grow to become large for sufficiently large angular momentum. The endpoint 
of this instability is expected to be a gyrating string where part of the original angular mo- 
mentum is carried by gyrations of the string (as opposed to spin). This type of instability 
has been proposed as a possible counterexample to the CSC (see also Section [631) . 

Superradiance instability 

Superradiance is a phenomenon that occurs in a scattering process when the incident 
wave is reflected back amplified. The first classic example of superradiant scattering, 
which later led to the notion of superradiant scattering in black hole space-times, was 
given by Zel'dovich in [1961 1197] . By examining what happens when scalar waves of the 
form e ~" lujt + tm <t ) scatter off a rotating absorbing cylinder, Zel'dovich concluded that the 
scattered wave is amplified provided that the frequency u of the incident wave satisfies 
the inequality 

u < mn . (8.1) 

In this relation f2 is the angular velocity of the rotating cylinder. Similar superradiant 
scattering occurs in the background of Kerr black holes [1961 11971 H98, 199, 200J for fre- 
quencies <jj still satisfying (|8.1|) with 0, the angular velocity of the black hole. 

An instability can occur if the black hole is surrounded by a reflecting wall that scatters 
the returning wave back towards the horizon. In that case the wave will bounce back and 
forth from the horizon extracting an exponentially growing amount of energy from the 
black hole until the pressure destroys the mirror. The combined system of the black 
hole and the mirror is known as Press and Teukolsky's black hole bomb [2U11 202J. A 
natural reflecting wall occurs in a variety of situations: a massive scalar field propagating 
in a Kerr background [2031 1204, 205J, a wave propagating around spinning black strings 
[1521 12061 1207] . Kerr-anti-deSitter black holes [208J, large radius doubly spinning black 
rings [209] , In general, the unstable modes can be scalar, electromagnetic or gravitational. 

Ergoregion instability 

Space-times with ergospheres and no horizon are unstable under perturbations of scalar 
and electromagnetic fields or gravitational waves [210] . The instability also occurs in the 
case of rotating stars with an ergoregion [2111 1212] and in smooth non-super symmetric 
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D1-D5-P geometries |213| . In the latter case the endpoint of the instability has been 
argued to be a smooth supersymmetric configuration. 



A detour on the stability of black rings 

The recent discovery of the black ring solution in five dimensions [15] demonstrates clearly 
that in higher dimensions new qualitative features can appear which are not present in four 
dimensions. Most notably, one finds that there are black holes with non-spherical horizon 
topologiejf^ and that conventional notions of black hole uniqueness do not apply. This 
realization has made the five dimensional black ring a protagonist in the arena of higher 
dimensional black hole physics. Some important partial results on the classical stability of 
the five-dimensional black ring have been obtained. Here we review the emerging picture 
beginning with a quick reminder of the major features of black rings. We will see that 
black rings can exhibit a combination of the above-mentioned instabilities and more. We 
will focus mostly on neutral black rings. 

The neutral black ring metric has been discussed in a series of papers [T5 l 12161 Ill5j (f° r 
a recent review of the subject with a complete list of references see [217] ). The resulting 
one-parameter family of solutions is conveniently parameterized by a single variable v in 
the range < v < 1. In terms of v the dimensionless angular momentum and horizon 
area are given by the expressions 



J = TT7^7T7^ = — o ' a H 



,_27«J'_(l + ^ ___S./» A„ _ 2y/ ^r^_ ( 8 . 2) 



32G M 3 ' n 16 V vr (GMf/ 2 

Compare these with the corresponding expressions for the five dimensional Myers-Perry 
black hole 

a H = 2j2-— . (8.3) 



1 + v 1 V 1 + v 

The situation is summarized for fixed mass in Fig. [TJ] along with the corresponding curve 
for the Myers-Perry black hole. This diagram has two striking features. First, there is 
a range of j where three objects with the same mass and spin coexist (one Myers-Perry 
black hole and two black rings). Secondly, the black ring curve exhibits a cusp at v = | 
which is the meeting point of two branches of black rings. One branch has < v < 1/2 
and has been called the branch of thin black rings, because it includes the regime v ~ 
of very thin (highly spinning) rings. The other branch with 1/2 < v < 1 has been called 
the branch of fat black rings and terminates at the singular solution with v = 1. 

The classical stability of these solutions has been studied with a combination of meth- 
ods (qualitative and semi-quantitative) in [El \2M IH3 EUl [22H dm [2221 E21 and has led 
to the following picture (this picture has been recently summarized and further elaborated 
upon in [223J to which we refer the interested reader for more details): 



50 The horizon topology of black rings is S 2 x S 1 . In general, topological arguments allow for a large 
class of horizon topologies in higher dimensions. For recent discussions of the relevant theory see 214, 215 . 
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Figure 14: Rescaled horizon area a# versus angular momentum squared j 2 for five- 
dimensional Myers- Perry black holes and black rings. The green curve represents the 
Myers-Perry black hole, the red curve (for v > ^) the branch of fat black rings and the 
solid black curve (for < v < 5) the branch of thin black rings. 

• At large angular momentum the thin black ring becomes identical to a highly boosted 
black string (see Section which is known to exhibit a Gregory-Laflamme instabil- 
ity. This suggests |15|. I116j that the thin black ring branch suffers from a Gregory- 
Laflamme instability quite possibly driving the system towards a configuration of 
two or more spinning black holes. There are some similarities here with the case of 
ultra-spinning Myers-Perry black holes in six or higher dimensions. 

• Black rings can become unstable under radial perturbations. The analysis of [116J 
indicates that fat black rings are radially unstable while thin black rings are radially 
stable. The black ring with minimal spin is also argued to be unstable. This type 
of instability appears to mesh nicely with the topological arguments of [218, 1219] 
which are based on the Poincare method. 

• The angular momentum of black rings is bounded from below. By absorbing a 
counter-spinning null-particle the black ring can decrease its angular momentum j. 
For the lowest-spinning black ring this process indicates an instability. On the other 
hand, it is impossible for a Myers- Perry black hole or a fat black ring to overspin 
by absorbing a co-rotating null particle or to spontaneously decay by emitting null 
particles without violating the area law [223J. 

• Doubly spinning large radius (i.e. v — > 0) thin black rings have been argued to suffer 
from superradiance instabilities [209J. 

The addition of charges can improve the stability of black rings. For instance, supersym- 
metric black rings |224j are expected to be stable. Black rings with dipole charges |115j 
are expected to suffer from similar instabilities as the neutral ones [223J. 
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8.2 Future directions and open problems 

We conclude with a number of important issues and questions for future research. 

Pure gravity: static solutions asymptoting to M d x S 1 

• At present the complete non-uniform branch is known numerically for d = 4 and 
d = 5 [301 H2] and recently also higher dimensions 6 < d < 10 have been studied 
[33] . It would be interesting to compute it above the critical dimension, i.e. d > 13. 
For the black hole branch, the present status is that we know numerically the entire 
branch only for d = 4 and d = 5 [44, 45, 46J. Having the numerical data for other 
dimensions would be interesting as well. 

• If possible, it would be interesting to find the analytic form of the non-uniform 
and black hole branches, for example using the ansatz (|3.1ip . A first step in this 
direction would be to extend the known first order analytical results to second order. 
For the d = 4 black hole this has been done in Ref. [H] , while [33] gives the corrected 
thermodynamics for all d, but not the metric. 

• For the neutral black string the study of the classical stability revealed the presence 
of the GL instability. It would be highly interesting if one could study the classical 
stability of other branches of solutions, such as the non-uniform black string, the 
localized black hole, or a black hole attached to a Kaluza-Klein bubble. 

• The bubble-black hole sequences are only known in five and six dimensions. It is 
interesting to ask if solutions for bubble-black hole sequences exist for D > 7, and if 
so, whether they can be related to the solutions of [50] by a map similar to the one 
mapping five- to six-dimensional solutions. For spaces with more than one Kaluza- 
Klein circle, i.e. asymptotics Ai d x T k - 2 , it is easy to construct solutions describing 
regular bubble-black hole sequences [50]. On the other hand, finding bubble-black 
hole sequences in A4 d x S 1 may be difficult since one cannot use the generalized 
Weyl ansatz for such solutions, but we do not see any physical obstructions to their 
existence. Should these higher dimensional bubble-black hole sequences exist, it 
would be interesting to see what type of horizon topologies they exhibit. 

• In connection with the above points, also the general question of the d-dependence 
of the phase diagram is an important one. For example, the non-uniform branch 
is known to have a critical dimension beyond which the slope changes sign. Like- 
wise one could imagine new phases appearing (or other phases disappearing) as the 
number of dimensions increases. 

• The instability of the uniform black string implies that a naked singularity may be 
formed when the horizon of the black string pinches off. As originally pointed out by 
Gregory and Laflamme [21 [3] , this would entail a violation of the Cosmic Censorship 
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Hypothesis. The results of [3U1 113 EES E3] suggest that, for certain dimensions, the 
localized black hole is the only solution with higher entropy than that of the uniform 
black string (for masses where the black string is unstable). Thus, the endpoint of 
the instability seems to be the localized black hole which means that the horizon 
should pinch off in the classical evolution. On the other hand, in [T7j it was argued 
that the horizon cannot pinch off in finite affine parameter, thus suggesting that 
it is impossible for the black string horizon to pinch off. A way to reconcile these 
two results is if the horizon pinches off in infinite affine parameter (see also |73j). 
Recently, the numerical analysis of [7U [76] indicates that this indeed is the case. 
If this is correct it would be interesting to examine the implications for the Cosmic 
Censorship Hypothesis. 

• A similar formation of a naked singularity (now for static configurations) has been 
predicted at the merger point of the non-uniform and localized black hole branch in 
the phase diagram of black strings. Direct evidence for the existence of the merger 
point has been found with a numerical study in [52] (see also [Ml [27]). It would be 
very interesting to explore the physics of the merger point transition further. 

Pure gravity: static solutions asymptoting to M d x J\f 

Static solutions in pure gravity asymptoting to A4 d x S 1 exhibit a rich structure. An even 
richer structure is expected for static solutions asymptoting to the more general manifold 
M d x M where J\f is any compact Ricci-flat manifold. It would be interesting to explore 
the N dependence of the phase structure of solutions in certain examples, e.g. with J\f 
being an n-torus T n or N being a Calabi-Yau manifold. Uniformly smeared solutions 
are still expected on general grounds to exhibit Gregory-Laflamme instabilities, but it 
would be further interesting to see if the phase diagram of the general case exhibits new 
qualitative features which are not present in the simplest N = S 1 case. A recent discussion 
of the Gregory-Laflamme instability for M = T n appeared in [37] . For some preliminary 
comments on the case of a general compact Ricci flat manifold N see Section 14.31 

Static black brane solutions in supergravity 

Many of the above pure gravity neutral static solutions can be easily lifted to corresponding 
solutions in type IIA or type IIB ten-dimensional supergravity where one can also make 
contact with D-brane configurations and gauge theories. It would be interesting to explore 
the full phase diagram of static solutions in ten-dimensional supergravity independently 
and see whether it admits new phases which cannot be obtained through a map from the 
neutral pure gravity ones. In this connection, we note that in Ref. [141] a new stable phase 
of non-uniform black strings was found. 
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Black hole microscopies 



Another interesting direction is to make contact with microscopic calculations of the en- 
tropy of black holes. While this is difficult for neutral Kaluza-Klein black holes, charging 
up the solutions by the boost/U-duality map and taking the near-extremal limit opens 
up the possibility of applying microscopic conjectures. This was successfully applied in 
Ref. [77] where five-dimensional three-charge black holes on a a circle were obtained by 
extending the map of Ref. [57] to include three charges. As a result a rich phase structure 
of three-charge black holes, including a new phase of three-charge black holes that are 
non-uniformly distributed on the circle, was found via a map that relates them to the 
phases of five-dimensional neutral Kaluza-Klein black holes. Moreover, for three-charge 
black holes localized on the transverse circle, it was shown that in a partial extremal limit 
with two charges sent to infinity and one finite, the first correction to the finite entropy 
is in agreement with the microscopic entropy. This was achieved by taking into account 
that the number of branes shift as a consequence of the interactions across the transverse 
circle. This matching was extended in Ref. [78] to second order. Furthermore, in Ref. [78] 
a simple microscopic model was proposed that reproduces most of the features of the 
phase diagram, including the new non-uniform phase. It would be interesting to pursue 
this microscopic description further (in this connection see e.g. the reviews \225\ [226J). 

Connection to thermodynamics 

A natural and intuitive connection between classical and thermodynamic instabilities of 
black brane solutions is expected to hold in general. The conjecture underlying this relation 
(Correlated Stability Conjecture) has passed a large number of tests, but is known to fail 
in certain cases as we saw in section [H It would be interesting to formulate a refinement 
of the conjecture that survives all known tests. Then one could hope to find a general 
proof for its validity. 

Other compactified solutions 

It would also be interesting to examine the existence of other classes of solutions with a 
compactified direction. For example in Ref. [227] a supersymmetric rotating black hole in a 
compactified space-time was found and charged black holes in compactified space-times are 
considered in Ref. [228] . In another direction, new solutions with Kaluza-Klein boundary 
conditions for anti-de-Sitter space times have recently been constructed in Refs. |229^ f230j. 

Black rings 

For black rings in five dimensions (and for that matter also for other higher dimensional 
black holes like the rotating Myers-Perry solutions) it is important to develop further the 
explicit linear stability analysis. The main obstacle in this enterprize is of course the 
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complexity of the equations in higher dimensions. The ultimate goal in that respect is a 
better understanding of the full phase diagram. Beyond that it would be very exciting 
to see if there are black ring solutions in D > 5 dimensions or other solutions with the 
more complex horizon topology allowed by topological arguments in higher dimensions 

[2m EE]. 

Braneworld black holes 

The higher dimensional black holes and branes described in this review appear naturally 
in the discussion of the braneworld model of large extra dimensions [531 154] . In other 
braneworld models such as the one proposed by Randall and Sundrum [2311 1232] the 
geometry is warped in the extra direction and the discovery of black hole solutions in this 
context has proven more difficult. A black string solution is also allowed in this model and 
is expected to have a Gregory-Laflamme instability as in the case of uniform black strings 
in KK space-times. Localized black hole solutions are anticipated as end-products of the 
Gregory-Laflamme instability of black strings, but have not been constructed explicitly 
so far. A partial list of references includes [233} 12341 12351 1236] (see also the recent work 
in [237, 238J and references therein). It should be noted that exact black hole solutions 
localized on the four-dimensional braneworld have been constructed by the authors of 
[23911220] . 

Acknowledgments 

We thank Ofer Aharony, Micha Berkooz, Jan de Boer, Poul Henrik Damgaard, Os- 
car Dias, Roberto Emparan, Dan Gorbonos, Shinji Hirano, Gary Horowitz, Veronika 
Hubeny, Barak Kol, Kristjan Kristjansson, Hideaki Kudoh, Finn Larsen, Luis Lehner, 
Ernesto Lozano-Tellechea, Don Marolf, Rob Myers, Mukund Rangamani, Simon Ross, 
Peter R0nne, Evgeny Sorkin and Toby Wiseman for useful discussions related to the ma- 
terial presented in this work. We thank Henriette Elvang for useful discussions and for 
collaboration on one of the topics presented in this work. We thank Hideaki Kudoh and 
Toby Wiseman for providing us with the numerical results of Ref. [46] , and Burkhard Klei- 
haus, Jutta Kunz and Eugen Radu for kindly providing the numerical results of Ref.[32j. 
We thank Ruth Gregory for permission to reprint Figure [TJ TH and NO would also like 
to thank the KITP for hospitality during the program "Scanning new horizons: GR be- 
yond 4 dimensions" , while part of this work was completed. The work of TH and NO is 
partially supported by the European Community's Human Potential Programme under 
contract MRTN-CT-2004-005104 'Constituents, fundamental forces and symmetries of the 
universe'. TH would like to thank the Carlsberg Foundation for support. VN acknowl- 
edges partial financial support by the EU under the contracts MEXT-CT-2003-509661, 
MRTN-CT-2004-005104 and MRTN-CT-2004-503369. 



100 



A The Gregory-Laflamme mode 



In this appendix we consider in more detail the Gregory-Laflamme mode [2, 3J for the met- 
ric (|2.9p that describes a uniform black string in D = d + 1 dimensions! 51 ! The instability 
mode is given by a metric perturbation h^ v that solves the perturbed Einstein equations 
of motion around the black string metric (12. 9|) . The precise form of the perturbation h^ u 
appears in the main text in eqs. (12. lip and (12.121) . In (12. 12H ifi, r/, x and k are all functions 
of the variable 



rk 
ro 



(A.l) 



and this will be the variable with respect to which we take derivatives. The function / is 
given in terms of x as 



/ = 1 



rj.d—3 



(A.2) 



Following |89 pi5H I35j we work here in a transverse (V^/i^ = 0), traceless (g^h^ = 0) 
gauge. The tracelessness condition is 



(d-2)K + x + ^ = . 



(A.3) 



The transversality conditions are 

n 



ij) + (d- 3)(1 - f)rj + / ((d - 2)r? + xrf) = , 



(A.4) 



-2-xr, + (d - 3)(1 - f)( X ~ 1>) + 2/ ((d - l) x + V + *x') = 



The four independent Einstein equations take the form 

k 

n 



x% + (d - 3)(1 - /)r? + £/ 2 ((d - 2)V/ + xV>") 



A: 
+/ 



2(d - 2)t? - x^V + 2x7?' + A(d - 3)(1 - /)(-%' + 



9. 



-(d-3)(l-/)x-( X + V) + 2/ 



x^((d - l)x + V + x (x + ip')) + ??x 2 



(A.5) 



0, (A.6) 



■'-I S + /)x-(d-3)(l-/)^77 



A: 2 



+/ 



-2x 2 ^t/ + (d - 3)(1 - f)x( X ' - tf) + xf{d X ' + 2i// + x X ") 



(A.7) 



, 



/ 



H - x 2 + 2(d - 3)(1 - /)) + X { ~ x 2 + 2(d - l)(d - 3)(1 - /)) 



. 2 (d - 2)^7? + (d - 3)(1 - /)x( X ' + 



+ r 



2(d-l)(d-3)x 



(A. 



+2(d - 3)V + x(3(d - 2)x' + (d - 2)^' + x(x" + <)) 



^ee also ^89 151, 35__ and the more recent [87] for related work on this. 
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Combining the gauge conditions ()A.3[) - (|A.4j) with the Einstein equations ()A.5[) - (|A.8|) one 
can derive Eq. (|2.13p which is a second order differential equation for ij) of the form 



<(ar) + Qd(x)if/(x) +Vd(x)ip(x) = . 



(A.9) 



Qd and are d-dependent rational functions of x, which we summarize here for com- 
pleteness 

Qd{x) = f^x' 1 \f 3 k 4 (39 + lOd + 3d 2 ) + n 2 (k 2 (3 - d) 2 - 4ft 2 ) 
+k 2 f 2 ( - 2fc 2 (41 - 16d + d 2 ) - 3{d + 5)(5 - 3d)fl 2 ) 
-k 2 f(k 2 (3 - df + 2(-5 - V2d + 5d 2 )0 2 )l 1 

x fc 4 / 4 (36 + lid + 6d 2 - d 3 ) - 3(3 - d)n 2 (k 2 (3 - d) 2 - 4ft 2 ) 
+k 2 f 2 {k 2 (3 - d)(64 - 23d + d 2 ) + ft 2 (-431 + 405d - 181d 2 + 31d 3 )) 
+k 2 f(k 2 (-178 + 28d - 2d 2 + Ad 3 ) + H 2 (135 - 94d + 51d 2 - 8d 3 )) 
f(2k 4 (3 - d) 3 + k 2 n 2 (3 - d)(163 - 127d + 26d 2 ) + 4(-7 + 2d)ft 4 )] (A.IO) 



3(5 + d)(5-3d)ft 2 

-l 

x 



2k & f{l-df{3-df 



V d {x) = k~ 2 f- 2 x- 2 f 3 k\39 + Wd + 3d 2 ) + n 2 (k 2 {3 - d) 2 - 4ft 2 ) 
+k 2 f 2 (-2k 2 (4l-l6d + d 2 ^ 1 JWc o^c^2\ 

-k 2 f(k 2 (3-d) 2 + 2(-5-l2d + 5d 2 )n 2 ) 
+k i f 4 (2k 2 (-3 + d)(-16 + 2ld - lid 2 + 2d 3 ) + (-39 - lOd - 3d 2 )x 2 ) 
-(3-d) 2 (l-d)(37-7d)ft 2 ) 

-2k A f 3 {k 2 {{-3 + d)(-25 + 21d - Id 2 + d 3 ) + (-41 + I6d - d 2 )x 2 ) 
+((3 - d)(-239 + 273d - 91d 2 + 9d 3 ) + (-21 + 18d + 7d 2 )x 2 )ft 2 ) 
+2fk 2 n 2 (k 2 (3 - d) 2 (35 - 14d + d 2 + x 2 ) 
+(2(3 - d)(d - 2) + (57 - 56d + 13d 2 )x 2 )ft 2 ) 
+ft 2 (/fc 4 (3 - d) 4 - A: 2 (3 - d) 2 (-2 + 7x 2 )ft 2 + 4x 2 ft 4 ) 

+k 2 f 2 (k A (3 - d) 2 (8 - 2d + x 2 ) - 2k 2 {(-3 + d)(-302 + 279d - 80d 2 + 7d 3 ) 
-6(d - 2) 2 x 2 )ft 2 + (2(d - 3)(d - 1) + (9 + 22d - 19d 2 )x 2 )ft 4 )] . (A.ll) 

The threshold unstable mode obeys the differential equation (]A.9P with $7 = 0, which 
is equivalent to the second order equation for the scalar \ that appears in (2.10) of [35] 

-2x 2 (ff"-f 2 )-x(d-2)ff + 2df- 



- fx" + 



x(xf> - 2f) 



X 



+ 



x 2 ff" + x 



2(d-l)ff" -(d + 2)f 2 +4//' 



x(xf - 2f) 



X 



-rlx 



(A.12) 



The differential equations governing the threshold unstable mode have appeared in the 
literature in various forms depending on the choice of gauge. A complete list of the gauge 
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choices that have been employed so far in the literature can be found in appendix A of 
[87j . Differential equations for the general time-dependent Gregory-Laflamme mode have 
appeared also in [2] (see eq. 10, which is a second order differential equation for h tr ) and 
in [56] (see eq. 12). 
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